TA的每日心情 | 衰 2019-11-19 15:29 |
|---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
MATLAB的复化柯特斯公式,提示我索引超出矩阵维度,但是我改不来& h0 F4 I k% z" c% g
function cotes(a,b)/ k' t7 ]" K( C. v/ t7 ^$ q3 s: m/ T
n=4;2 u5 {" d/ T- |+ S1 I% G1 E
for i=1:5' f2 L$ S0 G4 `
x=a b-a)/4*n:b;
( W* f3 e! x) qm=4*n+1;, x* x: ^' z8 s! h0 P) h
h=(b-a)/n;1 p# U( z3 Q& F* b7 ~' w. d: K
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));
: ~# t5 B: _6 Z, y: F2 K# Mfprintf('s(%d))=%f\n',n,s);
# s2 h# S9 d3 hn=n*4;" I& _! Z/ s# z! ? V. I a( h3 ]+ V
end
& H7 |' G) ^1 v2 ^3 H6 s请各位大佬帮我看看??? |
|