Anonymous user
Polynomial regression: Difference between revisions
m
→{{header|PARI/GP}}
m (J: minor simplification) |
|||
Line 863:
<lang parigp>V=[1,6,17,34,57,86,121,162,209,262,321]~;
M=matrix(#V,3,i,j,(i-1)^(j-1));Polrev(matsolve(M~*M,M~*V))</lang>
<small>Code thanks to [http://pari.math.u-bordeaux.fr/archives/pari-users-1105/msg00006.html Bill Allombert]</small>▼
{{out}}
<pre>3*x^2 + 2*x + 1</pre>
▲<small>Code thanks to [http://pari.math.u-bordeaux.fr/archives/pari-users-1105/msg00006.html Bill Allombert]</small>
Least-squares polynomial fit in its own function:
<lang parigp>lsf(X,Y,n)=my(M=matrix(#X,n,i,j,X[i]^(j-1))); Polrev(matsolve(M~*M,M~*Y~))
lsf([0..10], [1,6,17,34,57,86,121,162,209,262,321], 3)</lang>
| |||