Program that calculates the solutions to a 2nd degree equation of complex numbers

noelia.f90

PROGRAM Q_1 

    implicit none

   

    complex :: a, b, c !! coefficients

    complex :: x1, x2 !! solutions

   

    write (*,*) "Enter a, b, c complex values"

    read (*,*) a, b, c

    write (*,*) "Using values a=", a, "b=", b, "c=", c

    if (a==0) then

        if (b==0) then

            write (*,*) "There is no solution"

        else

            write (*,*) "There is only one solution x1=", -c/b

        end if

    else

        x1 = (-b + sqrt (b**2 - 4*a*c) ) / (2*a)

        x2 = (-b - sqrt (b**2 - 4*a*c) ) / (2*a)
        
        write (*,*) "There are two solutions x1=", x1, "and x2=", x2

    end if
    
    read(*,*)
    write (*,*) "Program terminated"

end program Q_1