2.3.m

利用龙格库塔法（等步长

clear all
t0=0;tend=2*pi;
u0=[0]';


%[time,result]=ode23(@fvdp1,[t0,tend],u0);
[time,result]=ode23(@fvd,[t0,tend],u0);
N=length(result);
A=abs(result);
for j=1:N   
   iv(j)=1*10^(-14)*(exp((result(j)-5/sqrt(2))/0.0259)-1);
   B(j)=abs(iv(j));
end
A=A-5/sqrt(2);
plot(A,B,'*');