changed r/rs to r**2/rs**2 in filament model