# HG changeset patch # User Nicholas R. Jankowski # Date 1714715722 14400 # Node ID b5fb87bf2008dab127265b3e077249ec683337c6 # Parent 09dd34fe53be7940e9ab206a36c79872b855c408# Parent 8e555ddce8001c6ce82da7fb6d29d9b035f59837 maint: Merge stable to default. diff -r 09dd34fe53be -r b5fb87bf2008 scripts/ode/ode23s.m --- a/scripts/ode/ode23s.m Fri May 03 01:39:49 2024 -0400 +++ b/scripts/ode/ode23s.m Fri May 03 01:55:22 2024 -0400 @@ -277,9 +277,9 @@ %! ## Demo function: stiff Van Der Pol equation %! fcn = @(t,y) [y(2); 10*(1-y(1)^2)*y(2)-y(1)]; %! ## Calling ode23s method -%! tmr = tic (); +%! tic () %! [vt, vy] = ode23s (fcn, [0 20], [2 0]); -%! toc (tmr); +%! toc () %! ## Plotting the result %! plot (vt,vy(:,1),'-o'); @@ -289,9 +289,9 @@ %! ## Calling ode23s method %! odeopts = odeset ("Jacobian", @(t,y) [0 1; -20*y(1)*y(2)-1, 10*(1-y(1)^2)], %! "InitialStep", 1e-3) -%! tmr = tic (); +%! tic () %! [vt, vy] = ode23s (fcn, [0 20], [2 0], odeopts); -%! toc (tmr); +%! toc () %! ## Plotting the result %! plot (vt,vy(:,1),'-o'); @@ -300,9 +300,9 @@ %! fcn = @(t,y) [y(2); 100*(1-y(1)^2)*y(2)-y(1)]; %! ## Calling ode23s method %! odeopts = odeset ("InitialStep", 1e-4); -%! tmr = tic (); +%! tic () %! [vt, vy] = ode23s (fcn, [0 200], [2 0]); -%! toc (tmr); +%! toc () %! ## Plotting the result %! plot (vt,vy(:,1),'-o'); @@ -312,9 +312,9 @@ %! ## Calling ode23s method %! odeopts = odeset ("Jacobian", @(t,y) [0 1; -200*y(1)*y(2)-1, 100*(1-y(1)^2)], %! "InitialStep", 1e-4); -%! tmr = tic (); +%! tic () %! [vt, vy] = ode23s (fcn, [0 200], [2 0], odeopts); -%! toc (tmr); +%! toc () %! ## Plotting the result %! plot (vt,vy(:,1),'-o');