Solución: Utilice Matlab para dibujar dos cargas puntuales con el mismo signo, ambas con Q (Q > 0) y separadas por 2a. Dibuje las líneas de campo eléctrico y las líneas equipotenciales.

clc

alpha=5.2182;

n=1.55;

e=0;

f= 0;

para l=1:50

syms x f1 f2

si l==1

f1=sin(x )./x-cos(x);

f1dot=diff(f1,x,l);

f2=(sin(x)+i*cos(x)) ./x-cos(x)+i*sin(x);

f2dot=diff(f2,x,l);

else

f1 =(pi*x/2).*besselj(l+0.5,x);

f1dot=(pi*x/2).*besselj(l-0.5,x)-l*(pi *x/2).*besselj(l+0.5,x)./x;

f2=(pi*x/2).*(besselj(l+0.5,x)-i*besselj (l+0.5,x));

f2dot=(pi*x/2).*(besselj(l-0.5,x)-i*bessely(l-0.5,x))-l *(pi*x/2).*(besselj(l+0.5,x)-i*bessely(l+0.5,x))./x;

end

g1=inline(vectorizar(f1), x);

g1dot=inline(vectorizar(f1dot), x);

g2=inline(vectorizar(f2), x);

g2dot=inline(vectorize(f2dot), x);

a1=(g1dot(n*alpha).*g1(alpha)-n*g1(n*alpha) .*g1dot(alpha))./(g1dot(n*alpha).*g2(alpha)-n*g1(n*alpha).*g2dot(alpha));

b1=(n *g1punto(n*alfa).*g1(alfa)-g1(n*alfa).*g1punto(alfa))./(n*g1punto(n*alfa).*g2(alfa)-g1(n*alfa ).*g2punto(alfa));

e=e+(2*l+1)*(abs(a1).^2+abs(b1).^2);

f=f+(2*l+1)*real(a1+b1);

end

kext=2/alpha.^2.*f

ksca=2/alpha.^2.*e