a75698b4c8a4ca8e1dc58bade7b6909f05656142
[scilab.git] / scilab / modules / xcos / examples / solvers / integDoPri.sce
1 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
2 // Copyright (C) 2012 - Scilab Enterprises - Paul Bignier
3 //
4 // This file must be used under the terms of the CeCILL.
5 // This source file is licensed as described in the file COPYING,
6 // which you should have received as part of this distribution.
7 // The terms are also available at
8 // http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
9
10 // Run with exec("SCI/modules/xcos/examples/solvers/integDoPri.sce");
11
12 // Import the diagram and augment the ending time
13 loadScicos();
14 loadXcosLibs();
15 importXcosDiagram("SCI/modules/xcos/examples/solvers/ODE_Example.xcos");
16 scs_m.props.tf = 30000;
17
18 solverName=["BDF/Newton", "BDF/Functional", "Adams/Newton", "Adams/Functional", "Dormand-Price"];
19
20 for solver=1:5
21
22  // Select the solver
23  scs_m.props.tol(6) = solver;
24
25  // Set max step size if DoPri
26  if (solver == 5) then scs_m.props.tol(7) = 0.01; end
27
28  // Start the timer, launch the simulation and display time
29  tic();
30  try scicos_simulate(scs_m, 'nw'); catch disp(lasterror()); end;
31  t = toc();
32  disp(t, "Time for " + solverName(solver) + " :");
33
34 end