clear
res=load('ruther.dat');

[x,y,z]=sphere;
hold on
plot(res(:,2),res(:,3))
