PRO emd2

r0 = 1.05
FOR i = 0, 10 DO	$
   BEGIN
      r = r0 + i * 0.1
      gcm2 = emdsun (r)
      Nelec = gcm2 / 2.0e-24
      print, 'r: ', r, '   Nelec: ', nelec, '   emass: ', gcm2
   END
END