---

ok, espero te sirva

% Integrales definidas
% Realizado por Ing E. Porto
clear;clc
%Ingreso de datos
disp(' Regla de trapecio')
disp('Integrales definidas')
fx=input('digite la funcion f(x) =','s');
a=input('limite inferior = ');
b=input('limite superior = ');
tol=input('tolerancia = ');
%Condiciones iniciales
err(1)=100;ns=1;exito=0;
%Calculo de la integral
while exito==0
h=(b-a)/ns;
x=a:h:b;
y=eval(fx);
Iaprox(ns)=h/2*(2*sum(y)-y(1)-y(ns+1));
if ns>1
err(ns)=abs((Iaprox(ns)-Iaprox(ns-1))/Iaprox(ns))*100; %calculo del error
if err(ns)<tol
exito=1;
break;
end
end
ns=ns+1;
end
%Presentacion de resultados
n=1:ns;
fprintf('\n');
disp([' segm' ' integ' ' error'])
disp([n' Iaprox' err'])
fprintf('\n se alcanzo la convergencia en %g segmentos ',ns);
fprintf('\n la integral aproximada es %g \n',Iaprox(ns));

---

quítale el punto y coma a los renglones

a(j,j)=a(j,j)-m(j,i)*a(i,j);
b(j)=b(j)-b(i)*m(j,i);

y te mostrará las operaciones paso a paso. Saludos. Un poco tarde, pero más vale tarde que nunca.

---

Pronto lo estaré programando; es muy sencillo.