here is my code:
这是我的代码:
C=@(k) [k,k,2.*k;3,2.*k,5;1,k,k];
AV=@(k,t) [3*t, 6, 9]*C(k)*[3*t ;6 ;9];
avaint=@(k,a,b) quadgk(@(k) AV(k,t),a,b);
AVAR=@(t) avaint(t,0,87600);
C=@(k) [k,k,2.*k;3here is my code:
这是我的代码:
C=@(k) [k,k,2.*k;3,2.*k,5;1,k,k];
AV=@(k,t) [3*t, 6, 9]*C(k)*[3*t ;6 ;9];
avaint=@(k,a,b) quadgk(@(k) AV(k,t),a,b);
AVAR=@(t) avaint(t,0,87600);
C=@(k) [k,k,2.*k;3