differential_coefficient.m

自己编的微分方程求解的程序,没有调用matlab的函数

clear;
close all;
dt=0.01;
u=5;
y1(1)=1;
y2(1)=-1;
for i=1:1000
    y1(i+1)=y2(i)*dt+y1(i);
    y2(i+1)=[u*(1-y1(i)^2)*y2(i)-y1(i)]*dt+y2(i);
end
i=1:1000;
hfig1=figure;
set(hfig1,'color',[1 1 0.5],'position',[0.2*1024 0.2*768 0.6*1024 0.6*768],'name','Differential Coefficient','menubar','none','numbertitle','off');
plot(i,y1(i))
title('Programming by myself','color','r','fontsize',15,'fontname','Time New Roman');