TA的每日心情 | 衰 2019-11-19 15:29 |
---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
MATLAB的复化柯特斯公式,提示我索引超出矩阵维度,但是我改不来/ W8 J- U+ c: k6 ~0 P& [' Y5 M- Z3 W" x
function cotes(a,b)! h' K. a5 ]# I4 Y
n=4; D7 D9 x+ H2 g
for i=1:54 m. \" r) J, r4 v1 {
x=a b-a)/4*n:b;
% u/ i9 b$ n: V# l( C! R! N( am=4*n+1;
7 a% q, \8 j3 Yh=(b-a)/n;- L: s+ |9 I9 |. a( L8 l+ R; j) C& ?
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));
5 {1 f( w, ]) m" w$ e3 ~6 jfprintf('s(%d))=%f\n',n,s);
7 p& q# P& v s) ^+ v$ xn=n*4;( ]6 j8 c/ X% Q2 f B1 f% q2 x
end0 U: H) E% Y) h* N9 l
请各位大佬帮我看看??? |
|