Calculo del rotacional de un campo vectorial en coordenadas esfericas con SAGE

Solapas principales


style="display:inline-block;width:336px;height:280px"
data-ad-client="ca-pub-3713658242215030"
data-ad-slot="9156825571">

Calculo del rotacional del campo vectorial  en coordenadas esfericas V=(k/r, 0, 0) con SAGE:
 
var('r th ph k')
def rotacionalesfericas(F):
    assert(len(F) == 3)
    return vector([diff(F[2]*sin(th),th)/(r*sin(th))-diff(F[1],ph)/(r*sin(th)), diff(F[0],ph)/(r*sin(th))-diff(r*F[2],r)/r, diff(r*F[1],r)/r-diff(F[0],ph)/r]) 
 
F1=vector([k/r,0,0])
rotacionalesfericas(F1)
 
 

Resultado:

 
\nabla \times \bar{V} = \left(0,\,0,\,0\right)   El campo vectorial es irrotacional
 
Español

Añadir nuevo comentario

Plain text

  • No se permiten etiquetas HTML.
  • Las direcciones de las páginas web y las de correo se convierten en enlaces automáticamente.
  • Saltos automáticos de líneas y de párrafos.
Pin It