Thanks for your advice. I modified my integral as this:
syms x y z
ePp0=@(x,y,z) exp((y+2/3/(2*pi^2)*(1+2*pi^24*pi^2*z+2*pi^2*z.^2cos(2*pi*z)+2*pi*(1z).*sin(2*pi*z)))/2).*besseli(0,sqrt(2/3/(2*pi^2)*(1+2*pi^24*pi^2*z+2*pi^2*z.^2cos(2*pi*z)+2*pi*(1z).*sin(2*pi*z)).*y)).*exp((x+2/3/(2*pi^2)*(1+2*pi^2*z.^2cos(2*pi*z)2*pi*z.*sin(2*pi*z)))/2).*besseli(0,sqrt(2/3/(2*pi^2)*(1+2*pi^2*z.^2cos(2*pi*z)2*pi*z.*sin(2*pi*z)).*x));
Pp1=1/4*1/0.001*int(int(int(ePp0(x,y,z),z,0,0.001),y,0,x),x,0,100);
ans=VPA(Pp1,4);
But still I can not get the final answer. Could you please help me with this? Thanks a lot!
"Roger Stafford" wrote in message <ith0c7$gpg$1@newscl01ah.mathworks.com>...
> "Yue " <violinsonata5@gmail.com> wrote in message <itg580$4qq$1@newscl01ah.mathworks.com>...
> > Hi! Can someone help me on a triple integral?
> >
> > I have something like this:
> > integrand is f(x,y,z)=exp(yzcos(z))*besseli(0,sqrt((z+cos(z))*y)*exp(xzsin(z))*besseli(0,sqrt((z+sin(z))*x)
> >
> > int over z from 0 to b
> > int over y from 0 to Inf
> > int over x from 0 to y
> > I've tried using triplequad and quad2d and int but always get error messages.
> > Thanks a lot in advance!
>           
> 1. The parentheses don't balance in your expression for f(x,y,z). There are two more left parentheses than right ones.
>
> 2. The way f(x,y,z) is written doesn't support receiving x as a vector. The documentation states, "fun(x,y,z) must accept a vector x and scalars y and z, and return a vector of values of the integrand." There should therefore be a dot on the multiplication between the second 'exp' and the second 'besseli'.
>
> 3. The limits of integration should all be finite numbers.
>
> Roger Stafford
