由于x=0时,K1为无穷大(inf),y也为无穷大,在使用trapz函数时就认为b是无效值(nan)。解决的方法是人为把x的0值改为eps值(eps系统默认为无穷小量)。即 x=(eps:0.0001:1)。运行结果值,b=0.42596