A. vb 辛普森法求积分 我要具体能运行的代码
Private Sub Command1_Click()
Dim 下限 As Single, 上限 As Single, 允许误差 As Single, 实际误差 As Single, 分条 As Integer
下限 = Val(InputBox("请输入积分下限"))
上限 = Val(InputBox("请输入积分上限"))
If 下限 > 上限 Then
MsgBox "下限高于上限。怀疑输入手误。程序自动交换"
t = 下限
下限 = 上限
上限 = t
End If
允许误差 = Val(InputBox("请输入允许误差"))
Print "积分下限："; 下限
Print "积分上限："; 上限
Print "允许误差为："; 允许误差
Print "积分结果为："; 积分(下限, 上限, 允许误差, 实际误差, 分条)
Print "划分梯形个数："; 分条
Print "截断误差为："; 实际误差
End Sub
Private Function 积分(下限 As Single, 上限 As Single, 误差要求 As Single, 误差 As Single, 块数 As Integer)
Dim p As Single
块数 = 2
积分值 = 0
步长 = (上限 - 下限) / 块数
p = 函数(下限) + 函数(上限)
d = 函数(下限 + 步长)
Do
s = (p + 4 * d) * 步长 / 3
误差 = Abs(s - 积分值)
If 误差 <= 误差要求 And 块数 > 100 Then Exit Do
块数 = 块数 * 2
步长 = (上限 - 下限) / 块数
积分值 = s
p = p + 2 * d
d = 0
For i = 1 To 块数 - 1 Step 2
d = d + 函数(下限 + i * 步长)
Next i
Loop
积分 = s
End Function
Private Function 函数(自变量 As Single) '这个函数的内容是要你自己修改的
函数 = 自变量 * 自变量
End Function
B. 在MATLAB用辛普森法则计算x^2+3x在区间[1,2]上面的积分值
用辛普森法则计算，可以得到 x≈ 6.8333。实现代码
function simpson_main()
clc,close all
f=inline('x.^2+3*x','x');
simpson(1,2,10,f)
辛普森公式代码内容：
h = (b - a) / n;
x = linspace(a,b,2*n+1);
y1 = feval(f,x);
y1(2:2:2*n) = 4 * y1(2:2:2*n);
y1(3:2:2*n-1) = 2 * y1(3:2:2*n-1);
y = h / 6 * sum(y1);
C. matlab中如何用复合辛普森公式求二重积分
代码如下：
function q=DblSimpson(f,a,A,b,B,m,n)
if(m==1 && n==1) %辛普森公式
q=((B-b)*(A-a)/9)*(subs(sym(f),findsym(sym(f)),{a,b})+...
subs(sym(f),findsym(sym(f)),{a,B})+...
subs(sym(f),findsym(sym(f)),{A,b})+...
subs(sym(f),findsym(sym(f)),{A,B})+...
4*subs(sym(f),findsym(sym(f)),{(A-a)/2,b})+...
4*subs(sym(f),findsym(sym(f)),{(A-a)/2,B})+...
4*subs(sym(f),findsym(sym(f)),{a,(B-b)/2})+...
4*subs(sym(f),findsym(sym(f)),{A,(B-b)/2})+...
16*subs(sym(f),findsym(sym(f)),{(A-a)/2,(B-b)/2}));
else %复合辛普森公式
q=0;
for i=0:n-1
for j=0:m-1
x=a+2*i*(A-a)/2/n;
y=b+2*j*(B-b)/2/m;
x1=a+(2*i+1)*(A-a)/2/n;
y1=b+(2*j+1)*(B-b)/2/m;
x2=a+2*(i+1)*(A-a)/2/n;
y2=b+2*(j+1)*(B-b)/2/m;
q=q+subs(sym(f),findsym(sym(f)),{x,y})+...
subs(sym(f),findsym(sym(f)),{x,y2})+...
subs(sym(f),findsym(sym(f)),{x2,y})+...
subs(sym(f),findsym(sym(f)),{x2,y2})+...
4*subs(sym(f),findsym(sym(f)),{x,y1})+...
4*subs(sym(f),findsym(sym(f)),{x2,y1})+...
4*subs(sym(f),findsym(sym(f)),{x1,y})+...
4*subs(sym(f),findsym(sym(f)),{x1,y2})+...
16*subs(sym(f),findsym(sym(f)),{x1,y1});
end
end
end
q=((B-b)*(A-a)/36/m/n)*q;
叫我雷锋！
D. 鲍勃·辛普森的公司职位
XTOEnergy Inc.
总部地点：美国
上市地点：纽约证交所股票代码：XTO美国克洛斯提柏石油公司(XTOEnergy)
克洛斯提柏石油公司（XTOEnergyInc.，NYSE:XTO）是
美国财富500强能源生产企业，公司成立于1986年。
洛斯提柏石油公司参与石油和气体物产的收购、发展、开发和勘探及石油和天然气的生产，处理，行销和运输。公司增长主要通过被证明的石油和天然气储备的收购，及发展和开发活动和战略收购额外的已获得物产临近的权益。美国克洛斯提柏石油公司并没有真正发现新的石油和天然气资源。相反，这家位于沃思堡的公司在首席执行官鲍勃·辛普森的领导下，通过所谓的“野猫勘探”（wildcatexploitation）实现了繁荣发展。公司的产量已经连续13年实现增长。
自1993年以每股13美元的价格公开上市后，公司的股价已经拆分了7次，目前已经涨到每股51.38美元。2006年，公司的营业额和净收入分别比2005年增长了30%和61%。目前，这股增长势头还未显露出任何终止的迹象。XTO已探明的资源储量主要集中在得克萨斯、新墨西哥、阿肯色以及邻近各州，这些储量已经达到有史以来的最高值，而且全部做好了开采准备。
收购
XTOEnergy同意以41.9亿美元的价格收购HuntPetroleum
2008年6月10日天然气生产商XTOEnergyInc.(XTO)已同意以41.9亿美元现金加股票的价格收购私人持股的HuntPetroleum。该交易或将帮助平息Hunt家族内部的矛盾。
董事长兼首席执行长BobR·impson称收购的资产将对该公司业务构成天然的补充，这是决定公司未来和为股东创造价值的关键时刻。
