Double integration of an array function

2 Ansichten (letzte 30 Tage)
Athira T Das
Athira T Das am 17 Aug. 2022
Kommentiert: Athira T Das am 17 Aug. 2022
clc; clear all; close all;
lambda=532*1e-9;
w=0.35*1e-3;
R=3500*1e-3;
alpha=40;
rho=789;
cp=2440;
beta=750*1e-6;
mew=1.36*1e-6;
nt=4*1e-4;
K=0.17;
d=2100*1e-3;
h=1.3*1e-3;
P=40*1e-3;
L=10*1e-3;
D=K/(rho*cp);
tc=(w^2)/(4*D);
theta=(nt*alpha*P*L)/(lambda*K);
U=sqrt((2*P)/(pi*w*w));
J=(1i*pi*w*w)/(lambda*d);
k=2*pi/lambda;
deltaT=1;
g=9.8;
ti=1;
syms t
xx =linspace(-30,30);
yy =linspace(-30,30);
v_x = (beta*g*deltaT*pi*(h^2))/(16*mew);
F1 = (((exp(-2.*((((xx-(v_x.*t)).^2) + (yy.^2))./((1+((2.*t)./(tc))).*(w.^2)))))-(exp(-2.*(((v_x.*t).*2)./((1+((2.*t)./tc)).*(w.^2))))))./(1+((2.*t)./tc)));
fun1 = matlabFunction(F1,'Vars',[t]);
I1=integral(fun1,0,ti,'ArrayValued',true);
del_phi=I1.*(theta./(2.*tc))
del_phi = 1×100
-0.0017 -0.0017 -0.0017 -0.0017 -0.0017 -0.0017 -0.0017 -0.0017 -0.0017 -0.0017 -0.0017 -0.0017 -0.0017 -0.0017 -0.0017 -0.0017 -0.0017 -0.0017 -0.0017 -0.0017 -0.0017 -0.0017 -0.0017 -0.0017 -0.0017 -0.0017 -0.0017 -0.0017 -0.0017 -0.0017
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
syms x y
F2 = exp(-((x.^2)+(y.^2))./(w.^2)).*exp(1i.*(((((x.^2)+(y.^2))./(2.*R)).*(-k)) + del_phi)).*exp(-1i.*k.*(x.*xx+y.*yy)./d);
fun2 = matlabFunction(F2,'Vars',[x,y])
fun2 = function_handle with value:
@(x,y)[exp(x.*1.687214099672284e+8i+y.*1.687214099672284e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*1.653128966345571e+8i+y.*1.653128966345571e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*1.619043833018858e+8i+y.*1.619043833018858e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*1.584958699692146e+8i+y.*1.584958699692146e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*1.550873566365433e+8i+y.*1.550873566365433e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*1.51678843303872e+8i+y.*1.51678843303872e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*1.482703299712007e+8i+y.*1.482703299712007e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*1.448618166385294e+8i+y.*1.448618166385294e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*1.414533033058582e+8i+y.*1.414533033058582e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*1.380447899731869e+8i+y.*1.380447899731869e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*1.346362766405156e+8i+y.*1.346362766405156e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*1.312277633078443e+8i+y.*1.312277633078443e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*1.27819249975173e+8i+y.*1.27819249975173e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*1.244107366425017e+8i+y.*1.244107366425017e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6).*exp(x.*1.210022233098305e+8i+y.*1.210022233098305e+8i),exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6).*exp(x.*1.175937099771592e+8i+y.*1.175937099771592e+8i),exp(x.*1.141851966444879e+8i+y.*1.141851966444879e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*1.107766833118166e+8i+y.*1.107766833118166e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*1.073681699791453e+8i+y.*1.073681699791453e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*1.039596566464741e+8i+y.*1.039596566464741e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*1.005511433138028e+8i+y.*1.005511433138028e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*9.71426299811315e+7i+y.*9.71426299811315e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*9.373411664846022e+7i+y.*9.373411664846022e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*9.032560331578894e+7i+y.*9.032560331578894e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*8.691708998311766e+7i+y.*8.691708998311766e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*8.350857665044638e+7i+y.*8.350857665044638e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*8.01000633177751e+7i+y.*8.01000633177751e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*7.669154998510382e+7i+y.*7.669154998510382e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*7.328303665243254e+7i+y.*7.328303665243254e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*6.987452331976126e+7i+y.*6.987452331976126e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*6.646600998708997e+7i+y.*6.646600998708997e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*6.305749665441869e+7i+y.*6.305749665441869e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*5.964898332174741e+7i+y.*5.964898332174741e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*5.624046998907613e+7i+y.*5.624046998907613e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*5.283195665640485e+7i+y.*5.283195665640485e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*4.942344332373357e+7i+y.*4.942344332373357e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6).*exp(x.*4.601492999106229e+7i+y.*4.601492999106229e+7i),exp(x.*4.260641665839101e+7i+y.*4.260641665839101e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6).*exp(x.*3.919790332571973e+7i+y.*3.919790332571973e+7i),exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6).*exp(x.*3.578938999304845e+7i+y.*3.578938999304845e+7i),exp(x.*3.238087666037717e+7i+y.*3.238087666037717e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*2.897236332770589e+7i+y.*2.897236332770589e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*2.556384999503461e+7i+y.*2.556384999503461e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*2.215533666236332e+7i+y.*2.215533666236332e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*1.874682332969204e+7i+y.*1.874682332969204e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*1.533830999702076e+7i+y.*1.533830999702076e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6).*exp(x.*1.192979666434948e+7i+y.*1.192979666434948e+7i),exp(x.*8.521283331678202e+6i+y.*8.521283331678202e+6i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*5.112769999006921e+6i+y.*5.112769999006921e+6i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*1.70425666633564e+6i+y.*1.70425666633564e+6i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-1.70425666633564e+6i-y.*1.70425666633564e+6i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-5.112769999006921e+6i-y.*5.112769999006921e+6i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-8.521283331678202e+6i-y.*8.521283331678202e+6i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6).*exp(x.*-1.192979666434948e+7i-y.*1.192979666434948e+7i),exp(x.*-1.533830999702076e+7i-y.*1.533830999702076e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-1.874682332969204e+7i-y.*1.874682332969204e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-2.215533666236332e+7i-y.*2.215533666236332e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-2.556384999503461e+7i-y.*2.556384999503461e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-2.897236332770589e+7i-y.*2.897236332770589e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-3.238087666037717e+7i-y.*3.238087666037717e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6).*exp(x.*-3.578938999304845e+7i-y.*3.578938999304845e+7i),exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6).*exp(x.*-3.919790332571973e+7i-y.*3.919790332571973e+7i),exp(x.*-4.260641665839101e+7i-y.*4.260641665839101e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6).*exp(x.*-4.601492999106229e+7i-y.*4.601492999106229e+7i),exp(x.*-4.942344332373357e+7i-y.*4.942344332373357e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-5.283195665640485e+7i-y.*5.283195665640485e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-5.624046998907613e+7i-y.*5.624046998907613e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-5.964898332174741e+7i-y.*5.964898332174741e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-6.305749665441869e+7i-y.*6.305749665441869e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-6.646600998708997e+7i-y.*6.646600998708997e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-6.987452331976126e+7i-y.*6.987452331976126e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-7.328303665243254e+7i-y.*7.328303665243254e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-7.669154998510382e+7i-y.*7.669154998510382e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-8.01000633177751e+7i-y.*8.01000633177751e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-8.350857665044638e+7i-y.*8.350857665044638e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-8.691708998311766e+7i-y.*8.691708998311766e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-9.032560331578894e+7i-y.*9.032560331578894e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-9.373411664846022e+7i-y.*9.373411664846022e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-9.71426299811315e+7i-y.*9.71426299811315e+7i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-1.005511433138028e+8i-y.*1.005511433138028e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-1.039596566464741e+8i-y.*1.039596566464741e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-1.073681699791453e+8i-y.*1.073681699791453e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-1.107766833118166e+8i-y.*1.107766833118166e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-1.141851966444879e+8i-y.*1.141851966444879e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6).*exp(x.*-1.175937099771592e+8i-y.*1.175937099771592e+8i),exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6).*exp(x.*-1.210022233098305e+8i-y.*1.210022233098305e+8i),exp(x.*-1.244107366425017e+8i-y.*1.244107366425017e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-1.27819249975173e+8i-y.*1.27819249975173e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-1.312277633078443e+8i-y.*1.312277633078443e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-1.346362766405156e+8i-y.*1.346362766405156e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-1.380447899731869e+8i-y.*1.380447899731869e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-1.414533033058582e+8i-y.*1.414533033058582e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-1.448618166385294e+8i-y.*1.448618166385294e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-1.482703299712007e+8i-y.*1.482703299712007e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-1.51678843303872e+8i-y.*1.51678843303872e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-1.550873566365433e+8i-y.*1.550873566365433e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-1.584958699692146e+8i-y.*1.584958699692146e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-1.619043833018858e+8i-y.*1.619043833018858e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-1.653128966345571e+8i-y.*1.653128966345571e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6),exp(x.*-1.687214099672284e+8i-y.*1.687214099672284e+8i).*exp(x.^2.*-1.687214099672284e+6i-y.^2.*1.687214099672284e+6i-1.742405633866455e-3i).*exp(x.^2.*(-8.16326530612245e+6)-y.^2.*8.16326530612245e+6)]
q = integral2(fun2,-Inf,Inf,-Inf,Inf,'ArrayValued',true)
Error using integral2ParseArgs
Invalid argument name 'ArrayValued'. Name must be 'AbsTol', 'RelTol', or 'Method'.

Error in integral2 (line 103)
opstruct = integral2ParseArgs(isImproper,varargin{:});

Antworten (1)

Torsten
Torsten am 17 Aug. 2022
"integral2" has no "ArrayValued" option - only "integral" has.
But since the function comes from a "MatlabFunction" command, I guess you don't need this option in your call to "integral2".
  3 Kommentare
Torsten
Torsten am 17 Aug. 2022
For a given (x,y) pair, fun2 must return a scalar value z = fun2(x,y).
In your case, it returns a vector of size 1x100.
Please state your problem in a mathematical description - for sure, xx and yy cannot be arrays.
Athira T Das
Athira T Das am 17 Aug. 2022
That's right. Let me recheck the function.

Melden Sie sich an, um zu kommentieren.

Kategorien

Mehr zu MATLAB finden Sie in Help Center und File Exchange

Community Treasure Hunt

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

Start Hunting!

Translated by