Trzeba zrobić funkcję (coś o aproksymacji) :P . Mam w zeszycie napisane.


x=0:0.05:4;
y1=-cos(x)+0.2*sin(50*x)+0.3*cos(40*x-2);
y2=-cos(x)-0.2*sin(50*x)-0.3*cos(40*x+2);
X=[x x];
Y=[y1 y2];
plot(X,Y,'.b');
p=polyfit(X,Y,1)
YA=polyval(p,X)
plot(X,Y,'.b');
hold on;
plot(X,YA,'.r');



if nargin==0
N=1:10;
end

for n=N