function []=make_iso(X,Y,Z,C,val)

fv=isosurface(X,Y,Z,C,val);

h=figure;
p=patch(fv);
set(p,'FaceColor','red','Edgecolor','none')
axis equal
view(30,30)
axis tight
camlight('headlight')
set(h,'Renderer','zbuffer');
lighting phong
grid on
set(gca,'fontsize',18,'fontname','times new roman')
xlabel('x')
ylabel('y')
zlabel('z')
drawnow

figfilename=strcat('isosurface_val','.fig');
saveas(h,figfilename)

end