from scipy.integrate import quad import numpy as np # 定義された関数 N = int(input()) # N の具体値を設定 f = lambda x: np.sqrt(x**3 + N**3) # 数値積分 result, error = quad(f, 0, 1) print(result)