Funkce potenciální energie

U=(1/2)(Subscript[k, 1])((((Sqrt[(((a+x))^2)+(y^2)])-(Subscript[a, 0])))^2)+(1/2)(Subscript[k, 1])((((Sqrt[(((a-x))^2)+(y^2)])-(Subscript[a, 0])))^2)+(1/2)(Subscript[k, 2])((((Sqrt[(x^2)+(((a-y))^2)])-(Subscript[a, 0])))^2)+(1/2)(Subscript[k, 2])((((Sqrt[(x^2)+(((a+y))^2)])-(Subscript[a, 0])))^2)

Derivace potenciálu -- spočtení druhé derivace v bodech x = 0 a y = 0

(Subscript[U, 11])=D[U,x,x]/.{x->0,y->0}
(Subscript[U, 12])=D[U,x,y]/.{x->0,y->0}
(Subscript[U, 21])=(Subscript[U, 12])
(Subscript[U, 22])=D[U,y,y]/.{x->0,y->0}

Matice potenciální energie

mU=MatrixForm[FullSimplify[{{(Subscript[U, 11]),(Subscript[U, 12])},{(Subscript[U, 21]),(Subscript[U, 22])}},Assumptions->{a>0}]]