Hello,
I'm trying to solve the triple integral in the attachment. I am so lost right now, because the constant A is actually a vector. Any help would be greatly appreciated. Thanks.
MATLAB: Solving a triple integral
triple integral
Related Solutions
Something like this?
psf = @(x,y,z)exp(-2*(x.^2+y.^2)-2*z.^2);fun1 = @(k,e,x,y,z)((e.*psf(x,y,z)).^k).*exp(-e.*psf(x,y,z))/factorial(k);% Make a function that takes a scalar k and a scalar e and returns the
% integral. Can use -100,100 limits (faster), or expressions involving k.
% Can use different tolerances.
scalar_k_scalar_e_fun = @(k,e)integral3(@(x,y,z)fun1(k,e,x,y,z),-inf,inf,-inf,inf,-inf,inf,'Abstol',1e-4,'RelTol',1e-3);% Make the latter function work with an array input for e, to facilitate
% plotting.
scalar_k_array_e_fun = @(k,e)arrayfun(@(e)scalar_k_scalar_e_fun(k,e),e);e = 0:0.1:1;k = 3:5;q = zeros(length(k),length(e));for i = 1:length(k) q(i,:) = scalar_k_array_e_fun(k(i),e);endhold onfor i = 1:length(k) plot(e,q(i,:));endhold off
Best Answer