積分変数をベクトルと​して与える際の重積分​を行う際の関数とIn​tegral2の扱い​方

のように与えた場合の，x1, x2についての積分になります．

fun =  @(x, y) [x y] * [1 2; 3 4] * [x; y]
q = integral2(fun, 0, 1, 0, 1)　%積分範囲はa=c=0, b=d=1として記述

を行うと，以下のようなエラーになります．

エラー:  * 
内部行列の次元は一致しなければなりません。

エラー: LiveEditorEvaluationHelperESectionEval>@(x1,x2)[x1,x2]*[1,2;3,4]*[x1;x2]

エラー: integral2Calc>integral2t/tensor (line 228)
        Z = FUN(X,Y);  NFE = NFE + 1;

エラー: integral2Calc>integral2t (line 55)
[Qsub,esub] = tensor(thetaL,thetaR,phiB,phiT);

エラー: integral2Calc (line 9)
    [q,errbnd] = integral2t(fun,xmin,xmax,ymin,ymax,optionstruct);

エラー: integral2 (line 106)
    Q = integral2Calc(fun,xmin,xmax,yminfun,ymaxfun,opstruct);

式を展開し，

fun =  @(x1, x2) x1.^2 + 4 .*x2 .^2 + 5 .* x1 .* x2
q = integral2(fun, 0, 1, 0, 1)

>>
ans =
   2.916666666666940

のようにすれば問題なく積分が可能なのですが，より複雑な式の積分を行いたいので，どのように関数を定義してあげればベクトルのまま積分関数（integralなど）に入力することが可能か教えていただきたいです．

Answer 1

[Hernia Baby]

integral2は行列の非積分関数には適用できないみたいです。

integralであれば'ArrayValue'オプションをtrueにすることでできそうです。

参考：using integral2 on a matrix with function handler components

ということでintegralでやってみました。

clc,clear;
fun =  @(x, y) [x y] * [1 2; 3 4] * [x; y];
q = integral(@(y) integral(@(x) fun(x,y),0,1,'ArrayValued',true),0,1,'ArrayValued',true);

format long
q

q = 
   2.916666666666667

[Mutsumi Yoshino]

過去の質問を追いきれていませんでした．教えていただきありがとうございました．

また，解決策も提示していただき，ありがとうございました．