一个mathematica程序添加作图语句Clear[x,y,n,h,S1,S2,S3,S4,i];y[1] = 0; h = 0.1;n = 100Do[x[i] = 1.0 + (i - 1)*h,{i,1,n + 1}]f[x_,y_] = 2.5-0.8 y;Do[{S1 = f[x[i],y[i]],S2 = f[x[i] + h/2,y[i] + h*S1/2],S3 = f[x[i] + h/2,y[i] + h*S2/2],S4 = f[x
Clear[x,y,n,h,S1,S2,S3,S4,i];
y[1] = 0; h = 0.1;
n = 100
Do[x[i] = 1.0 + (i - 1)*h,{i,1,n + 1}]
f[x_,y_] = 2.5-0.8 y;
Do[{S1 = f[x[i],y[i]],S2 = f[x[i] + h/2,y[i] + h*S1/2],
S3 = f[x[i] + h/2,y[i] + h*S2/2],S4 = f[x[i] + h,y[i] + h*S3],
y[i + 1] = y[i] + (S1 + 2 S2 + 2 S3 + S4)*h/6},{i,1,n}]
Do[Print[x[i]," ",y[i]],{i,1,n + 1}]
需要在后面做出图形,求问需要添加的语句.
你要是想画出 x 与 y 的图像可以用 ListPlot 函数,需要添加的语句如下:
ListPlot[Table[{x[i], y[i]}, {i, n + 1}]]输出图像:

要想将点连接起来可以加上选项 Joined .
更多图像选项及细节调整,请参考函数 ListPlot 的用法及例子.
函数y=x-根号下(1-2x)的值域是?令根号下(1-2x)=t则t≥0,x=-(t²-1)/2函 2020-04-26 …
设{X(t),t>=0}是正交增量过程,X(0)=0,V是标准正态随机变量,若对任意的t>=0,X 2020-05-13 …
求过曲线y=∫(0,x)(t-1)(t-2)dt上点(0,0)处切线方程其中0是下限,x是上限, 2020-07-22 …
已知:dx(t)/dt=y(t)x(t),初始条件为x(0)=0,y(t)不恒等于0.x(t)必定 2020-07-23 …
∫(t:0->x)t^2*e^(-t^2)d(t^2)求该定积分,X为一常数,希望有每步详细步骤 2020-07-23 …
变上限积分求导f(x)=∫(0,x)(x-t)^2*f(t)dt如何对x求导?原题:函数f(x)满 2020-07-31 …
这道高数的上限积分求到怎么求,说明白点!F(x,0)(x-t)f‘(x)dt的导数是多少!F(x, 2020-07-31 …
∫(0→π)√(sin∧3x-sin∧5x)dx正确答案是4/5,如果令sinx=t,x=0→t= 2020-08-01 …
mathematica求解这个微分方程DSolve[{x[t]x''[t]+x'[t]^2-gx[ 2020-08-02 …
周期函数的定义限定∀x∈D,∃T≠0,x±T∈D吗?周期函数的一个必要条件是:∀x∈D,∃T≠0,f 2021-01-04 …