3楼: Originally posted by 月只蓝 at 2016-03-31 11:03:15
以下代码要完整复制到一个m文件中，再运行。
function qiu_jifen
clear all;clc
format long
jifen=quadgk(@objfun,3.4,10.5)

function f=objfun(t)
T=150.92*t.^(-0.184);
f=exp(-755.8335./(T+70)).*exp ...

jifen=quadgk(@(t)exp(-755.8335./((150.92*t.^(-0.184))+70)).*exp(-70000*(143+(150.92*t.^(-0.184)))./2./((150.92*t.^(-0.184)).^2)./(143-(150.92*t.^(-0.184)))),3.4,10.5)

function qiu_jifen
clear all;clc
format long
jifen=quadgk(@objfun,3.4,10.5)

function f=objfun(t)
T=150.92*t.^(-0.184);
f=exp(-755.8335./(T+70)).*exp(-70000*(143+T)./2./(T.^2)./(143-T));