program PotencialAnillo
  implicit none
!Física Computacional (Horario 2:00-3:00 pm)
!Tarea 10
!Tema Potencial Eléctrico
!Nalleli Iridian Avila Garcia

  real(8), parameter :: pi = 3.14159265358979323846
  real(8), parameter :: k = 8.9875517873681764e9 ! Constante de Coulomb
  real(8) :: R, Q, r, theta, phi
  real(8) :: V
  integer :: n
  real(8) :: Ksq, K, dphi, integral

  ! Parámetros del anillo
  R = 5.0 ! Radio del anillo en metros
  Q = 10.0 ! Carga total en Coulombs
  r = 1.0 ! Valor de r en metros
  theta = pi / 2.0 ! Ángulo theta en radianes

  ! Inicialización del potencial
  V = 0.0

  ! Cálculo de K(K)
  Ksq = (4.0 * R * r * sin(theta)) / (R**2 + r**2 + 2.0 * R * r * sin(theta))
  K = sqrt(Ksq)

  ! Aproximación numérica de la integral
  n = 10000 ! Número de puntos para la aproximación
  dphi = pi / (2.0 * real(n))

  integral = 0.0
  do n = 0, n - 1
    phi = dphi * real(n)
    integral = integral + dphi / sqrt(1.0 - Ksq * sin(phi)**2)
  end do

  ! Cálculo del potencial
  V = (2.0 * k * Q / sqrt(R * r * sin(theta))) * integral

  ! Imprimir el resultado
  print *, "Potencial en el punto (r=", r, ", theta=", theta, ", phi=", phi, "): V =", V

end program PotencialAnillo