edsODE(3,[0; 1; 2]) help ode45 [t,y] = ode45(@edsODE,[1 3],[0; 0; 2]); [t,x] = ode45(@edsODE,[1 3],[0; 0; 2]); clear y x(end,1) format long x(end,1) plot(t,x) legend('y','v','w') figure plot(t,x(:,1)) help slash x 53/(log(10)/log(2))