close all
clear all
fid_polinomio=fopen("polinomio.txt","r");
fid_integrale=fopen("integrale.txt","w");

while not(feof(fid_polinomio))
    linea=fgets(fid_polinomio)
    A=str2num(linea)
    k=find(A | not(A));  % vettore degli indici (tutti)
    B=zeros(1,length(k)+1);
    B(2:end)=A./k %Coefficienti di P'(x)
    fprintf(fid_integrale,'%d  ', B);
    fprintf(fid_integrale,'\n');
    fprintf('%d  ', B);
    fprintf('\n');
end

% A=[1 4 6 3.3 9 -6 0 6] %coefficienti di P(x)


fclose(fid_polinomio);