Untitled

function ft = RegresjaKwadratowaWart2(x,y,t) % ok
    n=size(x,2); % ilosc punktow danych
    A=zeros(3);
    B=zeros(3,1);
    for i=1:n %obliczamy kolejne wyrazy odpowiednich sum i dodajemy je
        A(1,1)=A(1,1)+x(i)^4;
        A(1,2)=A(1,2)+x(i)^3;
        A(1,3)=A(1,3)+x(i)^2;
        A(2,1)=A(2,1)+x(i)^3;
        A(2,2)=A(2,2)+x(i)^2;
        A(2,3)=A(2,3)+x(i);
        A(3,1)=A(3,1)+x(i)^2;
        A(3,2)=A(3,2)+x(i);
        A(3,3)=A(3,3)+1;
        B(1,1)=B(1,1)+x(i)^2*y(i);
        B(2,1)=B(2,1)+x(i)*y(i);
        B(3,1)=B(3,1)+y(i);
    end
    rozw=eliminacja_Gaussa_cw(A,B); % rozwiazaniem ukladu rownan sa wspolczynniki krzywej
    ft=rozw(1,1).*t.^2+rozw(2,1).*t+rozw(3,1); % wartosc w punktach