Ответ:
Python @3

import cmath

def solve_quadratic(a, b, c):

   # Вычисляем дискриминант

   D = (b**2) - (4*a*c)

   # Находим два решения

   sol1 = (-b-cmath.sqrt(D))/(2*a)

   sol2 = (-b+cmath.sqrt(D))/(2*a)

   return sol1, sol2

# Коэффициенты

a = 1

b = 16

c = 15

sol1, sol2 = solve_quadratic(a, b, c)

print("Решения квадратного уравнения {0}x^2 + {1}x + {2} являются {3} и {4}".format(a, b, c, sol1, sol2))

Эта программа сначала вычисляет дискриминант, а затем находит 2 решения квадратного уравнения.