clear all
clc
f=@(x,y)([y(2);2*x*y(2)/(1+x^2)]);
[X,Y]=ode45(f,0:100,[1 3]);
y=dsolve('D2y-2*x*Dy/(1+x^2)','y(0)=1','Dy(0)=3','x');
y=@(x)(eval(y));
YY=y(0:100);
plot(0:100,Y(:,1),'o',0:100,YY)
legend('数值解','解析解')
内容全部来源于网络收集，如有侵权，请联系网站删除：QQ：24596024