No error when I simulate the code.
What are you trying to do?
r=50; d=50; H0=50; l=5; b=5; R=1.33;
a=(2*l*H0/pi)^0.5;
Omega=1;
va=(3*(pi)^(3/2))/(16*sqrt(2))*(Omega*sqrt(l*H0))/(1+pi/4*R*(l/b));
J=@(x,k)(besselj(0,k.*a.*x)).*sqrt(1-x.^2).*x;
B=@(k)integral(@(x)J(x,k),0,1)
J0=@(k)besselj(0,k.*r);
F=@(t)Omega.*sin(Omega.*t);
B1=@(k,t)J0(k).*B(k).*F(t).*k
E=@(t)integral(@(k)B1(k,t),0,10)
Command Window
B =
@(k)integral(@(x)J(x,k),0,1)
B1 =
@(k,t)J0(k).*B(k).*F(t).*k
E =
@(t)integral(@(k)B1(k,t),0,10)
>>
Best Answer