其中
a0 =1/k[ ( b2* b4 - b3^2)* c1 + ( b2* b3 - b1 *b4) *c2 +( b1 *b3 - b2^2) c3 ] ,
a1 =1/k[ ( b2 *b3 - b1 *b4)* c1 + + ( n*b4 - b2^2) *c2 +( b1 *b2 - n*b3) c3) ]
a2 =1/k[ ( b1* b3 - b2^2)* c1 + ( b1 *b2 - nb3) *c2 +( n*b2 - b1^2) c3 ]
其中k=n*b2* b4 +2*b1* b2 *b3 - n*b3^2- b1^2*b4 - b2^3
b1 =Σx i ,b2 =Σxi^2 ,b3 =Σxi^3 ,b4 =Σxi^4
c1 =Σyi ,c2 =Σxi*yi ,c3 =Σxi^2 *yi
i=1 to n
a0 =1/k[ ( b2* b4 - b3^2)* c1 + ( b2* b3 - b1 *b4) *c2 +( b1 *b3 - b2^2) c3 ] ,
a1 =1/k[ ( b2 *b3 - b1 *b4)* c1 + + ( n*b4 - b2^2) *c2 +( b1 *b2 - n*b3) c3) ]
a2 =1/k[ ( b1* b3 - b2^2)* c1 + ( b1 *b2 - nb3) *c2 +( n*b2 - b1^2) c3 ]
其中k=n*b2* b4 +2*b1* b2 *b3 - n*b3^2- b1^2*b4 - b2^3
b1 =Σx i ,b2 =Σxi^2 ,b3 =Σxi^3 ,b4 =Σxi^4
c1 =Σyi ,c2 =Σxi*yi ,c3 =Σxi^2 *yi
i=1 to n
dim x(n+1)
b1=0
b2=0
For i=1 to n
b1=b1+x(i)
b2=b2+x(i)^2
.....
next i
k=n*b2* b4 +2*b1* b2 *b3 - n*b3^2- b1^2*b4 - b2^3
a0 =1/k[ ( b2* b4 - b3^2)* c1 + ( b2* b3 - b1 *b4) *c2 +( b1 *b3 - b2^2) c3 ] ,
a1 =1/k[ ( b2 *b3 - b1 *b4)* c1 + + ( n*b4 - b2^2) *c2 +( b1 *b2 - n*b3) c3) ]
a2 =1/k[ ( b1* b3 - b2^2)* c1 + ( b1 *b2 - nb3) *c2 +( n*b2 - b1^2) c3 ]