program electricForce integer,parameter :: nCharges = 3 ! numero de cargas real,dimension(2,nCharges) :: chargesPositions ! las columnas de la mat, referncia las cargas ! la fila 1 la pos X, la fila 2 la pos Y real,dimension(nCharges) :: charges ! cargas electricas de las cargas real,dimension(nCharges) :: electricForces ! Fuerza electrica sibre las cargas integer :: q1,q2 real :: num,dem ! ---instancio las variables iniciales --- charges(1:nCharges) = (/1.0,1.5,1.3/) chargesPositions(1,1:nCharges) = (/0.5,-0.5,0.0/) chargesPositions(2,1:nCharges) = (/1.0,1.0,-1.0/) do q1 = 1,nCharges electricForces(q1) = 0 do q2 = 1,nCharges if ( q1 /= q2) then num = charges(q1)*charges(q2) den = (chargesPositions(1,q1)-chargesPositions(1,q2))**2 + (chargesPositions(2,q1)-chargesPositions(2,q2))**2 electricForces(q1) = electricForces(q1) + num/den end if end do end do print *, "La suma de los módulos de la fuerza electrica para c/ carga es: ", electricForces end program electricForce