całki różnymi metodami

[aisak7]

Witam,piszę program w liczący całkę korzystając z kwadratur,i mam pewien problem.

Kod: Zaznacz cały

%program liczący całkę kwadraturą
function f=pomoc(a,b,n);
suma=0;
h=(b-a)*n; 
x=a:h:b;
y=exp(x);
for i=1:(n+1)
pole=h/6*(y(i)+4*y(i+1)+y(i+2));
suma=suma+pole
end 
endfunction

I właśnie tu chciałam zapytać co w nim jest nie tak,
co mogłabym w nim poprawić?

[octahedron]

Kod: Zaznacz cały

h=(b-a)*n; %dlugosc podstawy przedziału

a powinno chyba być / zamiast *

Kod: Zaznacz cały

h=(b-a)/n; %dlugosc podstawy przedziału

poza tym

Kod: Zaznacz cały

x=a:(b-a)/(2*n):b;
y=exp(x);
for i=1:2:(2*n-1)