PROGRAM radici
!Scopo: Questo programma calcola le radici di un’equazione
! a*x**2+ b*x + c = 0
! *** SEZIONE DICHIARATIVA *** ! 
IMPLICIT NONE
REAL :: a, b, c
REAL :: discriminante, parte_immag, parte_reale
REAL :: x1, x2
! Soluzioni dell’equazione
! *** SEZIONE ESECUTIVA *** !
WRITE (*,*) 'Inserire i tre coeefficienti A, B e C: '
READ (*,*) a, b, c
!Calcola il discriminante
discriminante = b**2 - 4 * a * c
IF (discriminante >0) THEN
   x1 = (-b + sqrt(discriminante)) / (2.*a)
   x2 = (-b - sqrt(discriminante)) / (2.*a)
   WRITE (*,*) 'Due radici reali: ', x1, x2
ELSE IF (discriminante == 0.) THEN
   x1 = ( -b ) / (2. * a)
   WRITE (*,*) 'due radici reali e coincidenti: ', x1
ELSE    ! Radici complesse
   parte_reale = ( -b ) / (2. * a)
   parte_immag = sqrt(abs(discriminante)) / (2.*a)
   WRITE (*,*) 'due radici immaginarie:', parte_reale, '+/- i', parte_immag
END IF
STOP 
END PROGRAM radici