Lokální extrémy funkce

1
y=var('y')
2
def loc_extr(FUN):
3
    
4
    pretty_print("Zadaná funkce:", FUN)
5
    sols=solve([diff(FUN,x),diff(FUN,y)],[x,y], solution_dict=True)
6
    real_sols = [sol for sol in sols if all(v in RR for v in sol.values())]
7
    fxx=diff(FUN,x,x)
8
    H=FUN.hessian().determinant()
9
    for s in real_sols:
10
        
11
        if H.subs(s) < 0:
12
            pretty_print("[",s[x],",",s[y],"] = sedlový bod", )
13
        elif H.subs(s) > 0:
14
            druha_derivace=fxx.subs(s)
15
            if fxx.subs(s) > 0:
16
                pretty_print("[",s[x],",",s[y],"] = lokální minimum", )
17
            else:
18
                pretty_print("[",s[x],",",s[y],"] = lokální maximum", )
19
        else:
20
            pretty_print("[",s[x],",",s[y],"] = nelze rozhodnout", )
21