MATLAB Answers

Double integral on matlabFunction

2 ビュー (過去 30 日間)
vinayak vijay
vinayak vijay 2020 年 10 月 20 日
コメント済み: vinayak vijay 2020 年 10 月 20 日
I am trying to integrate a function(Ke) with arguments (x,y). The code is given below:
clear all;
syms x y a0 a1 a2 a3 a4 a5 a6 a7 a8 a9 a10 a11 de1 de2 de3 de4 de5 de6 de7 de8 de9 de10 de11 de12
w = @(x,y)a0 +a1*x + a2*y + a3*x^2 + a4*x*y + a5*y^2 + a6*x^3 + a7*(x^2)*y + a8*x*y^2 +a9*y^3 + a10*(x^3)*y + a11*x*(y^3);
a=[a0 a1 a2 a3 a4 a5 a6 a7 a8 a9 a10 a11].';
A =[1 0 0 0 0 0 0 0 0 0 0 0;0 1 0 0 0 0 0 0 0 0 0 0;0 0 1 0 0 0 0 0 0 0 0 0;
1 0.5 0 0.25 0 0 0.125 0 0 0 0 0;0 1 0 1 0 0 0.75 0 0 0 0 0;0 0 1 0 0.5 0 0 0.25 0 0 0.125 0;
1 0.5 0.5 0.25 0.25 0.25 0.125 0.125 0.125 0.125 0.0625 0.0625;
0 1 0 1 0.5 0 0.75 0.5 0.25 0 0.375 0.125;
0 0 1 0 0.5 1 0 0.25 0.5 0.75 0.125 0.375;
1 0 0.5 0 0 0.25 0 0 0 0.125 0 0;
0 1 0 0 0.5 0 0 0 0.25 0 0 0.125;
0 0 1 0 0 1 0 0 0 0.75 0 0];
de = [de1 de2 de3 de4 de5 de6 de7 de8 de9 de10 de11 de12].';
H = ([1 x y x^2 x*y y^2 x^3 (x^2)*y x*(y^2) y^3 (x^3)*y x*(y^3)]*inv(A)).';
B = [diff(H(1),x,2) diff(H(2),x,2) diff(H(3),x,2) diff(H(4),x,2) diff(H(5),x,2) diff(H(6),x,2) diff(H(7),x,2) diff(H(8),x,2) diff(H(9),x,2) diff(H(10),x,2) diff(H(11),x,2) diff(H(12),x,2);
diff(H(1),y,2) diff(H(2),y,2) diff(H(3),y,2) diff(H(4),y,2) diff(H(5),y,2) diff(H(6),y,2) diff(H(7),y,2) diff(H(8),y,2) diff(H(9),y,2) diff(H(10),y,2) diff(H(11),y,2) diff(H(12),y,2);
diff(diff(H(1),x),y) diff(diff(H(2),x),y) diff(diff(H(3),x),y) diff(diff(H(4),x),y) diff(diff(H(5),x),y) diff(diff(H(6),x),y) diff(diff(H(7),x),y) diff(diff(H(8),x),y) diff(diff(H(9),x),y) diff(diff(H(10),x),y) diff(diff(H(11),x),y) diff(diff(H(12),x),y)];
E = 210e9;
h = 0.001;
nu = 0.3;
d = (E*h^3)/12*(1-nu^2);
D = [d nu*d 0; nu*d d 0; 0 0 0.5*(1-nu)*d];
Ke = matlabFunction(B.'*D*B,'Vars',[x y]);
K_e = integral2(Ke,0,0.5,0,0.5)
The error message is:
Error using reshape
Number of elements must not change. Use [] as
one of the size inputs to automatically
calculate the appropriate size for that
dimension.
Error in
symengine>@(x,y)reshape([(x.*4.8e+1+y.*9.6e+1-x.*y.*1.92e+2-2.4e+1).*(x.*1.22304e+3+y.*1.75812e+3-x.*y.*3.97488e+3-4.9686e+2)+(x.*9.6e+1+y.*4.8e+1-x.*y.*1.92e+2-2.4e+1).*(x.*1.75812e+3+y.*1.22304e+3-x.*y.*3.97488e+3-4.9686e+2)+(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),(x.*9.6e+1+y.*4.8e+1-x.*y.*1.92e+2-2.4e+1).*(x.*3.822e+2+y.*2.548e+2-x.*y.*7.644e+2-6.37e+2./5.0)+(x.*4.8e+1+y.*9.6e+1-x.*y.*1.92e+2-2.4e+1).*(x.*1.1466e+2+y.*7.644e+1-x.*y.*2.2932e+2-3.822e+1)+(x.*(-8.918e+1)+x.^2.*1.3377e+2+1.11475e+1).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0),(x.*4.8e+1+y.*9.6e+1-x.*y.*1.92e+2-2.4e+1).*(x.*2.548e+2+y.*3.822e+2-x.*y.*7.644e+2-6.37e+2./5.0)+(x.*9.6e+1+y.*4.8e+1-x.*y.*1.92e+2-2.4e+1).*(x.*7.644e+1+y.*1.1466e+2-x.*y.*2.2932e+2-3.822e+1)+(y.*(-8.918e+1)+y.^2.*1.3377e+2+1.11475e+1).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0),-(x.*4.8e+1+y.*9.6e+1-x.*y.*1.92e+2-2.4e+1).*(x.*1.22304e+3+y.*2.2932e+2-x.*y.*3.97488e+3-1.1466e+2)-(x.*9.6e+1+y.*4.8e+1-x.*y.*1.92e+2-2.4e+1).*(x.*1.75812e+3+y.*7.644e+2-x.*y.*3.97488e+3-3.822e+2)-(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),(x.*9.6e+1+y.*4.8e+1-x.*y.*1.92e+2-2.4e+1).*(x.*3.822e+2+y.*(6.37e+2./5.0)-x.*y.*7.644e+2-6.37e+2./1.0e+1)+(x.*4.8e+1+y.*9.6e+1-x.*y.*1.92e+2-2.4e+1).*(x.*1.1466e+2+y.*3.822e+1-x.*y.*2.2932e+2-1.911e+1)-(x.*4.459e+1-x.^2.*1.3377e+2).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0),-(x.*2.548e+2-x.*y.*7.644e+2).*(x.*4.8e+1+y.*9.6e+1-x.*y.*1.92e+2-2.4e+1)-(x.*7.644e+1-x.*y.*2.2932e+2).*(x.*9.6e+1+y.*4.8e+1-x.*y.*1.92e+2-2.4e+1)-(y.*(-8.918e+1)+y.^2.*1.3377e+2+1.11475e+1).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0),(x.*7.644e+2+y.*2.2932e+2-x.*y.*3.97488e+3).*(x.*4.8e+1+y.*9.6e+1-x.*y.*1.92e+2-2.4e+1)+(x.*2.2932e+2+y.*7.644e+2-x.*y.*3.97488e+3).*(x.*9.6e+1+y.*4.8e+1-x.*y.*1.92e+2-2.4e+1)+(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),-(y.*(6.37e+2./5.0)-x.*y.*7.644e+2).*(x.*9.6e+1+y.*4.8e+1-x.*y.*1.92e+2-2.4e+1)-(y.*3.822e+1-x.*y.*2.2932e+2).*(x.*4.8e+1+y.*9.6e+1-x.*y.*1.92e+2-2.4e+1)+(x.*4.459e+1-x.^2.*1.3377e+2).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0),-(x.*(6.37e+2./5.0)-x.*y.*7.644e+2).*(x.*4.8e+1+y.*9.6e+1-x.*y.*1.92e+2-2.4e+1)-(x.*3.822e+1-x.*y.*2.2932e+2).*(x.*9.6e+1+y.*4.8e+1-x.*y.*1.92e+2-2.4e+1)+(y.*4.459e+1-y.^2.*1.3377e+2).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0),-(x.*9.6e+1+y.*4.8e+1-x.*y.*1.92e+2-2.4e+1).*(x.*2.2932e+2+y.*1.22304e+3-x.*y.*3.97488e+3-1.1466e+2)-(x.*4.8e+1+y.*9.6e+1-x.*y.*1.92e+2-2.4e+1).*(x.*7.644e+2+y.*1.75812e+3-x.*y.*3.97488e+3-3.822e+2)-(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),-(y.*2.548e+2-x.*y.*7.644e+2).*(x.*9.6e+1+y.*4.8e+1-x.*y.*1.92e+2-2.4e+1)-(y.*7.644e+1-x.*y.*2.2932e+2).*(x.*4.8e+1+y.*9.6e+1-x.*y.*1.92e+2-2.4e+1)-(x.*(-8.918e+1)+x.^2.*1.3377e+2+1.11475e+1).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0),(x.*4.8e+1+y.*9.6e+1-x.*y.*1.92e+2-2.4e+1).*(x.*(6.37e+2./5.0)+y.*3.822e+2-x.*y.*7.644e+2-6.37e+2./1.0e+1)+(x.*9.6e+1+y.*4.8e+1-x.*y.*1.92e+2-2.4e+1).*(x.*3.822e+1+y.*1.1466e+2-x.*y.*2.2932e+2-1.911e+1)-(y.*4.459e+1-y.^2.*1.3377e+2).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0),(x.*2.4e+1+y.*1.6e+1-x.*y.*4.8e+1-8.0).*(x.*1.75812e+3+y.*1.22304e+3-x.*y.*3.97488e+3-4.9686e+2)+(x.*-1.6e+1+x.^2.*2.4e+1+2.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),(x.*2.4e+1+y.*1.6e+1-x.*y.*4.8e+1-8.0).*(x.*3.822e+2+y.*2.548e+2-x.*y.*7.644e+2-6.37e+2./5.0)+(x.*-1.6e+1+x.^2.*2.4e+1+2.0).*(x.*(-8.918e+1)+x.^2.*1.3377e+2+1.11475e+1),(x.*2.4e+1+y.*1.6e+1-x.*y.*4.8e+1-8.0).*(x.*7.644e+1+y.*1.1466e+2-x.*y.*2.2932e+2-3.822e+1)+(x.*-1.6e+1+x.^2.*2.4e+1+2.0).*(y.*(-8.918e+1)+y.^2.*1.3377e+2+1.11475e+1),-(x.*2.4e+1+y.*1.6e+1-x.*y.*4.8e+1-8.0).*(x.*1.75812e+3+y.*7.644e+2-x.*y.*3.97488e+3-3.822e+2)-(x.*-1.6e+1+x.^2.*2.4e+1+2.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),(x.*2.4e+1+y.*1.6e+1-x.*y.*4.8e+1-8.0).*(x.*3.822e+2+y.*(6.37e+2./5.0)-x.*y.*7.644e+2-6.37e+2./1.0e+1)-(x.*4.459e+1-x.^2.*1.3377e+2).*(x.*-1.6e+1+x.^2.*2.4e+1+2.0),-(x.*-1.6e+1+x.^2.*2.4e+1+2.0).*(y.*(-8.918e+1)+y.^2.*1.3377e+2+1.11475e+1)-(x.*7.644e+1-x.*y.*2.2932e+2).*(x.*2.4e+1+y.*1.6e+1-x.*y.*4.8e+1-8.0),(x.*2.2932e+2+y.*7.644e+2-x.*y.*3.97488e+3).*(x.*2.4e+1+y.*1.6e+1-x.*y.*4.8e+1-8.0)+(x.*-1.6e+1+x.^2.*2.4e+1+2.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),(x.*4.459e+1-x.^2.*1.3377e+2).*(x.*-1.6e+1+x.^2.*2.4e+1+2.0)-(y.*(6.37e+2./5.0)-x.*y.*7.644e+2).*(x.*2.4e+1+y.*1.6e+1-x.*y.*4.8e+1-8.0),(y.*4.459e+1-y.^2.*1.3377e+2).*(x.*-1.6e+1+x.^2.*2.4e+1+2.0)-(x.*3.822e+1-x.*y.*2.2932e+2).*(x.*2.4e+1+y.*1.6e+1-x.*y.*4.8e+1-8.0),-(x.*2.4e+1+y.*1.6e+1-x.*y.*4.8e+1-8.0).*(x.*2.2932e+2+y.*1.22304e+3-x.*y.*3.97488e+3-1.1466e+2)-(x.*-1.6e+1+x.^2.*2.4e+1+2.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),-(x.*-1.6e+1+x.^2.*2.4e+1+2.0).*(x.*(-8.918e+1)+x.^2.*1.3377e+2+1.11475e+1)-(y.*2.548e+2-x.*y.*7.644e+2).*(x.*2.4e+1+y.*1.6e+1-x.*y.*4.8e+1-8.0),(x.*2.4e+1+y.*1.6e+1-x.*y.*4.8e+1-8.0).*(x.*3.822e+1+y.*1.1466e+2-x.*y.*2.2932e+2-1.911e+1)-(y.*4.459e+1-y.^2.*1.3377e+2).*(x.*-1.6e+1+x.^2.*2.4e+1+2.0),(x.*1.6e+1+y.*2.4e+1-x.*y.*4.8e+1-8.0).*(x.*1.22304e+3+y.*1.75812e+3-x.*y.*3.97488e+3-4.9686e+2)+(y.*-1.6e+1+y.^2.*2.4e+1+2.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),(x.*1.6e+1+y.*2.4e+1-x.*y.*4.8e+1-8.0).*(x.*1.1466e+2+y.*7.644e+1-x.*y.*2.2932e+2-3.822e+1)+(x.*(-8.918e+1)+x.^2.*1.3377e+2+1.11475e+1).*(y.*-1.6e+1+y.^2.*2.4e+1+2.0),(x.*1.6e+1+y.*2.4e+1-x.*y.*4.8e+1-8.0).*(x.*2.548e+2+y.*3.822e+2-x.*y.*7.644e+2-6.37e+2./5.0)+(y.*-1.6e+1+y.^2.*2.4e+1+2.0).*(y.*(-8.918e+1)+y.^2.*1.3377e+2+1.11475e+1),-(x.*1.6e+1+y.*2.4e+1-x.*y.*4.8e+1-8.0).*(x.*1.22304e+3+y.*2.2932e+2-x.*y.*3.97488e+3-1.1466e+2)-(y.*-1.6e+1+y.^2.*2.4e+1+2.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),(x.*1.6e+1+y.*2.4e+1-x.*y.*4.8e+1-8.0).*(x.*1.1466e+2+y.*3.822e+1-x.*y.*2.2932e+2-1.911e+1)-(x.*4.459e+1-x.^2.*1.3377e+2).*(y.*-1.6e+1+y.^2.*2.4e+1+2.0),-(y.*-1.6e+1+y.^2.*2.4e+1+2.0).*(y.*(-8.918e+1)+y.^2.*1.3377e+2+1.11475e+1)-(x.*2.548e+2-x.*y.*7.644e+2).*(x.*1.6e+1+y.*2.4e+1-x.*y.*4.8e+1-8.0),(x.*7.644e+2+y.*2.2932e+2-x.*y.*3.97488e+3).*(x.*1.6e+1+y.*2.4e+1-x.*y.*4.8e+1-8.0)+(y.*-1.6e+1+y.^2.*2.4e+1+2.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),(x.*4.459e+1-x.^2.*1.3377e+2).*(y.*-1.6e+1+y.^2.*2.4e+1+2.0)-(y.*3.822e+1-x.*y.*2.2932e+2).*(x.*1.6e+1+y.*2.4e+1-x.*y.*4.8e+1-8.0),(y.*4.459e+1-y.^2.*1.3377e+2).*(y.*-1.6e+1+y.^2.*2.4e+1+2.0)-(x.*(6.37e+2./5.0)-x.*y.*7.644e+2).*(x.*1.6e+1+y.*2.4e+1-x.*y.*4.8e+1-8.0),-(x.*1.6e+1+y.*2.4e+1-x.*y.*4.8e+1-8.0).*(x.*7.644e+2+y.*1.75812e+3-x.*y.*3.97488e+3-3.822e+2)-(y.*-1.6e+1+y.^2.*2.4e+1+2.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),-(x.*(-8.918e+1)+x.^2.*1.3377e+2+1.11475e+1).*(y.*-1.6e+1+y.^2.*2.4e+1+2.0)-(y.*7.644e+1-x.*y.*2.2932e+2).*(x.*1.6e+1+y.*2.4e+1-x.*y.*4.8e+1-8.0),(x.*1.6e+1+y.*2.4e+1-x.*y.*4.8e+1-8.0).*(x.*(6.37e+2./5.0)+y.*3.822e+2-x.*y.*7.644e+2-6.37e+2./1.0e+1)-(y.*4.459e+1-y.^2.*1.3377e+2).*(y.*-1.6e+1+y.^2.*2.4e+1+2.0),-(x.*9.6e+1+y.*4.8e+1-x.*y.*1.92e+2-2.4e+1).*(x.*1.75812e+3+y.*1.22304e+3-x.*y.*3.97488e+3-4.9686e+2)-(x.*4.8e+1-x.*y.*1.92e+2).*(x.*1.22304e+3+y.*1.75812e+3-x.*y.*3.97488e+3-4.9686e+2)-(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),-(x.*9.6e+1+y.*4.8e+1-x.*y.*1.92e+2-2.4e+1).*(x.*3.822e+2+y.*2.548e+2-x.*y.*7.644e+2-6.37e+2./5.0)-(x.*4.8e+1-x.*y.*1.92e+2).*(x.*1.1466e+2+y.*7.644e+1-x.*y.*2.2932e+2-3.822e+1)-(x.*(-8.918e+1)+x.^2.*1.3377e+2+1.11475e+1).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0),-(x.*9.6e+1+y.*4.8e+1-x.*y.*1.92e+2-2.4e+1).*(x.*7.644e+1+y.*1.1466e+2-x.*y.*2.2932e+2-3.822e+1)-(x.*4.8e+1-x.*y.*1.92e+2).*(x.*2.548e+2+y.*3.822e+2-x.*y.*7.644e+2-6.37e+2./5.0)-(y.*(-8.918e+1)+y.^2.*1.3377e+2+1.11475e+1).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0),(x.*9.6e+1+y.*4.8e+1-x.*y.*1.92e+2-2.4e+1).*(x.*1.75812e+3+y.*7.644e+2-x.*y.*3.97488e+3-3.822e+2)+(x.*4.8e+1-x.*y.*1.92e+2).*(x.*1.22304e+3+y.*2.2932e+2-x.*y.*3.97488e+3-1.1466e+2)+(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),-(x.*9.6e+1+y.*4.8e+1-x.*y.*1.92e+2-2.4e+1).*(x.*3.822e+2+y.*(6.37e+2./5.0)-x.*y.*7.644e+2-6.37e+2./1.0e+1)-(x.*4.8e+1-x.*y.*1.92e+2).*(x.*1.1466e+2+y.*3.822e+1-x.*y.*2.2932e+2-1.911e+1)+(x.*4.459e+1-x.^2.*1.3377e+2).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0),(x.*7.644e+1-x.*y.*2.2932e+2).*(x.*9.6e+1+y.*4.8e+1-x.*y.*1.92e+2-2.4e+1)+(y.*(-8.918e+1)+y.^2.*1.3377e+2+1.11475e+1).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0)+(x.*4.8e+1-x.*y.*1.92e+2).*(x.*2.548e+2-x.*y.*7.644e+2),-(x.*2.2932e+2+y.*7.644e+2-x.*y.*3.97488e+3).*(x.*9.6e+1+y.*4.8e+1-x.*y.*1.92e+2-2.4e+1)-(x.*4.8e+1-x.*y.*1.92e+2).*(x.*7.644e+2+y.*2.2932e+2-x.*y.*3.97488e+3)-(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),(y.*(6.37e+2./5.0)-x.*y.*7.644e+2).*(x.*9.6e+1+y.*4.8e+1-x.*y.*1.92e+2-2.4e+1)-(x.*4.459e+1-x.^2.*1.3377e+2).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0)+(x.*4.8e+1-x.*y.*1.92e+2).*(y.*3.822e+1-x.*y.*2.2932e+2),(x.*3.822e+1-x.*y.*2.2932e+2).*(x.*9.6e+1+y.*4.8e+1-x.*y.*1.92e+2-2.4e+1)-(y.*4.459e+1-y.^2.*1.3377e+2).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0)+(x.*4.8e+1-x.*y.*1.92e+2).*(x.*(6.37e+2./5.0)-x.*y.*7.644e+2),(x.*9.6e+1+y.*4.8e+1-x.*y.*1.92e+2-2.4e+1).*(x.*2.2932e+2+y.*1.22304e+3-x.*y.*3.97488e+3-1.1466e+2)+(x.*4.8e+1-x.*y.*1.92e+2).*(x.*7.644e+2+y.*1.75812e+3-x.*y.*3.97488e+3-3.822e+2)+(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),(y.*2.548e+2-x.*y.*7.644e+2).*(x.*9.6e+1+y.*4.8e+1-x.*y.*1.92e+2-2.4e+1)+(x.*(-8.918e+1)+x.^2.*1.3377e+2+1.11475e+1).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0)+(x.*4.8e+1-x.*y.*1.92e+2).*(y.*7.644e+1-x.*y.*2.2932e+2),-(x.*9.6e+1+y.*4.8e+1-x.*y.*1.92e+2-2.4e+1).*(x.*3.822e+1+y.*1.1466e+2-x.*y.*2.2932e+2-1.911e+1)-(x.*4.8e+1-x.*y.*1.92e+2).*(x.*(6.37e+2./5.0)+y.*3.822e+2-x.*y.*7.644e+2-6.37e+2./1.0e+1)+(y.*4.459e+1-y.^2.*1.3377e+2).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0),(x.*2.4e+1+y.*8.0-x.*y.*4.8e+1-4.0).*(x.*1.75812e+3+y.*1.22304e+3-x.*y.*3.97488e+3-4.9686e+2)-(x.*8.0-x.^2.*2.4e+1).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),(x.*2.4e+1+y.*8.0-x.*y.*4.8e+1-4.0).*(x.*3.822e+2+y.*2.548e+2-x.*y.*7.644e+2-6.37e+2./5.0)-(x.*8.0-x.^2.*2.4e+1).*(x.*(-8.918e+1)+x.^2.*1.3377e+2+1.11475e+1),(x.*2.4e+1+y.*8.0-x.*y.*4.8e+1-4.0).*(x.*7.644e+1+y.*1.1466e+2-x.*y.*2.2932e+2-3.822e+1)-(x.*8.0-x.^2.*2.4e+1).*(y.*(-8.918e+1)+y.^2.*1.3377e+2+1.11475e+1),-(x.*2.4e+1+y.*8.0-x.*y.*4.8e+1-4.0).*(x.*1.75812e+3+y.*7.644e+2-x.*y.*3.97488e+3-3.822e+2)+(x.*8.0-x.^2.*2.4e+1).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),(x.*2.4e+1+y.*8.0-x.*y.*4.8e+1-4.0).*(x.*3.822e+2+y.*(6.37e+2./5.0)-x.*y.*7.644e+2-6.37e+2./1.0e+1)+(x.*8.0-x.^2.*2.4e+1).*(x.*4.459e+1-x.^2.*1.3377e+2),(x.*8.0-x.^2.*2.4e+1).*(y.*(-8.918e+1)+y.^2.*1.3377e+2+1.11475e+1)-(x.*7.644e+1-x.*y.*2.2932e+2).*(x.*2.4e+1+y.*8.0-x.*y.*4.8e+1-4.0),(x.*2.2932e+2+y.*7.644e+2-x.*y.*3.97488e+3).*(x.*2.4e+1+y.*8.0-x.*y.*4.8e+1-4.0)-(x.*8.0-x.^2.*2.4e+1).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),-(x.*8.0-x.^2.*2.4e+1).*(x.*4.459e+1-x.^2.*1.3377e+2)-(y.*(6.37e+2./5.0)-x.*y.*7.644e+2).*(x.*2.4e+1+y.*8.0-x.*y.*4.8e+1-4.0),-(x.*8.0-x.^2.*2.4e+1).*(y.*4.459e+1-y.^2.*1.3377e+2)-(x.*3.822e+1-x.*y.*2.2932e+2).*(x.*2.4e+1+y.*8.0-x.*y.*4.8e+1-4.0),-(x.*2.4e+1+y.*8.0-x.*y.*4.8e+1-4.0).*(x.*2.2932e+2+y.*1.22304e+3-x.*y.*3.97488e+3-1.1466e+2)+(x.*8.0-x.^2.*2.4e+1).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),(x.*8.0-x.^2.*2.4e+1).*(x.*(-8.918e+1)+x.^2.*1.3377e+2+1.11475e+1)-(y.*2.548e+2-x.*y.*7.644e+2).*(x.*2.4e+1+y.*8.0-x.*y.*4.8e+1-4.0),(x.*2.4e+1+y.*8.0-x.*y.*4.8e+1-4.0).*(x.*3.822e+1+y.*1.1466e+2-x.*y.*2.2932e+2-1.911e+1)+(x.*8.0-x.^2.*2.4e+1).*(y.*4.459e+1-y.^2.*1.3377e+2),-(x.*1.6e+1-x.*y.*4.8e+1).*(x.*1.22304e+3+y.*1.75812e+3-x.*y.*3.97488e+3-4.9686e+2)-(y.*-1.6e+1+y.^2.*2.4e+1+2.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),-(x.*(-8.918e+1)+x.^2.*1.3377e+2+1.11475e+1).*(y.*-1.6e+1+y.^2.*2.4e+1+2.0)-(x.*1.6e+1-x.*y.*4.8e+1).*(x.*1.1466e+2+y.*7.644e+1-x.*y.*2.2932e+2-3.822e+1),-(y.*-1.6e+1+y.^2.*2.4e+1+2.0).*(y.*(-8.918e+1)+y.^2.*1.3377e+2+1.11475e+1)-(x.*1.6e+1-x.*y.*4.8e+1).*(x.*2.548e+2+y.*3.822e+2-x.*y.*7.644e+2-6.37e+2./5.0),(x.*1.6e+1-x.*y.*4.8e+1).*(x.*1.22304e+3+y.*2.2932e+2-x.*y.*3.97488e+3-1.1466e+2)+(y.*-1.6e+1+y.^2.*2.4e+1+2.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),(x.*4.459e+1-x.^2.*1.3377e+2).*(y.*-1.6e+1+y.^2.*2.4e+1+2.0)-(x.*1.6e+1-x.*y.*4.8e+1).*(x.*1.1466e+2+y.*3.822e+1-x.*y.*2.2932e+2-1.911e+1),(y.*-1.6e+1+y.^2.*2.4e+1+2.0).*(y.*(-8.918e+1)+y.^2.*1.3377e+2+1.11475e+1)+(x.*1.6e+1-x.*y.*4.8e+1).*(x.*2.548e+2-x.*y.*7.644e+2),-(x.*1.6e+1-x.*y.*4.8e+1).*(x.*7.644e+2+y.*2.2932e+2-x.*y.*3.97488e+3)-(y.*-1.6e+1+y.^2.*2.4e+1+2.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),-(x.*4.459e+1-x.^2.*1.3377e+2).*(y.*-1.6e+1+y.^2.*2.4e+1+2.0)+(x.*1.6e+1-x.*y.*4.8e+1).*(y.*3.822e+1-x.*y.*2.2932e+2),-(y.*4.459e+1-y.^2.*1.3377e+2).*(y.*-1.6e+1+y.^2.*2.4e+1+2.0)+(x.*1.6e+1-x.*y.*4.8e+1).*(x.*(6.37e+2./5.0)-x.*y.*7.644e+2),(x.*1.6e+1-x.*y.*4.8e+1).*(x.*7.644e+2+y.*1.75812e+3-x.*y.*3.97488e+3-3.822e+2)+(y.*-1.6e+1+y.^2.*2.4e+1+2.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),(x.*(-8.918e+1)+x.^2.*1.3377e+2+1.11475e+1).*(y.*-1.6e+1+y.^2.*2.4e+1+2.0)+(x.*1.6e+1-x.*y.*4.8e+1).*(y.*7.644e+1-x.*y.*2.2932e+2),(y.*4.459e+1-y.^2.*1.3377e+2).*(y.*-1.6e+1+y.^2.*2.4e+1+2.0)-(x.*1.6e+1-x.*y.*4.8e+1).*(x.*(6.37e+2./5.0)+y.*3.822e+2-x.*y.*7.644e+2-6.37e+2./1.0e+1),(x.*4.8e+1-x.*y.*1.92e+2).*(x.*1.22304e+3+y.*1.75812e+3-x.*y.*3.97488e+3-4.9686e+2)+(y.*4.8e+1-x.*y.*1.92e+2).*(x.*1.75812e+3+y.*1.22304e+3-x.*y.*3.97488e+3-4.9686e+2)+(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),(x.*4.8e+1-x.*y.*1.92e+2).*(x.*1.1466e+2+y.*7.644e+1-x.*y.*2.2932e+2-3.822e+1)+(y.*4.8e+1-x.*y.*1.92e+2).*(x.*3.822e+2+y.*2.548e+2-x.*y.*7.644e+2-6.37e+2./5.0)+(x.*(-8.918e+1)+x.^2.*1.3377e+2+1.11475e+1).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0),(x.*4.8e+1-x.*y.*1.92e+2).*(x.*2.548e+2+y.*3.822e+2-x.*y.*7.644e+2-6.37e+2./5.0)+(y.*4.8e+1-x.*y.*1.92e+2).*(x.*7.644e+1+y.*1.1466e+2-x.*y.*2.2932e+2-3.822e+1)+(y.*(-8.918e+1)+y.^2.*1.3377e+2+1.11475e+1).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0),-(x.*4.8e+1-x.*y.*1.92e+2).*(x.*1.22304e+3+y.*2.2932e+2-x.*y.*3.97488e+3-1.1466e+2)-(y.*4.8e+1-x.*y.*1.92e+2).*(x.*1.75812e+3+y.*7.644e+2-x.*y.*3.97488e+3-3.822e+2)-(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),(x.*4.8e+1-x.*y.*1.92e+2).*(x.*1.1466e+2+y.*3.822e+1-x.*y.*2.2932e+2-1.911e+1)+(y.*4.8e+1-x.*y.*1.92e+2).*(x.*3.822e+2+y.*(6.37e+2./5.0)-x.*y.*7.644e+2-6.37e+2./1.0e+1)-(x.*4.459e+1-x.^2.*1.3377e+2).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0),-(y.*(-8.918e+1)+y.^2.*1.3377e+2+1.11475e+1).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0)-(x.*4.8e+1-x.*y.*1.92e+2).*(x.*2.548e+2-x.*y.*7.644e+2)-(x.*7.644e+1-x.*y.*2.2932e+2).*(y.*4.8e+1-x.*y.*1.92e+2),(x.*4.8e+1-x.*y.*1.92e+2).*(x.*7.644e+2+y.*2.2932e+2-x.*y.*3.97488e+3)+(y.*4.8e+1-x.*y.*1.92e+2).*(x.*2.2932e+2+y.*7.644e+2-x.*y.*3.97488e+3)+(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),(x.*4.459e+1-x.^2.*1.3377e+2).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0)-(x.*4.8e+1-x.*y.*1.92e+2).*(y.*3.822e+1-x.*y.*2.2932e+2)-(y.*4.8e+1-x.*y.*1.92e+2).*(y.*(6.37e+2./5.0)-x.*y.*7.644e+2),(y.*4.459e+1-y.^2.*1.3377e+2).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0)-(x.*4.8e+1-x.*y.*1.92e+2).*(x.*(6.37e+2./5.0)-x.*y.*7.644e+2)-(x.*3.822e+1-x.*y.*2.2932e+2).*(y.*4.8e+1-x.*y.*1.92e+2),-(x.*4.8e+1-x.*y.*1.92e+2).*(x.*7.644e+2+y.*1.75812e+3-x.*y.*3.97488e+3-3.822e+2)-(y.*4.8e+1-x.*y.*1.92e+2).*(x.*2.2932e+2+y.*1.22304e+3-x.*y.*3.97488e+3-1.1466e+2)-(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),-(x.*(-8.918e+1)+x.^2.*1.3377e+2+1.11475e+1).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0)-(x.*4.8e+1-x.*y.*1.92e+2).*(y.*7.644e+1-x.*y.*2.2932e+2)-(y.*4.8e+1-x.*y.*1.92e+2).*(y.*2.548e+2-x.*y.*7.644e+2),(x.*4.8e+1-x.*y.*1.92e+2).*(x.*(6.37e+2./5.0)+y.*3.822e+2-x.*y.*7.644e+2-6.37e+2./1.0e+1)+(y.*4.8e+1-x.*y.*1.92e+2).*(x.*3.822e+1+y.*1.1466e+2-x.*y.*2.2932e+2-1.911e+1)-(y.*4.459e+1-y.^2.*1.3377e+2).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0),-(y.*8.0-x.*y.*4.8e+1).*(x.*1.75812e+3+y.*1.22304e+3-x.*y.*3.97488e+3-4.9686e+2)+(x.*8.0-x.^2.*2.4e+1).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),(x.*8.0-x.^2.*2.4e+1).*(x.*(-8.918e+1)+x.^2.*1.3377e+2+1.11475e+1)-(y.*8.0-x.*y.*4.8e+1).*(x.*3.822e+2+y.*2.548e+2-x.*y.*7.644e+2-6.37e+2./5.0),(x.*8.0-x.^2.*2.4e+1).*(y.*(-8.918e+1)+y.^2.*1.3377e+2+1.11475e+1)-(y.*8.0-x.*y.*4.8e+1).*(x.*7.644e+1+y.*1.1466e+2-x.*y.*2.2932e+2-3.822e+1),(y.*8.0-x.*y.*4.8e+1).*(x.*1.75812e+3+y.*7.644e+2-x.*y.*3.97488e+3-3.822e+2)-(x.*8.0-x.^2.*2.4e+1).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),-(x.*8.0-x.^2.*2.4e+1).*(x.*4.459e+1-x.^2.*1.3377e+2)-(y.*8.0-x.*y.*4.8e+1).*(x.*3.822e+2+y.*(6.37e+2./5.0)-x.*y.*7.644e+2-6.37e+2./1.0e+1),-(x.*8.0-x.^2.*2.4e+1).*(y.*(-8.918e+1)+y.^2.*1.3377e+2+1.11475e+1)+(x.*7.644e+1-x.*y.*2.2932e+2).*(y.*8.0-x.*y.*4.8e+1),-(y.*8.0-x.*y.*4.8e+1).*(x.*2.2932e+2+y.*7.644e+2-x.*y.*3.97488e+3)+(x.*8.0-x.^2.*2.4e+1).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),(x.*8.0-x.^2.*2.4e+1).*(x.*4.459e+1-x.^2.*1.3377e+2)+(y.*8.0-x.*y.*4.8e+1).*(y.*(6.37e+2./5.0)-x.*y.*7.644e+2),(x.*8.0-x.^2.*2.4e+1).*(y.*4.459e+1-y.^2.*1.3377e+2)+(x.*3.822e+1-x.*y.*2.2932e+2).*(y.*8.0-x.*y.*4.8e+1),(y.*8.0-x.*y.*4.8e+1).*(x.*2.2932e+2+y.*1.22304e+3-x.*y.*3.97488e+3-1.1466e+2)-(x.*8.0-x.^2.*2.4e+1).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),-(x.*8.0-x.^2.*2.4e+1).*(x.*(-8.918e+1)+x.^2.*1.3377e+2+1.11475e+1)+(y.*8.0-x.*y.*4.8e+1).*(y.*2.548e+2-x.*y.*7.644e+2),-(x.*8.0-x.^2.*2.4e+1).*(y.*4.459e+1-y.^2.*1.3377e+2)-(y.*8.0-x.*y.*4.8e+1).*(x.*3.822e+1+y.*1.1466e+2-x.*y.*2.2932e+2-1.911e+1),-(x.*8.0-x.*y.*4.8e+1).*(x.*1.22304e+3+y.*1.75812e+3-x.*y.*3.97488e+3-4.9686e+2)+(y.*8.0-y.^2.*2.4e+1).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),(y.*8.0-y.^2.*2.4e+1).*(x.*(-8.918e+1)+x.^2.*1.3377e+2+1.11475e+1)-(x.*8.0-x.*y.*4.8e+1).*(x.*1.1466e+2+y.*7.644e+1-x.*y.*2.2932e+2-3.822e+1),(y.*8.0-y.^2.*2.4e+1).*(y.*(-8.918e+1)+y.^2.*1.3377e+2+1.11475e+1)-(x.*8.0-x.*y.*4.8e+1).*(x.*2.548e+2+y.*3.822e+2-x.*y.*7.644e+2-6.37e+2./5.0),(x.*8.0-x.*y.*4.8e+1).*(x.*1.22304e+3+y.*2.2932e+2-x.*y.*3.97488e+3-1.1466e+2)-(y.*8.0-y.^2.*2.4e+1).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),-(x.*4.459e+1-x.^2.*1.3377e+2).*(y.*8.0-y.^2.*2.4e+1)-(x.*8.0-x.*y.*4.8e+1).*(x.*1.1466e+2+y.*3.822e+1-x.*y.*2.2932e+2-1.911e+1),-(y.*8.0-y.^2.*2.4e+1).*(y.*(-8.918e+1)+y.^2.*1.3377e+2+1.11475e+1)+(x.*8.0-x.*y.*4.8e+1).*(x.*2.548e+2-x.*y.*7.644e+2),-(x.*8.0-x.*y.*4.8e+1).*(x.*7.644e+2+y.*2.2932e+2-x.*y.*3.97488e+3)+(y.*8.0-y.^2.*2.4e+1).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),(x.*4.459e+1-x.^2.*1.3377e+2).*(y.*8.0-y.^2.*2.4e+1)+(x.*8.0-x.*y.*4.8e+1).*(y.*3.822e+1-x.*y.*2.2932e+2),(y.*8.0-y.^2.*2.4e+1).*(y.*4.459e+1-y.^2.*1.3377e+2)+(x.*8.0-x.*y.*4.8e+1).*(x.*(6.37e+2./5.0)-x.*y.*7.644e+2),(x.*8.0-x.*y.*4.8e+1).*(x.*7.644e+2+y.*1.75812e+3-x.*y.*3.97488e+3-3.822e+2)-(y.*8.0-y.^2.*2.4e+1).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),-(y.*8.0-y.^2.*2.4e+1).*(x.*(-8.918e+1)+x.^2.*1.3377e+2+1.11475e+1)+(x.*8.0-x.*y.*4.8e+1).*(y.*7.644e+1-x.*y.*2.2932e+2),-(y.*8.0-y.^2.*2.4e+1).*(y.*4.459e+1-y.^2.*1.3377e+2)-(x.*8.0-x.*y.*4.8e+1).*(x.*(6.37e+2./5.0)+y.*3.822e+2-x.*y.*7.644e+2-6.37e+2./1.0e+1),-(x.*4.8e+1+y.*9.6e+1-x.*y.*1.92e+2-2.4e+1).*(x.*1.22304e+3+y.*1.75812e+3-x.*y.*3.97488e+3-4.9686e+2)-(y.*4.8e+1-x.*y.*1.92e+2).*(x.*1.75812e+3+y.*1.22304e+3-x.*y.*3.97488e+3-4.9686e+2)-(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),-(x.*4.8e+1+y.*9.6e+1-x.*y.*1.92e+2-2.4e+1).*(x.*1.1466e+2+y.*7.644e+1-x.*y.*2.2932e+2-3.822e+1)-(y.*4.8e+1-x.*y.*1.92e+2).*(x.*3.822e+2+y.*2.548e+2-x.*y.*7.644e+2-6.37e+2./5.0)-(x.*(-8.918e+1)+x.^2.*1.3377e+2+1.11475e+1).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0),-(x.*4.8e+1+y.*9.6e+1-x.*y.*1.92e+2-2.4e+1).*(x.*2.548e+2+y.*3.822e+2-x.*y.*7.644e+2-6.37e+2./5.0)-(y.*4.8e+1-x.*y.*1.92e+2).*(x.*7.644e+1+y.*1.1466e+2-x.*y.*2.2932e+2-3.822e+1)-(y.*(-8.918e+1)+y.^2.*1.3377e+2+1.11475e+1).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0),(x.*4.8e+1+y.*9.6e+1-x.*y.*1.92e+2-2.4e+1).*(x.*1.22304e+3+y.*2.2932e+2-x.*y.*3.97488e+3-1.1466e+2)+(y.*4.8e+1-x.*y.*1.92e+2).*(x.*1.75812e+3+y.*7.644e+2-x.*y.*3.97488e+3-3.822e+2)+(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),-(x.*4.8e+1+y.*9.6e+1-x.*y.*1.92e+2-2.4e+1).*(x.*1.1466e+2+y.*3.822e+1-x.*y.*2.2932e+2-1.911e+1)-(y.*4.8e+1-x.*y.*1.92e+2).*(x.*3.822e+2+y.*(6.37e+2./5.0)-x.*y.*7.644e+2-6.37e+2./1.0e+1)+(x.*4.459e+1-x.^2.*1.3377e+2).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0),(x.*2.548e+2-x.*y.*7.644e+2).*(x.*4.8e+1+y.*9.6e+1-x.*y.*1.92e+2-2.4e+1)+(y.*(-8.918e+1)+y.^2.*1.3377e+2+1.11475e+1).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0)+(x.*7.644e+1-x.*y.*2.2932e+2).*(y.*4.8e+1-x.*y.*1.92e+2),-(x.*7.644e+2+y.*2.2932e+2-x.*y.*3.97488e+3).*(x.*4.8e+1+y.*9.6e+1-x.*y.*1.92e+2-2.4e+1)-(y.*4.8e+1-x.*y.*1.92e+2).*(x.*2.2932e+2+y.*7.644e+2-x.*y.*3.97488e+3)-(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),(y.*3.822e+1-x.*y.*2.2932e+2).*(x.*4.8e+1+y.*9.6e+1-x.*y.*1.92e+2-2.4e+1)-(x.*4.459e+1-x.^2.*1.3377e+2).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0)+(y.*4.8e+1-x.*y.*1.92e+2).*(y.*(6.37e+2./5.0)-x.*y.*7.644e+2),(x.*(6.37e+2./5.0)-x.*y.*7.644e+2).*(x.*4.8e+1+y.*9.6e+1-x.*y.*1.92e+2-2.4e+1)-(y.*4.459e+1-y.^2.*1.3377e+2).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0)+(x.*3.822e+1-x.*y.*2.2932e+2).*(y.*4.8e+1-x.*y.*1.92e+2),(x.*4.8e+1+y.*9.6e+1-x.*y.*1.92e+2-2.4e+1).*(x.*7.644e+2+y.*1.75812e+3-x.*y.*3.97488e+3-3.822e+2)+(y.*4.8e+1-x.*y.*1.92e+2).*(x.*2.2932e+2+y.*1.22304e+3-x.*y.*3.97488e+3-1.1466e+2)+(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),(y.*7.644e+1-x.*y.*2.2932e+2).*(x.*4.8e+1+y.*9.6e+1-x.*y.*1.92e+2-2.4e+1)+(x.*(-8.918e+1)+x.^2.*1.3377e+2+1.11475e+1).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0)+(y.*4.8e+1-x.*y.*1.92e+2).*(y.*2.548e+2-x.*y.*7.644e+2),-(x.*4.8e+1+y.*9.6e+1-x.*y.*1.92e+2-2.4e+1).*(x.*(6.37e+2./5.0)+y.*3.822e+2-x.*y.*7.644e+2-6.37e+2./1.0e+1)-(y.*4.8e+1-x.*y.*1.92e+2).*(x.*3.822e+1+y.*1.1466e+2-x.*y.*2.2932e+2-1.911e+1)+(y.*4.459e+1-y.^2.*1.3377e+2).*(x.*-4.8e+1-y.*4.8e+1+x.^2.*9.6e+1+y.^2.*9.6e+1+4.0),-(y.*1.6e+1-x.*y.*4.8e+1).*(x.*1.75812e+3+y.*1.22304e+3-x.*y.*3.97488e+3-4.9686e+2)-(x.*-1.6e+1+x.^2.*2.4e+1+2.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),-(x.*-1.6e+1+x.^2.*2.4e+1+2.0).*(x.*(-8.918e+1)+x.^2.*1.3377e+2+1.11475e+1)-(y.*1.6e+1-x.*y.*4.8e+1).*(x.*3.822e+2+y.*2.548e+2-x.*y.*7.644e+2-6.37e+2./5.0),-(x.*-1.6e+1+x.^2.*2.4e+1+2.0).*(y.*(-8.918e+1)+y.^2.*1.3377e+2+1.11475e+1)-(y.*1.6e+1-x.*y.*4.8e+1).*(x.*7.644e+1+y.*1.1466e+2-x.*y.*2.2932e+2-3.822e+1),(y.*1.6e+1-x.*y.*4.8e+1).*(x.*1.75812e+3+y.*7.644e+2-x.*y.*3.97488e+3-3.822e+2)+(x.*-1.6e+1+x.^2.*2.4e+1+2.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),(x.*4.459e+1-x.^2.*1.3377e+2).*(x.*-1.6e+1+x.^2.*2.4e+1+2.0)-(y.*1.6e+1-x.*y.*4.8e+1).*(x.*3.822e+2+y.*(6.37e+2./5.0)-x.*y.*7.644e+2-6.37e+2./1.0e+1),(x.*-1.6e+1+x.^2.*2.4e+1+2.0).*(y.*(-8.918e+1)+y.^2.*1.3377e+2+1.11475e+1)+(x.*7.644e+1-x.*y.*2.2932e+2).*(y.*1.6e+1-x.*y.*4.8e+1),-(y.*1.6e+1-x.*y.*4.8e+1).*(x.*2.2932e+2+y.*7.644e+2-x.*y.*3.97488e+3)-(x.*-1.6e+1+x.^2.*2.4e+1+2.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),-(x.*4.459e+1-x.^2.*1.3377e+2).*(x.*-1.6e+1+x.^2.*2.4e+1+2.0)+(y.*1.6e+1-x.*y.*4.8e+1).*(y.*(6.37e+2./5.0)-x.*y.*7.644e+2),-(y.*4.459e+1-y.^2.*1.3377e+2).*(x.*-1.6e+1+x.^2.*2.4e+1+2.0)+(x.*3.822e+1-x.*y.*2.2932e+2).*(y.*1.6e+1-x.*y.*4.8e+1),(y.*1.6e+1-x.*y.*4.8e+1).*(x.*2.2932e+2+y.*1.22304e+3-x.*y.*3.97488e+3-1.1466e+2)+(x.*-1.6e+1+x.^2.*2.4e+1+2.0).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),(x.*-1.6e+1+x.^2.*2.4e+1+2.0).*(x.*(-8.918e+1)+x.^2.*1.3377e+2+1.11475e+1)+(y.*1.6e+1-x.*y.*4.8e+1).*(y.*2.548e+2-x.*y.*7.644e+2),(y.*4.459e+1-y.^2.*1.3377e+2).*(x.*-1.6e+1+x.^2.*2.4e+1+2.0)-(y.*1.6e+1-x.*y.*4.8e+1).*(x.*3.822e+1+y.*1.1466e+2-x.*y.*2.2932e+2-1.911e+1),(x.*8.0+y.*2.4e+1-x.*y.*4.8e+1-4.0).*(x.*1.22304e+3+y.*1.75812e+3-x.*y.*3.97488e+3-4.9686e+2)-(y.*8.0-y.^2.*2.4e+1).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),(x.*8.0+y.*2.4e+1-x.*y.*4.8e+1-4.0).*(x.*1.1466e+2+y.*7.644e+1-x.*y.*2.2932e+2-3.822e+1)-(y.*8.0-y.^2.*2.4e+1).*(x.*(-8.918e+1)+x.^2.*1.3377e+2+1.11475e+1),(x.*8.0+y.*2.4e+1-x.*y.*4.8e+1-4.0).*(x.*2.548e+2+y.*3.822e+2-x.*y.*7.644e+2-6.37e+2./5.0)-(y.*8.0-y.^2.*2.4e+1).*(y.*(-8.918e+1)+y.^2.*1.3377e+2+1.11475e+1),-(x.*8.0+y.*2.4e+1-x.*y.*4.8e+1-4.0).*(x.*1.22304e+3+y.*2.2932e+2-x.*y.*3.97488e+3-1.1466e+2)+(y.*8.0-y.^2.*2.4e+1).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),(x.*8.0+y.*2.4e+1-x.*y.*4.8e+1-4.0).*(x.*1.1466e+2+y.*3.822e+1-x.*y.*2.2932e+2-1.911e+1)+(x.*4.459e+1-x.^2.*1.3377e+2).*(y.*8.0-y.^2.*2.4e+1),(y.*8.0-y.^2.*2.4e+1).*(y.*(-8.918e+1)+y.^2.*1.3377e+2+1.11475e+1)-(x.*2.548e+2-x.*y.*7.644e+2).*(x.*8.0+y.*2.4e+1-x.*y.*4.8e+1-4.0),(x.*7.644e+2+y.*2.2932e+2-x.*y.*3.97488e+3).*(x.*8.0+y.*2.4e+1-x.*y.*4.8e+1-4.0)-(y.*8.0-y.^2.*2.4e+1).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),-(x.*4.459e+1-x.^2.*1.3377e+2).*(y.*8.0-y.^2.*2.4e+1)-(y.*3.822e+1-x.*y.*2.2932e+2).*(x.*8.0+y.*2.4e+1-x.*y.*4.8e+1-4.0),-(y.*8.0-y.^2.*2.4e+1).*(y.*4.459e+1-y.^2.*1.3377e+2)-(x.*(6.37e+2./5.0)-x.*y.*7.644e+2).*(x.*8.0+y.*2.4e+1-x.*y.*4.8e+1-4.0),-(x.*8.0+y.*2.4e+1-x.*y.*4.8e+1-4.0).*(x.*7.644e+2+y.*1.75812e+3-x.*y.*3.97488e+3-3.822e+2)+(y.*8.0-y.^2.*2.4e+1).*(x.*(-2.6754e+2)-y.*2.6754e+2+x.^2.*5.3508e+2+y.^2.*5.3508e+2+2.2295e+1),(y.*8.0-y.^2.*2.4e+1).*(x.*(-8.918e+1)+x.^2.*1.3377e+2+1.11475e+1)-(y.*7.644e+1-x.*y.*2.2932e+2).*(x.*8.0+y.*2.4e+1-x.*y.*4.8e+1-4.0),(x.*8.0+y.*2.4e+1-x.*y.*4.8e+1-4.0).*(x.*(6.37e+2./5.0)+y.*3.822e+2-x.*y.*7.644e+2-6.37e+2./1.0e+1)+(y.*8.0-y.^2.*2.4e+1).*(y.*4.459e+1-y.^2.*1.3377e+2)],[12,12])
Error in integral2Calc>integral2t/tensor (line
228)
Z = FUN(X,Y); NFE = NFE + 1;
Error in integral2Calc>integral2t (line 55)
[Qsub,esub] = tensor(thetaL,thetaR,phiB,phiT);
Error in integral2Calc (line 9)
[q,errbnd] =
integral2t(fun,xmin,xmax,ymin,ymax,optionstruct);
Error in integral2 (line 106)
Q =
integral2Calc(fun,xmin,xmax,yminfun,ymaxfun,opstruct);
Error in ASS1 (line 28)
K_e = integral2(Ke,0,0.5,0,0.5)
Need a solution for this error. Thanks

  3 件のコメント

Ameer Hamza
Ameer Hamza 2020 年 10 月 20 日
Your function K_e returns a matrix for scalar inputs. Try
K_e(1, 1)
How do you want to integrate this matrix?
Walter Roberson
Walter Roberson 2020 年 10 月 20 日
integral2 cannot be used with matrix integration.
vinayak vijay
vinayak vijay 2020 年 10 月 20 日
I need to evaluate the above statement where the integral is a double integral over x and y. The limits for x an y are (0,0.5). Could you advise a way to do this?

サインインしてコメントする。

回答 (0 件)

製品


リリース

R2020a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by