TA的每日心情 | 衰 2019-11-19 15:29 |
|---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
MATLAB的复化柯特斯公式,提示我索引超出矩阵维度,但是我改不来
2 @, t) ?9 j+ J% vfunction cotes(a,b)9 \2 d2 o" H2 {, @+ n9 L" |2 X: W
n=4;9 c4 x( T0 }9 d, K2 e' J% N3 @
for i=1:5
3 O+ y0 l3 ~4 J) O) D6 w3 {x=a b-a)/4*n:b;- n! V, f! f r
m=4*n+1;2 W8 f8 f# }, J
h=(b-a)/n;9 q: C% S) c$ r6 }3 t
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));
( w0 s7 W" d! r9 n% lfprintf('s(%d))=%f\n',n,s);
) {- q! T' t5 `1 Jn=n*4;
% L. L& n( x- ` |* a: Yend1 _, r4 m4 H4 F, q
请各位大佬帮我看看??? |
|