一、题目 二、代码 from scipy.integrate import quadimport numpy as np# 函数f(u)=1/2*u**2,故f`(u)=udef f_positive(upp_value):# 积分f+(u)中需要使用的函数is_over_zero = int(upp_value > 0)return is_over_zero * upp_valuedef