c12804d2183c486babd6e8531ef7821b8fe3b94f
[scilab.git] / scilab / modules / optimization / tests / nonreg_tests / bug_3161.dia.ref
1 // =============================================================================
2 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
3 // Copyright (C) 2008 - INRIA - Allan CORNET
4 //
5 //  This file is distributed under the same license as the Scilab package.
6 // =============================================================================
7 // <-- Non-regression test for bug 3161 -->
8 //
9 // <-- Bugzilla URL -->
10 // http://bugzilla.scilab.org/show_bug.cgi?id=3161
11 //
12 // <-- Short Description -->
13 // datafit crashes
14 function y=FF(x,p),y=p(1)*(x-p(2))+p(3)*x.*x,endfunction
15 X=[];Y=[];
16 pg=[34;12;14]; //parameter used to generate data
17 for x=0:.1:3, Y=[Y,FF(x,pg)+100*(rand()-.5)];X=[X,x];end
18 Z=[Y;X];
19 //The criterion function
20 function e=G(p,z)
21    y=z(1),x=z(2);
22    e=y-FF(x,p);
23 endfunction
24 //Solve the problem
25 p0=[3;5;10];
26 p = [];
27 err = [];
28 [p,err]=datafit(G,Z,p0);
29 if p == [] then bugmes();quit;end;
30 if err == [] then bugmes();quit;end;