如题所述
求解微分方程的函数是dsolve()。如要画微分方程的解函数,可以用plot()函数。例如:
求解dy/dx=1+y^2,并且当y(0)=1的解,画出y=f(x)的图形。
y=dsolve('Dy=1+y^2','y(0)=1','x')
x=-2*pi:pi/10:2*pi;
y1=eval(y);
plot(x,y1,'k-')
xlabel('x'),ylabel('y(x)');grid on