defvars(m):#CONVERTS STRINGS TO VARIABLES#m=m.split(',')mvar=[]forninm:mvar.append(n)foriinrange(len(m)):m[i]='m[i]'+'='+'var("'+m[i]+'")'exec(m[i])return(m)defineq(m,funct):#Now count sum of partial derivatives(=global derivative)foriinrange(len(m)):m[i]=diff(funct,m[i])==0returnmprintsolve(ineq(vars('x,y,z'),(x^2+y^2+z^2)^2-3*(x^3*y+y^3*z+z^3*x)),x,y,z)
[
[x == r5, y == r5, z == r5],
[x == 0, y == 0, z == 0]
]