%sphere.m   main program
nref=4
nvmax=12*4^nref
ntmax=20*4^nref
x=zeros(nvmax,3);
v=zeros(ntmax,3);
icosahedron
for nr=1:nref
  refine
  triplot
end
