1 #include "leastsquares.hh"
5 Least_squares::OK() const
7 assert (input.size() > 1);
9 for (int i=1; i < input.size(); i++)
10 dx += abs (input[i-1].x() - input[i].x ());
16 Least_squares::minimise (Real &coef, Real &offset)
23 for (int i=0; i < input.size();i++)
26 Real y = input[i].y();
37 Real den = (N*sqx - sqr (sx));
39 programming_error ("Least_squares::minimise(): Nothing to minimise");
40 coef = (N * sxy - sx*sy)/den;
41 offset = (sy - coef * sx)/N;