x=[0,1,2,3,4,5,6,7,8,9,10];
Fy=[10,5,15,10,20,15,25,20,30,25,35];

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("Bw=%f, Bk=%f, Sbp=%f, Sook=%f",bw, bk, sbp, sook);
 title(opis);

 disp("********************");

 else
 disp("Za duzy stopien wielomianu aproksymujacego");
 end
end