fvdb = @(t,y) [y(2); (1 - y(1)^2) * y(2) - y(1)]; opt = odeset ("OutputFcn", @odeplot, "RelTol", 1e-6); tic; sol = ode45 (fvdb, [0 10], [2; 0], opt); toc