Matlab - Ejercicio con integrales

n=4;

fx=@(x) sqrt(x);

li=1;

ls=2;

dx=(ls -li)/n;

x=linspace(li,ls,n+1);

L4=sum(fx(x(1:n)))*dx;

R4=sum(fx(x(2:n+1)))*dx;

exacto=integral(fx,li,ls);

n=8;

fx=@(x)  x.^2 - 2.*x + 3;

li=-2;

ls=3;

dx=(ls -li)/n;

x=linspace(li,ls,n+1);

L4=sum(fx(x(1:n)))*dx;

R4=sum(fx(x(2:n+1)))*dx;

exacto=integral(fx,li,ls);

n=[8 16 32 48];

Ln=[];

Rn=[];

fx=@(x)  x.^2 - 2.*x + 3;

li=-2;

ls=3;

exacto=integral(fx,li,ls)

for k=1:length(n)

dx=(ls -li)/n(k);

x=linspace(li,ls,n(k)+1);

Ln(k)=sum(fx(x(1:n(k))))*dx;

Rn(k)=sum(fx(x(2:n(k)+1)))*dx;

end