function [X,Y]= lint(a,b,r,rn);
X=linspace(0,r,100);
Y=a+b*X+(rand(1,100)-0.5)*rn;
end;
