给你介绍两种方法,一种是符号运算,一种是数值运算,代码分别如下:
% 符号运算
syms x;
f=0.04/((1.1*10^14)*(exp(-11000/(323+39.01*x))*0.04*0.04*(1-x)^2))
double(int(f,0,0.85))
% 数值运算
f=@(x) 0.04./((1.1*10^14).*(exp(-11000./(323+39.01*x)).*0.04*0.04.*(1-x).^2))
quad(f,0,0.85)
两种运算结果都是一样的,都是91.1109.