clear all;
clc;

syms x y z a b c N
A = [a;b;c];
R = [x;y;z];
AR = a*x+b*y+c*z;
AxR = cross(A,R);
r = sqrt(x^2+y^2+z^2);
A2 = (a^2+b^2+c^2);
AE = AR/r;

rot = curl(AR*AxR/(r^N));
rot_szamolt = ((3-N)*AR*A+(N*AE^2-A2)*R)/(r^N);

simplify(rot - rot_szamolt) % juhé, a különbség nullvektor