решение ДУ в maple

#25690

Решение задач

Не указан

Механика

Не указана

пост роить график решений х1(1) от k. начало программы: c:=proc(k) local F1,F2,func,y,X,cond1,cond2,cond; func := {x1(w), x2(w), x3(w), x4(w)}; y := diff(x1(w), w) = x2(w), diff(x2(w), w) = x3(w)/G(w), diff(x3(w), w) = x4(w), diff(x4(w), w) = k*x1(w); cond1 := x1(0) = 0, x2(0) = 0, x3(0) = 1, x4(0) = 0; cond2 := x1(0) = 0, x2(0) = 0, x3(0) = 0, x4(0) = 1; cond:=x1(0) = 0, x2(0) = 0, x3(1) = 1, x4(1) = 0; X := dsolve({cond, y}, func, numeric); F1 := dsolve({cond1, y}, func, numeric); F2 := dsolve({cond2, y}, func, numeric); #A1:=solve(F1) X(1),F1(1),F2(1); rhs(F1(1)[4])*A1+rhs(F2(1)[4])*A2=1, rhs(F1(1)[5])*A1+rhs(F2(1)[5])*A2=0; #A1*rhs(F1(1))+A2*rhs(F2(1)); end proc: