Anonymous user
Polynomial regression: Difference between revisions
m
→[[Polynomial_Fitting#ALGOL 68]]: comments
m (→[[Polynomial_Fitting#ALGOL 68]]: comments) |
|||
Line 60:
<lang algol>MODE FIELD = REAL;
MODE
VEC = [0]FIELD,
MAT = [0,0]FIELD;
Line 73:
[2 LWB in:2 UPB in,1 LWB in:1UPB in]FIELD out;
FOR i FROM LWB in TO UPB in DO
out[,i]:=in[i,]
OD;
out
);
COMMENT from http://rosettacode.org/wiki/Matrix_multiplication#ALGOL_68 END COMMENT
OP * = (VEC a,b)FIELD: ( # basically the dot product #
FIELD result:=0;
Line 112:
[,]FIELD lu = lu decomp(b, p, sign);
[LWB a:UPB a, 2 LWB a:2 UPB a]FIELD out;
FOR col FROM 2 LWB a TO 2 UPB a DO
out[,col] := lu solve(b, lu, p, a[,col]) [@LWB out[,col]]
OD;
Line 172:
print polynomial(d)
END # fitting #</lang>
Output:
<pre>
3x**2+2x+1
1.0848x**2+10.3552x-0.6164
</pre>
=={{header|C}}==
{{libheader|libgsl}}
| |||