24小时热门版块排行榜

返回列表

当前主题已经存档。

linwenwen-1

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 409
红花: 1
帖子: 122
在线: 55小时
虫号: 136978
注册: 2005-12-17
专业: 半导体材料

[交流] 【求助】求租一个计算式

请看附件的一个定积分，k=-20：0.1：20 用mat如何求出来？在线等，

回复此楼

» 猜你喜欢

1楼 2009-09-22 19:15:16

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

liang_2004

木虫 (正式写手)

始作俑者

应助: 16 (小学生)
金币: 1803.3
散金: 56
红花: 5
帖子: 377
在线: 124.4小时
虫号: 848455
注册: 2009-09-15
专业: 无机非金属类光电信息与功

★
linwenwen-1(金币+1):谢谢参与
nono2009(金币+0,VIP+0):谢谢参与。写出来后加BB。 9-23 06:07

这个式子直接不好写应该的变化一下
现在晚了明天试试呵呵

赞一下(7人)

回复此楼

走着，看看再说

2楼2009-09-22 23:32:07

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

change0618

铁杆木虫 (著名写手)

方丈大师

应助: 44 (小学生)
金币: 17724.5
红花: 17
帖子: 2413
在线: 546.7小时
虫号: 496517
注册: 2008-01-19
专业: 化学反应工程

★ ★ ★ ★
linwenwen-1(金币+1):谢谢参与
nono2009(金币+3,VIP+0):谢谢交流！欢迎常来。 9-23 13:11

function test
clear all
clc
format long
global k

kspan=-20:0.5:20;
y=zeros(length(kspan),1);
for i=1:length(kspan)
k=kspan(i);
y(i)=quadgk(@func,0,inf);
end
figure(1)
plot(kspan,y,'r*-')
xlabel('k')
ylabel('y=\int_0^{\infty}x^{0.5}/(1+e^{x-k})dx')

function f=func(x)
global k
f=x.^0.5./(1+exp(x-k));

赞一下(7人)

回复此楼

3楼2009-09-23 13:00:11

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