lemon-project-template-glpk
annotate deps/glpk/examples/cf12a.mod @ 10:5545663ca997
Configure GLPK build
| author | Alpar Juttner <alpar@cs.elte.hu> |
|---|---|
| date | Sun, 06 Nov 2011 21:42:23 +0100 (2011-11-06) |
| parents | |
| children |
| rev | line source |
|---|---|
| alpar@9 | 1 /* |
| alpar@9 | 2 |
| alpar@9 | 3 Curve fitting problem 12.11(a) H P Williams "Model Building in Mathematical Programming" |
| alpar@9 | 4 |
| alpar@9 | 5 Dr. H J Mackenzie |
| alpar@9 | 6 HARD software |
| alpar@9 | 7 hjm@hardsoftware.com |
| alpar@9 | 8 |
| alpar@9 | 9 2006-01-05 |
| alpar@9 | 10 |
| alpar@9 | 11 */ |
| alpar@9 | 12 |
| alpar@9 | 13 # set of points |
| alpar@9 | 14 |
| alpar@9 | 15 set I; |
| alpar@9 | 16 |
| alpar@9 | 17 # independent variable |
| alpar@9 | 18 |
| alpar@9 | 19 param x {i in I}; |
| alpar@9 | 20 |
| alpar@9 | 21 # dependent variable |
| alpar@9 | 22 |
| alpar@9 | 23 param y {i in I}; |
| alpar@9 | 24 |
| alpar@9 | 25 # output input values |
| alpar@9 | 26 |
| alpar@9 | 27 printf {i in I} "x = %.1f; y = %.1f\n", x[i], y[i]; |
| alpar@9 | 28 |
| alpar@9 | 29 # define equation variables |
| alpar@9 | 30 |
| alpar@9 | 31 var a; |
| alpar@9 | 32 |
| alpar@9 | 33 var b; |
| alpar@9 | 34 |
| alpar@9 | 35 var u {i in I}, >= 0; |
| alpar@9 | 36 |
| alpar@9 | 37 var v {i in I}, >= 0; |
| alpar@9 | 38 |
| alpar@9 | 39 # define objective function |
| alpar@9 | 40 |
| alpar@9 | 41 minimize error: sum {i in I} u[i] + sum {i in I} v[i]; |
| alpar@9 | 42 |
| alpar@9 | 43 # define equation constraint |
| alpar@9 | 44 |
| alpar@9 | 45 s.t. equation {i in I} : b * x[i] + a + u[i] - v[i] = y[i]; |
| alpar@9 | 46 |
| alpar@9 | 47 solve; |
| alpar@9 | 48 |
| alpar@9 | 49 printf "y = %.4fx + %.4f\n", b, a; |
| alpar@9 | 50 |
| alpar@9 | 51 /* |
| alpar@9 | 52 * |
| alpar@9 | 53 * DATA section |
| alpar@9 | 54 * |
| alpar@9 | 55 */ |
| alpar@9 | 56 |
| alpar@9 | 57 data; |
| alpar@9 | 58 |
| alpar@9 | 59 param : I : x y := |
| alpar@9 | 60 1 0 1 |
| alpar@9 | 61 2 0.5 0.9 |
| alpar@9 | 62 3 1 0.7 |
| alpar@9 | 63 4 1.5 1.5 |
| alpar@9 | 64 5 1.9 2 |
| alpar@9 | 65 6 2.5 2.4 |
| alpar@9 | 66 7 3 3.2 |
| alpar@9 | 67 8 3.5 2 |
| alpar@9 | 68 9 4 2.7 |
| alpar@9 | 69 10 4.5 3.5 |
| alpar@9 | 70 11 5 1 |
| alpar@9 | 71 12 5.5 4 |
| alpar@9 | 72 13 6 3.6 |
| alpar@9 | 73 14 6.6 2.7 |
| alpar@9 | 74 15 7 5.7 |
| alpar@9 | 75 16 7.6 4.6 |
| alpar@9 | 76 17 8.5 6 |
| alpar@9 | 77 18 9 6.8 |
| alpar@9 | 78 19 10 7.3 |
| alpar@9 | 79 ; |
| alpar@9 | 80 |
| alpar@9 | 81 end; |