clear;
clc;
f=@(x,a,b)(sin(x)-a-b.*x).^2;
g=@(p)quadl(@(x)f(x,p(1),p(2)),0,pi);
fminunc(g,[1 1])