TA的每日心情 | 衰
2019-11-19 15:29 |
|---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
MATLAB的复化柯特斯公式，提示我索引超出矩阵维度，但是我改不来
) A+ H& K% E6 V" a& t! V$ {function cotes(a,b)" V* F2 ~: i# N" u3 \6 S
n=4;" m* y" ]$ @0 N( k
for i=1:5: b0 U9 C N- \3 m! z
x=a b-a)/4*n:b;
" `9 U5 B) ?4 ]& _m=4*n+1;
# o* Y: Z3 T- c9 G. F, dh=(b-a)/n;
1 Q, B" ~% [% I is=(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));
" A# C" a' x% L1 [: G1 Kfprintf('s(%d))=%f\n',n,s);
3 b) I1 h) n" D# _; h& gn=n*4;6 k+ V0 r) h5 r" I
end
+ u! ?7 @* `% V0 l! V& c请各位大佬帮我看看？？？ |
|
/1 
发表于 2019-6-19 15:23
收藏
淘帖
支持！
反对！