方程解y的曲线一阶导数
clc,clear,close all
H=[0,20];
z0=[0; 1];
[x,z]=ode45('dzdx3',H,z0);
plot(x,z(:,1),'b--',x,z(:,2),'r--')
xlabel('轴\\it x'); ylabel('轴\\it y')
legend('是方程解y的曲线','是解y的一阶导数')
function dz=dzdx3(x,z)
dz(1)=z(2);
dz(2)=5*(1- 2*z(1)^4)* z(2)-7*z(1);
dz=[dz(1);dz(2)];
end
因篇幅问题不能全部显示，请点此查看更多更全内容