Вывести график зависимости x(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:
>
> G := proc (x) options operator, arrow; sqrt(x+1) end proc;
x -> sqrt(x + 1)
> k := 1;
1
c(k);
1.03829835612414167 A1 + 1.00726621369385327 A2 = 1,
0.150470140010300679 A1 + 1.03550157382713648 A2 = 0
> solve({1.03829835612414167*A1+1.00726621369385327*A2 = 1, .150470140010300679*A1+1.03550157382713648*A2 = 0});
{A1 = 1.121162791, A2 = -0.1629176878