x=[0,1,2,3,4,5,6,7,8,9,10];
 Fy=[15,10,20,15,25,20,30,20,35,30,40];
 if length(x)==length(Fy)
 n=length(x);
 stopien=4;
 if stopien<n-1
 A=polyfit(x,Fy,stopien);
 Wy=zeros(1,n);
 for i=1:n
 for j=1:stopien+1
 Wy(1,i)=Wy(1,i) + A(1,j)*x(1,i)^(stopien+1-j);
 end
 end
 plot(x,Fy,'rx');
 grid on;
 hold on;
 plot(x, Wy);
 disp("Sprawdzenie");
 disp("Blad wzgledny");
 bw=sum(abs(Fy-Wy))
 disp("Blad kwadratowy");
 bk=sum((Fy-Wy).^2)
 disp("Sredni blad procentowy");
 sbp=1/n*sum(abs((Fy-Wy)./Fy))*100
 disp("Srednie odchylenie od krzywej");
 sook=sqrt(sum((Fy-Wy).^2)/(n-stopien-1))
 opis=sprintf("B.w.=%f, B.k.=%f, S.b.p.=%f, sr.od.od.k.=%f",bw, bk, sbp, sook);
 title(opis);
 disp("----------------------------------------------------------------------");
 else
 disp("Za duzy stopien wielomianu aproksymujacego");
 end
 end
 clear;

 x=[0,1,2,3,4,5,6,7,8,9,10];
 Fy=[15,10,20,15,25,20,30,20,35,30,40];
 if length(x)==length(Fy)
 n=length(x);
 stopien=8;
 if stopien<n-1
 A=polyfit(x,Fy,stopien);
 Wy=zeros(1,n);
 for i=1:n
 for j=1:stopien+1
 Wy(1,i)=Wy(1,i) + A(1,j)*x(1,i)^(stopien+1-j);
 end
 end
 plot(x,Fy,'rx');
 grid on;
 hold on;
 plot(x, Wy);

 disp("Sprawdzenie");
 disp("Blad wzgledny");
 bw=sum(abs(Fy-Wy))
 disp("Blad kwadratowy");
 bk=sum((Fy-Wy).^2)
 disp("Sredni blad procentowy");
 sbp=1/n*sum(abs((Fy-Wy)./Fy))*100
 disp("Srednie odchylenie od krzywej");
 sook=sqrt(sum((Fy-Wy).^2)/(n-stopien-1))
 opis=sprintf("B.w.=%f, B.k.=%f, S.b.p.=%f, sr.od.od.k.=%f",bw, bk, sbp, sook);
 title(opis);
 disp("----------------------------------------------------------------------");
 else
 disp("Za duzy stopien wielomianu aproksymujacego");
 end
 end