Example 1.  Find the standard "least squares parabola"  [Graphics:Images/LeastSqPolyMod_gr_9.gif]  for the data points [Graphics:Images/LeastSqPolyMod_gr_10.gif].  
Use the subroutine LSParabola to find the line.  Compare with the line obtained with Mathematica's Fit procedure.

Solution 1.

[Graphics:../Images/LeastSqPolyMod_gr_11.gif]



[Graphics:../Images/LeastSqPolyMod_gr_12.gif]

[Graphics:../Images/LeastSqPolyMod_gr_13.gif]
[Graphics:../Images/LeastSqPolyMod_gr_14.gif]
[Graphics:../Images/LeastSqPolyMod_gr_15.gif]
[Graphics:../Images/LeastSqPolyMod_gr_16.gif]
[Graphics:../Images/LeastSqPolyMod_gr_17.gif]

Let's peek at the linear system that was solved.

[Graphics:../Images/LeastSqPolyMod_gr_18.gif]



[Graphics:../Images/LeastSqPolyMod_gr_19.gif]
[Graphics:../Images/LeastSqPolyMod_gr_20.gif]

[Graphics:../Images/LeastSqPolyMod_gr_21.gif]

[Graphics:../Images/LeastSqPolyMod_gr_22.gif]
[Graphics:../Images/LeastSqPolyMod_gr_23.gif]

[Graphics:../Images/LeastSqPolyMod_gr_24.gif]
[Graphics:../Images/LeastSqPolyMod_gr_25.gif]
[Graphics:../Images/LeastSqPolyMod_gr_26.gif]
[Graphics:../Images/LeastSqPolyMod_gr_27.gif]
[Graphics:../Images/LeastSqPolyMod_gr_28.gif]
[Graphics:../Images/LeastSqPolyMod_gr_29.gif]


Of course we want a graph.

[Graphics:../Images/LeastSqPolyMod_gr_30.gif]


[Graphics:../Images/LeastSqPolyMod_gr_31.gif]

[Graphics:../Images/LeastSqPolyMod_gr_32.gif]

[Graphics:../Images/LeastSqPolyMod_gr_33.gif]

[Graphics:../Images/LeastSqPolyMod_gr_34.gif]

The RMS error  [Graphics:../Images/LeastSqPolyMod_gr_35.gif] for this example is:

[Graphics:../Images/LeastSqPolyMod_gr_36.gif]

[Graphics:../Images/LeastSqPolyMod_gr_37.gif]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004