> f := proc(x) x^2 - x + 1; end;

                     f := proc(x) x^2 - x + 1 end

> rando := rand(-10^5 .. 10^5);

  rando := proc()
local t;
global _seed;
    _seed := irem(427419669081*_seed, 999999999989);
    t := _seed;
    irem(t, 200001) - 100000
end

> for i from 1 to 10 do
>     x[i] := i;
>     y[i] := f(i) + rando()/(10.0^5)*x[i]/10;
> od;
> 

                              x[1] := 1


                         y[1] := 1.031994000


                              x[2] := 2


                         y[2] := 2.975450000


                              x[3] := 3


                         y[3] := 7.166623000


                              x[4] := 4


                         y[4] := 13.28918000


                              x[5] := 5


                         y[5] := 21.13348500


                              x[6] := 6


                         y[6] := 30.98072800


                              x[7] := 7


                         y[7] := 42.73242500


                              x[8] := 8


                         y[8] := 57.05532800


                              x[9] := 9


                         y[9] := 73.83168100


                             x[10] := 10


                         y[10] := 91.09528000

> points := [];
> for i from 1 to 10 do
>     points := [op(points), [x[i], y[i]]];
> od:
> 
> p1 := plot(points, style=point,color = blue):
> 
> x_coords := []; y_coords:=[];
> for i from 1 to 10 do
>     x_coords := [op(x_coords), x[i]];
>     y_coords := [op(y_coords), y[i]];
> od:
> 

                             points := []


                            x_coords := []


                            y_coords := []

> model := stats[fit,leastsquare[[x,y],y=b*x + c]]([x_coords, y_coords]
> );

               model := y = 10.02535572 x - 21.01023907

> p2 := plot(op(2, model), x=0..10, color = red):
> 
> plots[display]({p1, p2});
> 

> model := stats[fit,leastsquare[[x,y],y=a*x^2 + b*x + c]]([x_coords,
> y_coords]
> );
> 

                                2
      model := y = 1.006190481 x  - 1.042739570 x + 1.125951517

> p2 := plot(op(2, model), x=0..10, color = red):
> 
> plots[display]({p1, p2});
> 

> 
> model := stats[fit,leastsquare[[x,y],y=d*x^3 + a*x^2 + b*x +
> c]]([x_coords, y_coords]
> );
> 
> p2 := plot(op(2, model), x=0..10, color = red):
> 
> plots[display]({p1, p2});
> 

                              3                2
  model := y = .002653205193 x  + .9624125954 x  - .8408306553 x

         + .8983065111


>