function figure_mesh(iso_vf,v,d,FColor,EColor)

% plot isosurface

hold off
p=patch(iso_vf);
set(p,'FaceColor',FColor,'EdgeColor',EColor);
view(140,20)
axis tight
camlight 
lighting gouraud
alpha(0.5)

% plot eigenvectors

c = [0,0,1;0,1,0;1,0,0];
hold on
for i=1:3
    quiver3(0,0,0,v(1,i),v(2,i),v(3,i),d(i,i)/8,'Color',c(i,:));
end
hold off
