lemon-project-template-glpk

comparison deps/glpk/examples/yacfs.mod @ 9:33de93886c88

find changesets by author, revision, files, or words in the commit message
Import GLPK 4.47
author Alpar Juttner <alpar@cs.elte.hu>
date Sun, 06 Nov 2011 20:59:10 +0100
parents
children
comparison
equal deleted inserted replaced
-1:000000000000 0:dbfd47dd9bdd
1 /*Yet Another Curve Fitting Solution
2
3 Obviously this solution produces the same answer
4 as examples/cflsq.mod
5
6 Nigel_Galloway@operamail.com
7 February 1st., 2009
8 */
9 set Sample;
10 param Sx {z in Sample};
11 param Sy {z in Sample};
12
13 var a;
14 var b;
15
16 equalz1 :sum{z in Sample} a*Sx[z]*Sx[z] + sum{z in Sample} b*Sx[z] = sum{z in Sample} Sy[z]*Sx[z];
17 equalz2 :sum{z in Sample} a*Sx[z] + sum{z in Sample} b = sum{z in Sample} Sy[z];
18
19 solve;
20
21 printf "\nbest linear fit is:\n\ty = %f %s %fx\n\n", b, if a < 0 then "-" else "+", abs(a);
22
23 data;
24
25 param:
26 Sample: Sx Sy :=
27 1 0 1
28 2 0.5 0.9
29 3 1 0.7
30 4 1.5 1.5
31 5 1.9 2
32 6 2.5 2.4
33 7 3 3.2
34 8 3.5 2
35 9 4 2.7
36 10 4.5 3.5
37 11 5 1
38 12 5.5 4
39 13 6 3.6
40 14 6.6 2.7
41 15 7 5.7
42 16 7.6 4.6
43 17 8.5 6
44 18 9 6.8
45 19 10 7.3
46 ;
47
48 end;