Calculo del rotacional del campo vectorial en coordenadas cilindricas V=(r*cos(φ), 2*r*cos(φ), 0)con SAGE:
var('r ph z')
def rotacionalcilindricas(F):
assert(len(F) == 3)
return vector([diff(F[2],ph)/r-diff(F[1],z), diff(F[0],z)-diff(F[2],r), diff(r*F[1],r)/r-diff(F[0],ph)/r])
F1=vector([r*sin(ph),2*r*cos(ph),0])
rotacionalcilindricas(F1)
Añadir nuevo comentario