lemon-project-template-glpk: deps/glpk/examples/cflsq.mod annotate

lemon-project-template-glpk

Import GLPK 4.47

author	Alpar Juttner <alpar@cs.elte.hu>
date	Sun, 06 Nov 2011 20:59:10 +0100
parents
children

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