Como calcular o rotacional (curl) - tiagosutter/symdel GitHub Wiki

Lembrando que antes de utilizar as funções no Octave é necessário carregar o pacote symdel, utilizando pkg load symdel

Os símbolos podem ter qualquer nome, eles devem ser passados para a função sym_curl no segundo argumento. Então para o caso onde os símbolos usados para o sistema esférico são, por exemplo, r, t (ao invés de theta) e fi (ao invés de phi), a função deve ser chamada como sym_curl(V, [r t fi], 'spherical')

Embora a escolha de nomes das variáveis seja algo pessoal é recomendável a utilização de nomes que façam sentido.

Em coordenadas cartesianas

% Declaração de símbolos
syms x y z

% Campo vetorial
V = [3*x^2; 2*y*z; z*x];

% Calcula o rotacional, armazenado-o na variável symCurl 
symCurl = sym_curl(V, [x y z], 'cartesian')

Em coordenadas cilíndricas

% Declaração de símbolos
syms rho phi z

% Campo vetorial
V = [rho^2; cos(phi); 0];

% Calcula o rotacional, armazenado-o na variável symCurl 
symCurl = sym_curl(V, [rho phi z ], 'cylindrical')

Em coordenadas esféricas

% Declaração de símbolos
syms r theta phi

% Campo vetorial
V = [r^(1/2), r*sin(theta), 0];

% Calcula o rotacional, armazenado-o na variável symCurl 
symCurl = sym_curl(V, [r theta phi], 'spherical')