> points := [];
> points := [op(points), [499, 5146.0]];
> points := [op(points), [3000, 41612.0]];
> points := [op(points), [5000, 75525.0]];
> points := [op(points), [9500, 158435.0]];
> points := [op(points), [20000,367527.0]];
> points := [op(points), [20000,350438.0]];
> points := [op(points), [20000,343021.0]];
> points := [op(points), [20000,333726.0]];
> points := [op(points), [20000,366864.0]];
> 
> 
> for i from 1 to nops(points) do
>    x[i] := points[i,1]; y[i] := points[i,2];
> od:
> 
> x_coords := []; y_coords:=[];
> for i from 1 to nops(points) do
>     x_coords := [op(x_coords), x[i]];
>     y_coords := [op(y_coords), y[i]];
> od:

                             points := []


                      points := [[499, 5146.0]]


              points := [[499, 5146.0], [3000, 41612.0]]


     points := [[499, 5146.0], [3000, 41612.0], [5000, 75525.0]]


  points := [[499, 5146.0], [3000, 41612.0], [5000, 75525.0],

        [9500, 158435.0]]


  points := [[499, 5146.0], [3000, 41612.0], [5000, 75525.0],

        [9500, 158435.0], [20000, 367527.0]]


  points := [[499, 5146.0], [3000, 41612.0], [5000, 75525.0],

        [9500, 158435.0], [20000, 367527.0], [20000, 350438.0]]


  points := [[499, 5146.0], [3000, 41612.0], [5000, 75525.0],

        [9500, 158435.0], [20000, 367527.0], [20000, 350438.0],

        [20000, 343021.0]]


  points := [[499, 5146.0], [3000, 41612.0], [5000, 75525.0],

        [9500, 158435.0], [20000, 367527.0], [20000, 350438.0],

        [20000, 343021.0], [20000, 333726.0]]


  points := [[499, 5146.0], [3000, 41612.0], [5000, 75525.0],

        [9500, 158435.0], [20000, 367527.0], [20000, 350438.0],

        [20000, 343021.0], [20000, 333726.0], [20000, 366864.0]]


                            x_coords := []


                            y_coords := []

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

    model := y = 1.397843910 x ln(x) + 3.887193643 x - 2200.954930

> subs(x=100000, op(2, model));

                 139784.3910 ln(100000) + 386518.4094

> evalf(%);
> 

                                         7
                           .1995845683 10

>