h = figure(1); set(h,'Position',[100 678 600 400]);

L = 25; Lx = L; Ly = L; Lz = L; N = 2^8; epsilon = 0.5;
x       = linspace(-Lx,Lx,N+1); x = x(1:end-1); dx = 2*Lx/N;
y       = linspace(-Ly,Ly,N+1); y = y(1:end-1); dy = 2*Ly/N;
z       = linspace(-Lz,Lz,N+1); z = z(1:end-1); dz = 2*Lz/N;
[x,y,z] = meshgrid(x,y,z);

colum = 3;

load Data/Older/TwoPlatesLoose_NoCharged.mat
phi = A{1,colum};
subplot(1,3,1)
p = patch(isosurface(x,y,z,phi,0.5));
isonormals(x,y,z,phi,p)
set(p,'FaceColor','red','EdgeColor','none');
view([1,.2,.3]); axis equal; axis off; axis tight; camlight;lighting gouraud

load Data/Older/TwoPlatesLoose_Charged_q1p2_q2p2.mat
phi = A{1,colum};
subplot(1,3,2)
p = patch(isosurface(x,y,z,phi,0.5));
isonormals(x,y,z,phi,p)
set(p,'FaceColor','red','EdgeColor','none');
view([1,.2,.3]); axis equal; axis off; axis tight; camlight;lighting gouraud


load Data/Older/TwoPlatesLoose_Charged_q1n2_q2p2.mat
phi = A{1,colum};
subplot(1,3,3)
p = patch(isosurface(x,y,z,phi,0.5));
isonormals(x,y,z,phi,p)
set(p,'FaceColor','red','EdgeColor','none');
view([1,.2,.3]); axis equal; axis off; axis tight; camlight;lighting gouraud

Published with MATLAB® R2012b