Xcos solvers : Sundials Dormand-Price 4(5) implementation
[scilab.git] / scilab / modules / xcos / examples / solvers / integLSodar.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/integLSodar.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 = ["LSodar", "BDF/Newton", "BDF/Functional", "Adams/Newton", "Adams/Functional"];
19
20 for solver = 0:4
21
22  // Select the solver
23  scs_m.props.tol(6) = solver;
24
25  // Start the timer, launch the simulation and display time
26  tic();
27  try scicos_simulate(scs_m, 'nw'); catch disp(lasterror()); end;
28  t = toc();
29  disp(t, "Time for " + solverName(solver+1) + " :");
30
31 end