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.'