function [p, fval] = solve_cvx_style(x, y)
cvx_begin
variable p(2);
minimize( sum(square(x*p(1) + p(2) - y)) );
cvx_end
fval = cvx_optval;
display(p);
display(fval);
end
