[MatLab] somma cicli for

[alvi90]

Ciao!

ho un problema da qui non riesco a venire fuori.
Ho due cicli for, questi cicli mi danno tutti e 9 risultati. Ora ho bisogno di sottrare dal ciclo for 1, il ciclo for 2, in modo da ottenere l'errore.
Il problema è che sottrando mi da 1 solo risultato, e non 9 come vorrei.
qualcuno sa perche??

     
for k=1:1:9;
    n=2.^k;
    h = (b-a)/n;
    sum_even = 0;
    for i = 1:n/2-1
        x(i) = a + 2*i*h;
        sum_even = sum_even + f(x(i));
    end
    sum_odd = 0;
    for i = 1:n/2
        x(i) = a + (2*i-1)*h;
        sum_odd = sum_odd + f(x(i));
    end
    integral = h*(f(a)+ 2*sum_even + 4*sum_odd +f(b))/3;
end


for k=1:1:9;
    n=2*((2.^k)-1)+1;
    h = (b-a)/n;
    sum_even = 0;
    for i = 1:n/2-1
        x(i) = a + 2*i*h;
        sum_even = sum_even + f(x(i));
    end
    sum_odd = 0;
    for i = 1:n/2
        x(i) = a + (2*i-1)*h;
        sum_odd = sum_odd + f(x(i));
    end
    integral2 = h*(f(a)+ 2*sum_even + 4*sum_odd +f(b))/3;
end
P=integral2-integral

a me interesserebe avere 9 valori di P!
Grazie mille a tutti!

[apatriarca]

P è uno scalare e non un vettore.. All'interno dei cicli dovresti scrivere qualcosa come integral(k) = ....