资 源 简 介
format long, format compact
[y,k]=quad("myerrf",-8,8)
[y,k]=quad("myerrf",-15,15)
f=inline("1/sqrt(2*pi)*exp(-x.^2/2)","x");
[y,k]=quad(f,-8,8)
[y,k]=quad8("myerrf",-8,8)
[y,k]=quad8(f,-8,8)
h_opt=odeset; x0=[-0.2;-0.7]; t_final=20;
mu=1; [t1,y1]=ode45("vdp_eq",[0,t_final],x0,h_opt,mu);
mu=2; [t2,y2]=ode45("vdp_eq",[0,t_final],x0,h_opt,mu);
clf,set(gcf,"position",[516 70 500 350])
plot(y1(:,1),y1(:,2),y2(:,1),y2(:,2))
