TA的每日心情 | 衰 2019-11-19 15:29 |
|---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
MATLAB的复化柯特斯公式,提示我索引超出矩阵维度,但是我改不来+ ]/ I/ ]0 U! q5 u
function cotes(a,b)9 `( J& e" Y- ~9 q$ s& {- X7 H% T+ |
n=4;
f3 V6 Y+ F- \, efor i=1:5
" u$ _* p7 T% ?9 ?( Ix=a b-a)/4*n:b;8 S7 \+ P# p& d3 R" h
m=4*n+1;4 x9 X( Q* Q& T7 w8 p4 \. l( n
h=(b-a)/n;9 P5 r" m, p5 j5 n* B1 K; W. @
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));, w7 L2 p1 G& {% T
fprintf('s(%d))=%f\n',n,s);3 m4 a: A+ C" p- d
n=n*4;
6 x& T6 U7 z/ pend
+ M5 a9 D1 x* |3 w请各位大佬帮我看看??? |
|