UC知道下面这个mathematica程序 哪位高手能改成 matlab 程序
来源：百度知道 编辑：UC知道 时间：2024/06/30 14:06:09
Clear[a,b,x,n,f];
a=Input["a="]
b=Input["b="]
f[x_]=3/20*Exp[-9(x-170)^2/800]/Sqrt[2Pi]
eps=0.0000001;
n=1;
h=b-a;
t1=(f[a]+f[b])*h/2;
h=h/2;
t2=t1/2+h*f[a+h];
er=t2-t1//N;
While[Abs[er]>eps,
Print["n=",2^n,"定积分值为",N[t2,10]];
Print["误差=",er];
h=h/2;
t1=t2;
n =n+1;
t2=t1/2+h*Sum[f[a+k*h],{k,1,2^n,2}];
er=t2-t1//N];
Print["n=",2^n,"定积分值为",N[t2,10]];
Print["误差=",er]
a=Input["a="]
b=Input["b="]
f[x_]=3/20*Exp[-9(x-170)^2/800]/Sqrt[2Pi]
eps=0.0000001;
n=1;
h=b-a;
t1=(f[a]+f[b])*h/2;
h=h/2;
t2=t1/2+h*f[a+h];
er=t2-t1//N;
While[Abs[er]>eps,
Print["n=",2^n,"定积分值为",N[t2,10]];
Print["误差=",er];
h=h/2;
t1=t2;
n =n+1;
t2=t1/2+h*Sum[f[a+k*h],{k,1,2^n,2}];
er=t2-t1//N];
Print["n=",2^n,"定积分值为",N[t2,10]];
Print["误差=",er]
mathematica程序有问题，运行结果狗屁不是.现在还不大看的懂
帮你纠正几个错误：
Clear[a, b, x, n, f];
a = Input["a="]
b = Input["b="]
f = 3/20*Exp[-9 (# - 170)^2/800]/Sqrt[2 Pi] &
eps = 0.0000001;
h = b - a;
t1 = (f[a] + f[b])*(b - a)/2;
h = h/2;
t2 = t1/2 + h*f[a + h];
er = t2 - t1 // N;
While[Abs[er] > eps, Print["n=", 2^n, "定积分值为", N[t2, 10]];
Print["误差=", er];
h = h/2;
t1 = t2;
n++;
t2 = t1/2 + h*Sum[f[a + k*h], {k, 1, 2^n, 2}];
er = t2 - t1 // N];
Print["n=", 2^n, "定积分值为", N[t2, 10]];
Print["误差=", er]