Merge remote-tracking branch 'origin/6.1'
[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 is released under the 3-clause BSD license. See COPYING-BSD.
5
6 // Import the diagram and augment the ending time
7 loadScicos();
8 loadXcosLibs();
9 importXcosDiagram("SCI/modules/xcos/examples/solvers/ODE_Example.zcos");
10 scs_m.props.tf = 30000;
11
12 solverName=["BDF/Newton", "BDF/Functional", "Adams/Newton", "Adams/Functional", "Dormand-Prince"];
13
14 for solver=1:5
15
16     // Select the solver
17     scs_m.props.tol(6) = solver;
18
19     // Set max step size if DoPri
20     if (solver == 5) then scs_m.props.tol(7) = 0.01; end
21
22     // Start the timer, launch the simulation and display time
23     tic();
24     try scicos_simulate(scs_m, "nw"); catch disp(lasterror()); end;
25     t = toc();
26     disp("Time for " + solverName(solver) + ":",t);
27
28 end