TA的每日心情 | 衰 2019-11-19 15:29 |
---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
MATLAB的复化柯特斯公式,提示我索引超出矩阵维度,但是我改不来9 ^2 d1 Z. I& `
function cotes(a,b)
+ Z3 S5 Y. n/ [& bn=4;
: ]' o# ^6 [2 I: @for i=1:56 ]" h% X. |; q1 Z5 O' X" o
x=a b-a)/4*n:b;' X2 B8 q% v! f; s
m=4*n+1;
; g, P5 x% d3 n" Hh=(b-a)/n;
1 X; N8 b# v: f4 y8 I' l9 ]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));
7 C* p% V8 V5 Efprintf('s(%d))=%f\n',n,s);
% G1 j! k: {6 c; v J+ _# B& Q8 In=n*4;
5 i' G: M/ d# z( s5 f3 Pend) G3 i- j8 z3 A' W
请各位大佬帮我看看??? |
|