A numerical scheme for initial compliance and creep response of a system