View difference between Paste ID: UQfWzGN3 and ycSUHbDL

-program riemann
+module functions
-implicit none
+contains
   function f(x)
-integer :: i,N
+    implicit none
-real :: integ,a,b,dx,x
+    real :: f
     integer,intent(in) :: x
-! zerando as variaveis
+
-integ = 0.
+    f = EXP(real(x))
-x=0.
+
-dx=0.
+  end function
-
+end module
-print *,"Digite os limites a,b e o numero N de intervalos:"
+
-print *,"Use N grande, ex: 1000"
+program trapezio
-read *,a,b,N
+  use functions
   implicit none
-dx = (b-a)/real(N)
+
   integer :: i, n, b, a
-!integ = 0.5*(f(a)+f(b))
+  real :: sumation, mean, deltax, integral
 
-do i = 1,N
+
-    x = a + i*dx
+  ! The value of the integral using the trapezium rule can be found using
-    integ = integ + f(x)*dx
+  ! integral = (b - a)*((f(a) +f(b))/2 + sumation_1_n-1 )/n
-enddo
+
   write(*,*) "type the limits b, a and the number of intervals"
-print *,'Integral = ', integ
+  read *, b, a, n
-
+
-stop '>>> Fim do Programa'
+  deltax = real(b - a)/real(n)
   mean = (f(a) + f(b))/2
- contains
+  sumation = 0
-real function f(x)
+
-real :: x
+  do i = 1, n-1
-! f(x) = x^2
+    sumation = sumation + f(i)
- f = x**2.
+  end do
-return
+
-end function f
+
   integral = deltax*(mean + sumation)
-end program riemann
+  write (*,*) "the value of the integral using the trapezoidal method is", integral
 
 end program trapezio