3826.sagews

def vars(m):  #CONVERTS STRINGS TO VARIABLES#
    m=m.split(',')
    mvar=[]
    for n in m:
        mvar.append(n)
    for i in range(len(m)):
        m[i]='m[i]'+'='+'var("'+m[i]+'")'
        exec(m[i])
    return(m)
def ineq(m,funct):    #Now count sum of partial derivatives(=global derivative)
    for i in range(len(m)):
        m[i]=diff(funct,m[i])==0
    return m
print solve(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] ]

print sin(3*pi*1/9)

1/2*sqrt(3)