program quadratic_equation
implicit none
real :: a, b, c
real :: d, root1, root2
d = b**2 - 4.0*a*c
read(*, *) a, b, c
root1 = (-b + sqrt(d)) / (2.0*a)
root2 = (-b - sqrt(d)) / (2.0*a)
print *, 'x1 = ', root1
print *, 'x2 = ', root2
end program quadratic_equation
cHJvZ3JhbSBxdWFkcmF0aWNfZXF1YXRpb24KaW1wbGljaXQgbm9uZQogICAgcmVhbCA6OiBhLCBiLCBjCiAgICByZWFsIDo6IGQsIHJvb3QxLCByb290MgogICAgZCA9IGIqKjIgLSA0LjAqYSpjCiAgICAKICAgIHJlYWQoKiwgKikgYSwgYiwgYwogICAgIHJvb3QxID0gKC1iICsgc3FydChkKSkgLyAoMi4wKmEpCiAgICAgICAgcm9vdDIgPSAoLWIgLSBzcXJ0KGQpKSAvICgyLjAqYSkKICAgICAgICBwcmludCAqLCAneDEgPSAnLCByb290MQogICAgICAgIHByaW50ICosICd4MiA9ICcsIHJvb3QyCiAgICAgICAgZW5kIHByb2dyYW0gcXVhZHJhdGljX2VxdWF0aW9uCiAgICA=