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