Metodo del trapecio

disp('Metodo del Trapecio');
a=input('valor de a:');
b=input('valor de b:');
n=input('valor de n:');
Ar=0;
h=(b-a)/n;
     for i=0:n-1
         xi=a+i*h;
         xii=a+(1+i)*h;
         S=1/(1-xi)+1/(1-xii);
         if (S<0)
             S=(-1)*S;
         else
             S=S*1;
         end
            
         Ar=Ar+ h/2*(S);
        
     end
fprintf('el valor del area es %.8f',Ar);