TA的每日心情 | 衰 2019-11-19 15:29 |
---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
MATLAB的复化柯特斯公式,提示我索引超出矩阵维度,但是我改不来 O4 C* |( a4 ^$ `
function cotes(a,b)+ B$ g# D/ r* E- l
n=4; m; R% }3 {( U, r- J m9 ]
for i=1:5
( I* A P- ?" E8 I& e, [! k' t( Ux=a b-a)/4*n:b;: @/ o: j# W; o9 }
m=4*n+1;3 v3 g. Q& D0 M- G
h=(b-a)/n;5 t4 s0 C ]6 u! n
s=(h/90)*(7*f1(a)+32*sum(f1(x(1:4:m-4)))+12*sum(f1(x(0:4:m-5)))+32*sum(f1(x(3:4:m-2)))+14*sum(f1(x(4:4:m-5)))+7*f1(b));
6 z0 i; }8 k- Q- n* Cfprintf('s(%d))=%f\n',n,s);
- q9 A& A' [ w* Fn=n*4;4 _0 E9 E4 m1 D' C) Z% H9 _. z. \$ g
end1 u( C! |& V, ^9 ?, [( x1 h% q9 t
请各位大佬帮我看看??? |
|