d(sum_(i=1)^n ((ax_i+b-y_i)^2))/db