clc;clear;

load ener1.data

ener=ener1;

len=length(ener)^(1/3);

len1=51;

max(ener(:,5))

min(ener(:,5))

cc=reshape(ener(:,5),len1,len1,len1);

% image(cc);figure(gcf);

ccccc=0:(len1-1);

[X,Y,Z] = meshgrid(ccccc,ccccc,ccccc);

% slice(X,Y,Z,cc,[25],25,[25])

%%

count=0;

for i=-62:5:175

   cla;

p = patch(isosurface(X,Y,Z,cc,i));

isonormals(X,Y,Z,cc,p)

daspect([1 1 1])

isocolors(X,Y,Z,flipdim(cc,2),p)

shading interp

grid off

axis(volumebounds(X,Y,Z,cc))

axis tight

view(3)

camlight

lighting phong

count=count+1;

name=['U=',num2str(i)];

legend(name,'North');

name=[num2str(count),'.jpg'];

saveas(gcf,name)

end