from scipy.integrate import quad import numpy as np N=int(input()) f = lambda x: np.sqrt(x**3 + N**3) # 数値積分 result, error = quad(f, 0, 1) print(result)