clear
F=load ('DANE7.PAS');
[n,p]=size(F);
    W=@(x) x^3-2*x^2+5*x-6; m=3;
    S=zeros(m+1,m+1);
    T=zeros(m+1,1);
    ms=0;
    potega=0;
    if(m<n)
        for i=1:m+1
            pom=ms;
            for j=1:m+1
                sum_s=0
                 for k=1:n
                   sum_s=sum_s+(F(k,1)^ms);
               end
                S(i,j)=sum_s;
                ms=ms+1;
            end
            ms=pom;
            ms=ms+1;
            sum_t=0;
            for k=1:n
                sum_t=sum_t+((F(k,1)^(i-1)*F(k,2)));
                end
          T(i,1)=sum_t;
        end
      
    end
    
    S
    T