clc; clear all; h=1/2;b2=1/2;l=1;
c=@(the,phi1) abs((cos(pi*h/l*cos(the).*sin(phi1))-cos(pi*h/l))./ \ (sqrt(1-cos(the).^2.*sin(phi1).^2))) \ *2.*abs(sin(pi*2*b2/l*cos(phi1)));
ezsurf(@(the,phi1) c(the,phi1).*cos(the).*sin(phi1),\ @(the,phi1) c(the,phi1).*sin(the).*sin(phi1),\ @(the,phi1) c(the,phi1).*cos(phi1), [-pi, pi, -pi/2, pi/2] ,30);
axis equal; #print("horDipolOPG.png", "-dpng");
:)