TA的每日心情 | 衰 2019-11-19 15:29 |
---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
MATLAB的复化柯特斯公式,提示我索引超出矩阵维度,但是我改不来
6 Q9 s2 }: Z1 p% O# U3 q$ Mfunction cotes(a,b)! l4 A6 p6 Q( r$ q4 I) m8 u7 ]
n=4;
! U" K. d$ Y- a& \, z5 ]for i=1:58 m" _5 S% t# L7 Z
x=a b-a)/4*n:b;) {9 O; _/ |' V& j2 J3 ^! N
m=4*n+1;
: X7 M# z) \: ^/ f+ eh=(b-a)/n;
8 s5 y: L! Y" B! N' Us=(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/ {. @% r; Q, T/ S
fprintf('s(%d))=%f\n',n,s);7 a! n) v2 ~# I8 h
n=n*4;) L6 o! ?" g7 z5 e
end
) ]( f+ k0 n- r* J9 k请各位大佬帮我看看??? |
|