如题所述
该常微分方程的数值解可以使用ode45函数求解得到。
>> fun=@(x,y)[y(2);-4/x-2*y(1)/x^2];
>> [x,y]=ode45(fun,[1,10],[2,-3])
>> plot(x,y)
>>xlabel('x');ylabel('y(x),y‘(x)');