0

No logro pasar este código de Python a Fortran. Estoy usando el método de Simpson para calcular una doble integral. Mi código en Python funciona bien, pero no logro que quede bien en Fortran. Este es mi código en Python

import math

def f(x, y):
    return math.exp(x + y)

def approximate_double_integral(f, a, b, c, d, n, m):
    hx = (b - a) / (2 * n)
    hy = (d - c) / (2 * m)

    integral_sum = 0

    for i in range(n):
        for j in range(m):
            x0 = a + 2 * i * hx
            x1 = a + (2 * i + 1) * hx
            x2 = a + (2 * i + 2) * hx

            y0 = c + 2 * j * hy
            y1 = c + (2 * j + 1) * hy
            y2 = c + (2 * j + 2) * hy

            integral_sum += f(x0, y0) + 4 * f(x1, y0) + f(x2, y0)
            integral_sum += 4 * (f(x0, y1) + 4 * f(x1, y1) + f(x2, y1))
            integral_sum += f(x0, y2) + 4 * f(x1, y2) + f(x2, y2)

    integral_value = (hx * hy / 9) * integral_sum
    return integral_value

a = 1.0
b = 3.0
c = 0.0
d = 4.0
m = 6
n = 4

result = approximate_double_integral(f, a, b, c, d, n, m)
print("El valor aproximado de la integral doble es:", result)

Añado ahora mi código en Fortran. He intentado cambiar los errores que salen pero me los vuelve a marcar. La mayoría es por algún argumento pero a pesar de agregarlo no corre.

program ApproximateDoubleIntegral
    implicit none
    real :: a, b, c, d, hx, hy
    integer :: n, m
    real :: integral_sum, integral_value

    ! Definimos la función f(x, y)
    interface
        real function f(x, y)
            real, intent(in) :: x, y
        end function f
    end interface

    ! Definimos los parámetros de la aproximación de la integral doble
    a = 1.0
    b = 3.0
    c = 0.0
    d = 4.0
    n = 4
    m = 6

    ! Calculamos los tamaños de los subintervalos
    hx = (b - a) / (2.0 * real(n))
    hy = (d - c) / (2.0 * real(m))

    ! Inicializamos la suma de la integral
    integral_sum = 0.0

    ! Calculamos la suma de la integral aproximada
    do i = 0, n - 1
        do j = 0, m - 1
            x0 = a + 2 * i * hx
            x1 = a + (2 * i + 1) * hx
            x2 = a + (2 * i + 2) * hx

            y0 = c + 2 * j * hy
            y1 = c + (2 * j + 1) * hy
            y2 = c + (2 * j + 2) * hy

            integral_sum = integral_sum + f(x0, y0) + 4.0 * f(x1, y0) + f(x2, y0)
            integral_sum = integral_sum + 4.0 * (f(x0, y1) + 4.0 * f(x1, y1) + f(x2, y1))
            integral_sum = integral_sum + f(x0, y2) + 4.0 * f(x1, y2) + f(x2, y2)
        end do
    end do

    ! Calculamos el valor de la integral aproximada
    integral_value = (hx * hy / 9.0) * integral_sum

    ! Imprimimos el valor de la integral aproximada
    write(*,*) "El valor aproximado de la integral doble es: ", integral_value
    pause

end program ApproximateDoubleIntegral

! Definimos la función f(x, y)
real function f(x, y)
    implicit none
    real, intent(in) :: x, y

    f = exp(x + y)

    return
end function f
3
  • 1
    Hola mario, puedes pasar tu codigo actual en Fortran? Asi lo podemos revisar y ver que error estas teniendo. Saludos y bienvenido!
    – leandroian
    Commented el 6 jul. 2023 a las 4:27
  • Puede ser que tengas un problema de precision? Por que no probas en definir tus reales como real(dp)? y cuando lo inicializas le pones el _dp Algo asi real(dp) :: a, b, c, d, hx, hy integer :: n, m real(dp) :: integral_sum, integral_value ! Definimos la función f(x, y) interface real function f(x, y) real, intent(in) :: x, y end function f end interface ! Definimos los parámetros de la aproximación de la integral doble a = 1.0_dp b = 3.0_dp c = 0.0_dp d = 4.0_dp n = 4 m = 6 Cuales son los errores que te salen?
    – leandroian
    Commented el 6 jul. 2023 a las 5:03
  • El motivo por el que tu programa no compila es que no sabe de qué tipo son las variables x0, x1, x2 y sus correspondientes y. Para que el programa corra basta con añadirlas en la declaración de las variables reales.
    – Jorge Luis
    Commented el 7 jul. 2023 a las 9:52

1 respuesta 1

0

Aqui esta la respuesta, faltaron definir las variables. tambien aumente la precision.PAUSE no era necesario. Se comparo el tiempo de ejecucion entre Python y Fortran.

codigo.f90

program ApproximateDoubleIntegral
    implicit none
    real*8 :: a, b, c, d, hx, hy
    integer :: n, m
    real*8 :: integral_sum, integral_value
!!!! falto defenir
    integer :: i,j
    real*4 :: x0,x1,x2,y0,y1,y2
!!!     
    ! Definimos la función f(x, y)
    interface
        real function f(x, y)
            real, intent(in) :: x, y
        end function f
    end interface
    ! Definimos los parámetros de la aproximación de la integral doble
    a = 1.d0 ; b = 3.d0 ; c = 0.d0 ; d = 4.d0
    n = 4 ; m = 6
    ! Calculamos los tamaños de los subintervalos
    hx = (b - a) / (2.d0 * real(n))
    hy = (d - c) / (2.d0 * real(m))
    ! Inicializamos la suma de la integral
    integral_sum = 0.d0
    ! Calculamos la suma de la integral aproximada
    do i = 0, n - 1
        do j = 0, m - 1
            x0 = a + 2 * i * hx
            x1 = a + (2 * i + 1) * hx
            x2 = a + (2 * i + 2) * hx
            y0 = c + 2 * j * hy
            y1 = c + (2 * j + 1) * hy
            y2 = c + (2 * j + 2) * hy
            integral_sum = integral_sum + f(x0, y0) + 4.d0 * f(x1, y0) + f(x2, y0)
            integral_sum = integral_sum + 4.0 * (f(x0, y1) + 4.d0 * f(x1, y1) + f(x2, y1))
            integral_sum = integral_sum + f(x0, y2) + 4.d0 * f(x1, y2) + f(x2, y2)
        end do
    end do
    ! Calculamos el valor de la integral aproximada
    integral_value = (hx * hy / 9.d0) * integral_sum

    ! Imprimimos el valor de la integral aproximada
    write(*,*) "El valor aproximado de la integral doble es: ", integral_value
end program ApproximateDoubleIntegral
! Definimos la función f(x, y)
real function f(x, y)
    implicit none
    real, intent(in) :: x, y
    f = exp(x + y)
    return
end function f

Comparacion de resultados y tiempo,

time python3 aprox.py

El valor aproximado de la integral doble es: 930.935806795878

real 0m0,045s

user 0m0,033s

sys 0m0,012s

gfortran -o integral codigo.f90

time ./integral

El valor aproximado de la integral doble es: 930.93585548577482

real 0m0,004s

user 0m0,001s

sys 0m0,004s

Tu Respuesta

By clicking “Publica tu respuesta”, you agree to our terms of service and acknowledge you have read our privacy policy.

¿No es la respuesta que buscas? Examina otras preguntas con la etiqueta o formula tu propia pregunta.