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)
E(10) Here the error!!
Best Answer