数值计算程序.利用亚当姆斯方法求解微分方程y'=y-2x/y ; y(0)=1 的值.0≤x≤1

来源：学生作业帮编辑：作业帮分类：数学作业时间：2024/04/28 21:40:17

给你编写matlab的程序,行不?
%y'=y-2x/y ; y(0)=1,0¡Üx¡Ü1
%first you have to give the old message of y1\
%y2\y3.So you can use the Runge-Kutta method.
h=0.01;
x=0:h:1;
N=length(x);%The iterations
y=zeros(1,N);y(1)=1;
dy=zeros(1,N);
dy_1=zeros(1,N);
dy_2=zeros(1,N);
y_1=zeros(1,N);
y_2=zeros(1,N);
k1(1)=y(1)-2*x(1)/y(1);
k2(1)=y(1)+h*k1(1)/2-2*(x(1)+h/2)/(y(1)+h*k1(1)/2);
k3(1)=y(1)+h*k2(1)/2-2*(x(1)+h/2)/(y(1)+h*k2(1)/2);
k4(1)=y(1)+h*k3(1)-2*x(2)/(y(1)+h*k3(1));
for n=2:1:4
y(n)=y(n-1)+h*(k1(n-1)+2*k2(n-1)+2*(k3(n-1))+k4(n-1))/6;
k1(n)=y(n)-2*x(n)/y(n);
k2(n)=y(n)+h*k1(n)/2-2*(x(n)+h/2)/(y(n)+h*k1(n)/2);
k3(n)=y(n)+h*k2(n)/2-2*(x(n)+h/2)/(y(n)+h*k2(n)/2);
k4(n)=y(n)+h*k3(n)-2*x(n)/(y(n)+h*k3(n));
end
%the next is the main program (Adams method)
for t0=1:1:4
dy(t0)=y(t0)-2*x(t0)/y(t0);
end
for t=5:1:N
y_1(t)=y(t-1)+h*(55*dy(t-1)-59*dy(t-2)+37*dy(t-3)-9*dy(t-4))/24;
dy_1(t)=y_1(t)-2*x(t)/y_1(t);
y_2(t)=y(t-1)+h*(9*dy_1(t)+19*dy(t-1)-5*dy(t-2)+dy(t-3))/24;
dy_2(t)=y_2(t)-2*x(t)/y_2(t);
dy(t)=dy_2(t);
end
程序有些小问题,现在没时间改了.有时间再和你交流.

MATLAB求方程y’= y - 2x/y,y(0) = 1 (0≤x≤1,h = 0.1) 的数值解,要求编写程序求解微分方程(y^2-1)dx+(y^2-y+2x)dy=0 急求解微分方程y'=(x-y+1)^2, 求解微分方程（xy^2+x）dx+(y-x^2y)dy=0,y(2)=1的通解用matlab求微分方程 y’=-y2+x,y(0)=1的数值解, 求解微分方程x(y^2-1)dx+y(x^2-1)dy=0 dy/dx+y/x=x^2,y(1)=0求解微分方程高数中微分方程求解求微分方程y'cos^2x+y-tanx=0的通解计算微分方程 y'+y-e^(-x)=0的通解 MATLAB 求微分方程数值解 y''-y'+y=0 y(1)=2,y(2)=0,x∈[0,20] MATLAB代码求解微分方程y'cosx+ysinx=0 求解微分方程dy/dx=y/（x+y的平方）求解微分方程：(1) 2yy‘‘=(y‘)^2+y^2 (2) yy‘‘+(y‘)^2+2x=0