https://pyclip.ru/en/tapes/FldHNQAEVcc/13-programming-with-modern-fortran-computing-discriminant