TA的每日心情 | 衰 2019-11-19 15:29 |
---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
MATLAB的复化柯特斯公式,提示我索引超出矩阵维度,但是我改不来: A I* q0 @ C5 T6 k
function cotes(a,b)
# ~. Q) H2 _3 @8 T( x6 ]n=4;
! S9 ^9 L c& H( m `for i=1:51 X1 Y5 E' s% E" d6 Z( m
x=a b-a)/4*n:b;( a2 n0 o8 m( X- B( J J9 n
m=4*n+1;
9 G! Q n4 {/ I& v8 m+ V* Xh=(b-a)/n;6 [- Z- |0 O( E7 U
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));. N& d- O d8 u! w, T& I9 u; W5 l e
fprintf('s(%d))=%f\n',n,s);( a! P" y! p, N- s
n=n*4;) R' e$ a! K2 F! q3 w
end' u9 C' g- [1 v; S3 ?! F7 A: G3 V) x: g
请各位大佬帮我看看??? |
|