function param = estimate_regression(x,y, param0)
    options = optimoptions(@fminunc,'GradObj','on');
    myfunc = @(p)(lserr(p,x,y));
    
    param = fminunc(myfunc,param0, options);