function res = ComSimpson(f,n,a,b)
n=input('please input"n"=');
a=input('please input"a"=');
b=input('please input"b"=');
B=input('please input"梁的宽度B"=');
F=input('please input"集中力F"=');
E=input('please input"弹性模量E"=');
format long;
if b < a
c = b;
b = a;
a = c;
end
h = (b-a)/n;
d = f(a);
for i = a+h:h:b-h
d = d + (2 * f(i));
end
for i = a+h/2:h:b-h/2
d = d + (4 * f(i));
end
d=d+f(b);
res=(d * h / 6);
res=F*res/B/E;
function f=f(x)
f=1/(-0.5*x+1);
end
ComSimpson
please input"n"=20
please input"a"=800
please input"b"=0
please input"B"=10
please input"F"=300
please input"E"=50000
输入参数的数目不足。
出错 ComSimpson (第 16 行)
d = f(a);