def nu(r):
'''Returns the stellar density function.'''
return 1/(r * (1 + (r/a))**3)
mass_int = lambda r: 4 * r**2 * nu(r)
print(mass_int(0))
이것은 1/r 용어가 분리되어 평가되기 때문에 아마도 제로로 나누기 오류를줍니다. 이 문제를 해결하는 유일한 방법은 정확한 대수 표현식을 형성하기 위해 sympy를 사용하고 있습니까? 터무니없는 것 같아.변수를 취소하는 함수에서 제로로 나누는 것을 피하십시오.
참으로, 아마도 try/except ZeroDivisionError가 도움이 될 것입니다 ... –