COIN-OR::LEMON - Graph Library

source: glpk-cmake/examples/cflsq.mod @ 1:c445c931472f

Last change on this file since 1:c445c931472f was 1:c445c931472f, checked in by Alpar Juttner <alpar@…>, 13 years ago

Import glpk-4.45

  • Generated files and doc/notes are removed
File size: 1.0 KB
RevLine 
[1]1/*Curve fitting problem by Least Squares
2  Nigel_Galloway@operamail.com
3  October 1st., 2007
4*/
5set Sample;
6param Sx {z in Sample};
7param Sy {z in Sample};
8
9var X;
10var Y;
11var Ex{z in Sample};
12var Ey{z in Sample};
13
14/* sum of variances is zero for Sx*/
15variencesX{z in Sample}: X + Ex[z] = Sx[z];
16zumVariancesX: sum{z in Sample} Ex[z] = 0;
17/* sum of variances is zero for Sy*/
18variencesY{z in Sample}: Y + Ey[z] = Sy[z];
19zumVariancesY: sum{z in Sample} Ey[z] = 0;
20
21solve;
22
23param b1 := (sum{z in Sample} Ex[z]*Ey[z])/(sum{z in Sample} Ex[z]*Ex[z]);
24printf "\nbest linear fit is:\n\ty = %f %s %fx\n\n", Y-b1*X, if b1 < 0 then "-" else "+", abs(b1);
25
26data;
27
28param:
29Sample:   Sx    Sy :=
30  1         0    1
31  2       0.5  0.9
32  3         1  0.7
33  4       1.5  1.5
34  5       1.9    2
35  6       2.5  2.4
36  7         3  3.2
37  8       3.5    2
38  9         4  2.7
39 10       4.5  3.5
40 11         5    1
41 12       5.5    4
42 13         6  3.6
43 14       6.6  2.7
44 15         7  5.7
45 16       7.6  4.6
46 17       8.5    6
47 18         9  6.8
48 19        10  7.3
49;
50
51end;
Note: See TracBrowser for help on using the repository browser.