function [err, grad] = lserr(x, y, param)
    err = mean( (param(1)*x + param(2) - y).^2 );
    
    if nargout == 2
        grad = lserr_gradient(param, x,y);
    end