TA的每日心情 | 衰
2019-11-19 15:29 |
|---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
MATLAB的复化柯特斯公式，提示我索引超出矩阵维度，但是我改不来
( X+ M* U2 g2 F# F8 x6 S3 M6 xfunction cotes(a,b)
! i5 F3 u+ [) G& r- rn=4;* K. n" e( u+ D, Z# U) O) k
for i=1:5$ o& {/ ]( h7 m$ z/ j( E) ?
x=a b-a)/4*n:b;
' _$ v$ |% z2 {: M0 n/ E( km=4*n+1;2 b; B3 C3 V0 s0 j9 N# t
h=(b-a)/n;6 D% `6 s( T! R9 _8 _ n+ [
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));
6 \0 ]- m" t% B0 ]fprintf('s(%d))=%f\n',n,s);
" A6 m T* L3 @n=n*4;3 e) Y( n8 R' c3 e- ^) w
end+ W& A5 A+ S4 E) |+ b( ?+ U9 j1 r
请各位大佬帮我看看？？？ |
|
/1 
发表于 2019-6-19 15:23
收藏
淘帖
支持！
反对！