[differential_equations] fix tests
[scilab.git] / scilab / modules / differential_equations / tests / nonreg_tests / bug_15523.tst
1 // =============================================================================
2 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3 // Copyright (C) 2018 - St├ęphane MOTTELET
4 //
5 //  This file is distributed under the same license as the Scilab package.
6 // =============================================================================
7
8 // <-- CLI SHELL MODE -->
9 // <-- NO CHECK REF -->
10
11 // <-- Non-regression test for bug 15523 -->
12 //
13 // <-- Bugzilla URL -->
14 // http://bugzilla.scilab.org/show_bug.cgi?id=15523
15 //
16 // <-- Short Description -->
17 //    %ODEOPTIONS(1)=2 didn't work with solvers 'rk' and 'rkf'
18 // =============================================================================
19
20 function ydot=f(t, y)
21     ydot=y^2-y*sin(t)+cos(t)
22 endfunction
23 %ODEOPTIONS=[2,0,0,%inf,0,2,500,12,5,0,-1,-1];
24 rtol=1e-5;
25 atol=1e-7;
26 yrkf=ode("rkf",0,0,%pi,rtol,atol,f);
27 yrk=ode("rk",0,0,%pi,rtol,atol,f);
28 assert_checktrue(size(yrkf,2)>1);
29 assert_checktrue(size(yrk,2)>1);
30 assert_checkalmostequal(yrk(2,:),sin(yrk(1,:)),1e-4,1e-4,"element")
31 assert_checkalmostequal(yrkf(2,:),sin(yrkf(1,:)),1e-4,1e-4,"element")