Question: Newmark Method

Hi everybody,


I'm an mechanical enginer and I have to solve numerically this kind of differential equation : m(d2z/dt2)+c(dz/dt) + k*z = F where z is the displacement.
The best resolution for this problem is to use a Newmark method (which is available here at page 10,d.d2k)


So I have written this method in Maple but nothing happen and I don't where the mistakes are. So if somebody could help me, that will be great.

My code is :


Please Wait...