How do you write a lowest degree polynomial curve fitting program in matlab.?

[Terra]

I know that matlab provides this in the toolbox but I'd like to see how the program actually works.