testmed.m

% script testmed
delete(gcf)
[x,y] = gendata;
[sol, chi2]=mlsq(x,y)
plot(x,y,'r+');
hold on
a = tan(sol(1));
b = -sol(2)/cos(sol(1));
plot(x, a*x+b)