Итогом должен быть график, зависимости k от x1(1)
часть программы:
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;
end proc: