69d4c17ee81c5da52b0bf44fe9982d238fe17d08
[scilab.git] / scilab / modules / optimization / tests / nonreg_tests / bug_2330.tst
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
8 // <-- Non-regression test for bug 2330 -->
9 //
10 // <-- Bugzilla URL -->
11 // http://www.scilab.org/cgi-bin/bugzilla_bug_II/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
18 //build the data to fit
19 //---------------------
20 function Xcalc=biexp(p,t)
21   Xcalc=p(1).*exp(-p(2).*t)+p(3).*exp(-p(4).*t)+p(5);
22 endfunction;
23 t=(0:100:36000)'; 
24 p=[0.1;0.0001;0.2;0.0002;0.3];
25 X=biexp(p,t);
26
27 //try to fit the data
28 //-------------------
29 //the error function
30 function e=myerf(p,X,t),e=X-biexp(p,t),endfunction
31 // the initial point
32 p0=[0.01;0.001;0.01;0.001;0.1];
33 //call datafit
34 [pr,err]=datafit(list(myerf,t),X,p0);
35 if err>=5d-6 then pause,end