Algoritmo da Regra dos Trapézios para resolução de Integrais

Início
entre com os limites inferior e superior da integral [x0,xn]
entre com o numero de subintervalos
Se n = 0 Então
      Escreva(‘Divisão por zero’)
    Senão
         Se n < 0 Então
        Escreva(‘Intervalo Inválido’)
           Senão
   Início
        h = (xn – x0) / n
        x = x0 + h
        soma = 0
Para i = 1 Até n – 1 Faça
    Início
     soma=soma + f(x)
           x =x + h
Fim-Para
   R= h * (( f(x0) + f(xn) ) / 2 + soma)
      Escreva(‘O resultado da integral da função f é ’, R)
Fim-Senão
Fim.
------------------------------------
Variáveis utilizadas no algoritmo:
· Reais: h, x, soma, R;
· Inteiras: i, n.

!!!!!!!!!!!!!!!!!!!!!! programa !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!11
!Universidade Estadual do Ceará
!Curso de Física
!Computação Aplicada a Física I
!Prof. Dr. Carlos Santana
!Programa baseado no algoritmo do livro Análise Numerica - Richad l. Burden e J. Douglas Faires integral de trapezio
 program trapezio_composta
 real h,x0,xn,soma,r
 integer n
   print*,'entre com o valor das extremidades'
 read*,x0,xn
5 print*,'entre com o numero de subintervalos'
 read*,n
 if(n==0)then
    print*,'erro, divisão por zero'
    print*,' tente novamente '
    print*
        goto 5
 endif
 if(n<0 p="" then="">    print*,'erro, intervalo invalido'
    print*,' tente novamente '
    print*
        goto 5
 endif
  h=(xn-x0)/float(n)
 x=x0+h
 soma=0
 do i=1,n-1
     soma=soma+f(x)
     x=x+h
 enddo
 r=h*((f(x0)+f(xn))/2.+soma)
 write(*,*)"integral =",r
 stop
 end
 real function f(x)
 !exemplo
 !x0=1  xn=4
  f=1/x*exp(x/2)
  return
 end