版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

wxy尘埃

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 518.6
散金: 20
红花: 3
帖子: 117
在线: 57小时
虫号: 1140882
注册: 2010-11-07
性别: MM
专业: 极地科学

[求助] matlab积分求解

对于下列积分式，如何用matlab来求解，不胜感激！
x.^4*exp(x)/(exp(x)-1).^2 对x积分

回复此楼

» 猜你喜欢

40+中年人失业后，还能找到好工作吗？已经有15人回复
二审4个多月了还是Awaiting Reviewer Scores 已经有5人回复
申请科研项目纯靠实力能成功吗？已经有12人回复
用NBS溴化时溶剂除了CCl4外还有什么已经有6人回复
极限拉扯，如何抉择已经有15人回复
OpenAI又攻克数学难题了，离AI做有机合成还有多远？已经有13人回复
课题组会帮你通过考核么？已经有7人回复
招收博士生一名已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

微笑吧，为你的曾经拥有~~~

1楼 2012-05-01 20:39:46

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

lidong

金虫 (小有名气)

光明骑士

应助: 1 (幼儿园)
金币: 2569
散金: 20
帖子: 231
在线: 70.1小时
虫号: 1204940
注册: 2011-02-17
性别: GG
专业: 固体力学

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
wxy尘埃: 金币+10, ★★★★★最佳答案 2012-05-04 20:07:08

%龙贝格外推法求函数积分，结果矩阵左下角数值为精度最高的求积结果
function bgl
a=input('a=');          %输入积分下限
b=input('b=');          %输入积分上限
n=input('n=');             %输入外推次数（n=3~5即可满足一般精度要求）
digits(10);
h=vpa(b-a);
y1=f(a);
y2=f(b);
t(1,1)=vpa(h/2*(y1+y2));
for i=2:n
s=0;
for k=1:2^(i-2)
      y0=f(a+2*(2*k-1)*(b-a)/2^i);
      s=s+y0;
end
t(1,i)=vpa(t(1,(i-1))/2+h/2*s);
h=h/2;
end
for m=2:n
for j=1

n+1-m)
c=4^(m-1);
t(m,j)=vpa(c/(c-1)*t((m-1),(j+1))-1/(c-1)*t((m-1),j));
end
end
t
end

function y=f(x)
if x==0
y=0;
else
y=vpa(x.^4*exp(x)/(exp(x)-1).^2);
end
end

赞一下

回复此楼

条条大路通罗马，总有一条是属于我的

5楼2012-05-02 16:39:56

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 8 个回答

lurencyj

木虫 (著名写手)

应助: 159 (高中生)
金币: 2869.2
散金: 520
红花: 8
沙发: 10
帖子: 1244
在线: 148.3小时
虫号: 888093
注册: 2009-10-29
性别: GG
专业: 凝聚态物性I:结构、力学和

【答案】应助回帖

★
感谢参与，应助指数 +1
xzhdty: 金币+1, 专家考核, 谢谢了 2012-05-02 06:42:47

工具而已，何必执着.
Mathematica:

CODE:

Integrate[x^4 Exp[x]/(Exp[x] - 1)^2, x]

-x^4 - x^4/(-1 + e^x) + 4 x^3 Log[1 - e^x] + 12 x^2 PolyLog[2, e^x] - 

 24 x PolyLog[3, e^x] + 24 PolyLog[4, e^x]

[ 发自手机版 http://muchong.com/3g ]

赞一下

回复此楼

很女子很弓虽大

2楼2012-05-01 22:36:50

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖