枫树中非线性DE的泰勒级数与数值解
Taylor series vs numeric solution of nonlinear DE in maple
我想绘制两个图:DE 的数值解和给定的 DE 的泰勒级数逼近。我有
de := diff(y(x), x) = x+y(x)-y(x)^2;
cond := y(0) = -1, (D(y))(0) = 1;
stp := 0.1e-1;
a, b := -5, 30;
numpts := floor((b-a)/stp+1);
p := dsolve({cond, de}, y(x), numeric, stepsize = stp, output = listprocedure);
绘图 eval
给出了奇怪的垂直线,而我希望获得的绘图似乎像 x -> ∞
一样振荡。对于泰勒级数,我试过 f:=[seq(taylor(y(x),x=i,n),i=-5..30 by stp)];
但似乎不会以这种方式工作。我能用它做什么?为什么我的剧情和预期的不一样?
restart;
kernelopts(version);
Maple 2018.0, X86 64 LINUX, Mar 9 2018, Build ID 1298750
de := diff(y(x), x) = x+y(x)-y(x)^2:
cond := y(0) = -1, (D(y))(0) = 1:
stp := 0.1e-1:
a, b := -5, 30:
numpts := floor((b-a)/stp+1):
p := dsolve({cond, de}, y(x), numeric, stepsize = stp,
output = listprocedure):
Y:=eval(y(x),p);
Y := proc(x) ... end;
plot(Y, 0..20);
Order:=10:
S := convert(rhs(dsolve({cond, de}, {y(x)}, series)),polynom);
plot([S, Y(x)], x=0..1.5);
Order:=40:
S := convert(rhs(dsolve({cond, de}, {y(x)}, series)),polynom):
plot([S, Y(x)], x=0..2.0);
我想绘制两个图:DE 的数值解和给定的 DE 的泰勒级数逼近。我有
de := diff(y(x), x) = x+y(x)-y(x)^2;
cond := y(0) = -1, (D(y))(0) = 1;
stp := 0.1e-1;
a, b := -5, 30;
numpts := floor((b-a)/stp+1);
p := dsolve({cond, de}, y(x), numeric, stepsize = stp, output = listprocedure);
绘图 eval
给出了奇怪的垂直线,而我希望获得的绘图似乎像 x -> ∞
一样振荡。对于泰勒级数,我试过 f:=[seq(taylor(y(x),x=i,n),i=-5..30 by stp)];
但似乎不会以这种方式工作。我能用它做什么?为什么我的剧情和预期的不一样?
restart;
kernelopts(version);
Maple 2018.0, X86 64 LINUX, Mar 9 2018, Build ID 1298750
de := diff(y(x), x) = x+y(x)-y(x)^2:
cond := y(0) = -1, (D(y))(0) = 1:
stp := 0.1e-1:
a, b := -5, 30:
numpts := floor((b-a)/stp+1):
p := dsolve({cond, de}, y(x), numeric, stepsize = stp,
output = listprocedure):
Y:=eval(y(x),p);
Y := proc(x) ... end;
plot(Y, 0..20);
Order:=10:
S := convert(rhs(dsolve({cond, de}, {y(x)}, series)),polynom);
plot([S, Y(x)], x=0..1.5);
Order:=40:
S := convert(rhs(dsolve({cond, de}, {y(x)}, series)),polynom):
plot([S, Y(x)], x=0..2.0);