sympy.lambdify(var, func, "numpy")はfuncが定数関数だと配列を返してくれないので： def numpi(var, func): if func.is_constant(): return lambda var: np.full_like(var, func) else: return sp.lambdify(var, func, "numpy")