b74edf0b63dcc10cdc1d1638ff2b153089c27d43
[scilab.git] / scilab / modules / optimization / tests / nonreg_tests / bug_2330.dia.ref
1 // =============================================================================
2 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3 // Copyright (C) 2007-2008 - INRIA - Serge STEER <serge.steer@inria.fr>
4 //
5 //  This file is distributed under the same license as the Scilab package.
6 // =============================================================================
7 // <-- CLI SHELL MODE -->
8 // <-- Non-regression test for bug 2330 -->
9 //
10 // <-- Bugzilla URL -->
11 // http://bugzilla.scilab.org/show_bug.cgi?id=2330
12 //
13 // <-- Short Description -->
14 //    datafit does not work equivalently on WinXP and Linux. With a given dataset and
15 //    same routines centered on datafit function, it works perfectly on Linux and
16 //    partially (some data are fitted some others not) on WinXP.
17 //build the data to fit
18 //---------------------
19 function Xcalc=biexp(p,t)
20   Xcalc=p(1).*exp(-p(2).*t)+p(3).*exp(-p(4).*t)+p(5);
21 endfunction;
22 t=(0:100:36000)';
23 p=[0.1;0.0001;0.2;0.0002;0.3];
24 X=biexp(p,t);
25 //try to fit the data
26 //-------------------
27 //the error function
28 function e=myerf(p,X,t),e=X-biexp(p,t),endfunction
29 // the initial point
30 p0=[0.01;0.001;0.01;0.001;0.1];
31 //call datafit
32 [pr,err]=datafit(list(myerf,t),X,p0);
33 if err>=5d-4 then bugmes();quit;end