'X1 x0 h0 vdp' presentation slideshows

function [x]=pointmilieu(a,t,f) N=length(t); n=length(a); x=zeros(n,N); x(:,1)=a; for i=2:N

function [x]=pointmilieu(a,t,f) N=length(t); n=length(a); x=zeros(n,N); x(:,1)=a; for i=2:N h=t(i)-t(i-1); x(:,i)=x(:,i-1)+h*f(t(i-1)+.5*h,x(:,i-1)+.5*h*f(t(i-1),x(:,i-1))); end endfunction. function tracemilieu(a,t) getf vdp.sci getf pointmilieu.sci x=pointmilieu(a,t,vdp);

168 views • 3 slides