Program restored to include both isotropic and proplyd cases

the program makes theoretical curves for 50 beta values from
beta=0.001 to beta = 0.08