TA的每日心情 | 衰 2019-11-19 15:29 |
---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
MATLAB的复化柯特斯公式,提示我索引超出矩阵维度,但是我改不来$ N8 U) v& R: d- x
function cotes(a,b)
7 k+ _/ M- A6 ^5 y' \( i- ~ ~# X' Vn=4;! k# j2 U2 G( Y" D
for i=1:5
" S4 l% b7 F. a! ]' @: tx=a b-a)/4*n:b; O" [7 x5 }1 G& a3 ?; i
m=4*n+1;
8 `4 d7 p5 x# Hh=(b-a)/n;
9 z6 H8 O. g6 Ks=(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));
& g, G. a- B# L2 S" afprintf('s(%d))=%f\n',n,s);
% @2 A+ e G a5 k% H) Vn=n*4;
: D1 c3 `5 ]: h: r" [ g8 ~end
; |6 w e: q7 Y3 k- N请各位大佬帮我看看??? |
|