Advanced Search

Question: i have problem in dsolve numeric

Question:i have problem in dsolve numeric

kambiz1199 170 Maple
ode dsolve numeric differential-equations
February 13 2018
0

hi . i want to get d^2 z[1] and d^2z[2] but this method give only dz[1] and dz[2] ....

 

restart;
q1 := 9045*(diff(z[1](t), t$2))+863728.*z[1](t) = -1963*sin(20*t);
q2 := 3454*(diff(z[2](t), t$2))+2^6*z[2](t) = -286*sin(20*t);
icy := seq([z[i](0) = 0, (D(z[i]))(0) = 0], i = 1 .. 2, 1);
so := dsolve({q1, q2, seq(icy[i][], i = 1 .. 2)}, numeric, method = classical[abmoulton], corrections = 2);
so(1);
[t = 1., z[1](t) = 0.114730517257337e-2, diff(z[1](t), t) = 0.192255592847066e-1, z[2](t) = -0.186463122818374e-3, diff(z[2](t), t) = 0.320601041204689e-2]
Please Wait...