
D'après le cours, on sait que si $a+b=S$ et $ab=P$ alors $a$ et $b$ sont racines de $x^2-Sx+P$, dont le discriminant est $\Delta=S^2-4P$.

On peut alors proposer la fonction suivante:
\begin{python}
def f(S,P):
    delta = S*S - 4*P
    if delta > 0:
        return (S - delta**0.5)/2 , (S + delta**0.5)/2
    elif delta == 0:
        return S/2
    else:
        return 'Pas de solution.'
\end{python}
