
function dy = vdp1000(t,y)
为规定的范式,其中y是代表函数的解,y(1)就是需要求的解y,y(2)是y的微分y’,相当于y(n)’=y(n+1)
dy = zeros(2,1);
作为函数返回值的列向量,里面存储的是解函数的微分,dy和y的区别在于相同索引的dy和y差一个微分d,相当于dy(1)=y(1)’=y(2)
dy(1)=y(2);
dy(1)存储的的是y’,可以简单记忆dy(1)=y’,dy(2)=y”,y(n)=dy(n-1)=y的n-1个’
dy(2)=1000*(1-y(1)^2)*y(2)-y(1);
注意到方恒的左边都是关于dy(n)的表达式,右边都是关于y(n)的表达式
end