6楼: Originally posted by daijiantao at 2016-04-10 23:15:42
如果我只有积分上限变化（tH），其余都不变，能基于上面的程序帮忙给一个程序吗，最好可以出图？...

function qiu_jifen
clear all;clc
format long
tH_0=100;    %积分上限tH的最小值
tH_n=500;    %积分上限tH的最大值
num=tH_n-tH_0+1;     %积分上限tH总的变化数目

tH_data=linspace(tH_0,tH_n,num);
for i=1:length(tH_data)
    tH=tH_data(i);
    dc(i)=quadgk(@obj,0,tH);
end
[tH_data' dc']
figure(1)
plot(tH_data,dc,'b-')

function f=obj(t)
M=0.230873094;
N=0.025914445;
n=0.86;
exp=114.8;
wf0=46;
p1=M;
p2=N;
p3=n;
p4=exp;
p5=wf0;
f=(((p1./(1+t)-p2.*log10(1+t)./(1+t)).^p3./p4)-(p1./(1+t)-p2.*log10(1+t)./(1+t))./p5);

2楼: Originally posted by 月只蓝 at 2016-04-07 13:51:21
“数值解分别计算一个一个案例”，对于不同的案例，积分式什么参数会改变？