@acer 

yeah got it. Thanks a lot once again.

@acer 

yes it is copy(x) and  the mutate procedure need to return y. 

@mmcdara 

I am getting the same answer. So what is the problem in writng the peicewisae function as i did.

@Carl Love 

Thank you

@mmcdara 

I am using maple 13. Is that the reason for this difference

@Preben Alsholm 

Why only first expression peicewise function on rhs (i.e for t<>8*n) only considered avoiding t=8*n

@sand15 

Is this correct. please reply

problem30.mw

@sand15 

firstl i assume n as positive integer and solve the system of differential equation for n. I repeat we need solution for general n not for a specific n.

problem30.mw

@Kitonum 

in your case you are only considering n=1. when t=8, n=1and t=16 then n=2 and so on. My main concern is how to deal with this while plotting against t.

@sand15 

as you can see in the problem t=n*8 and t not equal to n*8 where n can take 1,2,3.... So when you plot each against t we cannot set n a specific value. As i said it can take 1,2,3.... i.e t=8,t=16,t=24 and so on. The rhs of each diffrential equation is peicewise so at t=8*n and otherwise you will have different solution. 

@vv 

But when i run your code i get diffrent result. I am using maple  13.

Download K_(2).mw

@Mariusz Iwaniuk 

Thanks. But I don't want to put value for p. It need to be p as such.

@Mariusz Iwaniuk 

Thank you. I dont understand why others confirmed that it won't converge.

@Mariusz Iwaniuk 

I am waiting for your reply

@vv

Please see below code what is the difference between yours and mine. Thank you.

p7.mw