Xcos solvers : Sundials Dormand-Price 4(5) implementation
[scilab.git] / scilab / modules / xcos / examples / solvers / integRK.sce
index a8e93f3..52822aa 100644 (file)
@@ -1,10 +1,10 @@
 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-// Copyright (C) 2009 - Paul Bignier
+// Copyright (C) 2012 - Scilab Enterprises - Paul Bignier
 //
 // This file must be used under the terms of the CeCILL.
-// This source file is licensed as described in the file COPYING, which
-// you should have received as part of this distribution.  The terms
-// are also available at
+// This source file is licensed as described in the file COPYING,
+// which you should have received as part of this distribution.
+// The terms are also available at
 // http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
 
 // Run with exec("SCI/modules/xcos/examples/solvers/integRK.sce");
@@ -17,18 +17,18 @@ scs_m.props.tf = 30000;
 
 solverName=["BDF/Newton", "BDF/Functional", "Adams/Newton", "Adams/Functional", "Runge-Kutta"];
 
-for solver=0:4
+for solver=1:5
 
  // Select the solver
- scs_m.props.tol(6) = solver+1;
+ scs_m.props.tol(6) = solver;
 
  // Set max step size if Runge-Kutta
- if ((solver+1) == 5) scs_m.props.tol(7) = 0.01;
+ if (solver == 5) scs_m.props.tol(7) = 0.01;
 
  // Start the timer, launch the simulation and display time
  tic();
  try scicos_simulate(scs_m, 'nw'); catch disp(lasterror()); end;
  t = toc();
- disp(t, "Time for " + solverName(solver+1) + " :");
+ disp(t, "Time for " + solverName(solver) + " :");
 
 end