% E3DTH - program to plot grdE2 3 D for Thesis for Interdigitated
% elctrodes  D. Bakewell 20:30 hrs Fr/9/2/2001 
clear all; desktop;  d=10; d1=d; d2=2*d; Vo=1; 
Ximax=300; Yimax=10; Xmn=-d/2; Ymn=d/100; Xmx=4*d-d/2; Ymx=d; 
Ex_=[];Ey_=[]; 
for yi=1:1:Yimax+1; y=Ymn+(yi-1)*(Ymx-Ymn)/Yimax;
 for xi=1:1:Ximax+1; x=Xmn+(xi-1)*(Xmx-Xmn)/Ximax;
 xh=pi*x/d2+pi/4; yh=pi*y/d2; % h = 'hat'
 Ex=(2.*Vo/(pi*d1))*(atan(sin(xh)/sinh(yh))-atan(cos(xh)/sinh(yh)));
 Ey=(Vo/(pi*d1))*(log(cosh(yh)+cos(xh))-log(cosh(yh)-cos(xh))+...
    log(cosh(yh)+sin(xh))-log(cosh(yh)-sin(xh)));
 Dpos=cosh(2*yh)+cos(2*xh); Dneg=cosh(2*yh)-cos(2*xh);
 Exx=(2.*Vo*sinh(yh)/(d1^2))*(cos(xh)/Dneg+sin(xh)/Dpos); Eyy=-Exx;
 Exy=(2.*Vo*cosh(yh)/(d1^2))*(cos(xh)/Dpos-sin(xh)/Dneg); Eyx=Exy;
 Ex_(xi,yi)=Ex; Ey_(xi,yi)=Ey; Exx_(xi,yi)=Exx; Exy_(xi,yi)=Exy; 
 dxE2_(xi,yi)=2.*(Ex*Exx+Ey*Eyx); dyE2_(xi,yi)=2.*(Ex*Exy+Ey*Eyy);
 end % end for xi
end % end for yi
choice=2;
switch choice
   case 1
E2=Ex_.^2+Ey_.^2; 
surf(Ex_'); dummy=menu('Plot of Ex     ','Press to continu');close;
surf(Ey_'); dummy=menu('Plot of Ey     ','Press to continu');close;
surf(Exx_'); dummy=menu('Plot of Exx    ','Press to continu');close;
surf(Exy_'); dummy=menu('Plot of Exy    ','Press to continu');close;
surf(dxE2_');dummy=menu('Plot of dxE2   ','Press to continu');close;
surf(dyE2_');dummy=menu('Plot of dyE2   ','Press to continu');close;
   case 2
xi=1:1:Ximax+1; yi=1:1:Yimax+1; x=Xmn+(xi-1)*(Xmx-Xmn)/Ximax;
y=Ymn+(yi-1)*(Ymx-Ymn)/Yimax; xh=pi*x/d2+pi/4; yh=pi*y/d2;
hx=surf(x,y,dxE2_'); light('Position',[-10,-20,5]); light('Position',[-10,-20,0]);
light('Position',[-10,0,0]); light('Position',[0,0,-5]); lighting flat;
set(gca,'DataAspectRatio',[50 50 1]);
set(hx,'FaceColor','interp','EdgeColor','k') % interp
shading flat; % get(hx)
dummy=menu('Plot of dxE2   ','Press to continu');close;
hy=surf(x,y,dyE2_'); light('Position',[-10,-20,5]); light('Position',[-10,-20,0]);
light('Position',[-10,0,0]); light('Position',[80,30,5]); lighting flat;
set(gca,'DataAspectRatio',[50 50 1]); 
set(hy,'FaceColor','interp','EdgeColor','interp') %
shading flat; %get(hy)
dummy=menu('Plot of dyE2   ','Press to continu');close;
dxE2Vec(xi,1)=dxE2_(xi,1);plot(x,dxE2Vec); box off;
dummy=menu('dxE2 @ yi=1  ','Press to continu');close;
dyE2Vec(xi,1)=dyE2_(xi,1);plot(x,dyE2Vec); box off;
dummy=menu('dyE2 @ yi=1  ','Press to continu');close;
	end % end switch choice
cd e:
close;