how to plot double integral of bessel function?

I want to plot the below written function but it showing an error "' int' for input arguments of type 'double'. "

k = 10^4;
z = 0:10^(-6):10^(-5);
y = z.*besselj(0,(k*z));
a = int(y,z);
b = a/z;
c = b^0.5*(-266.46);
d = int(c,z);
plot(z,d);

Answer 1

[Mikhail]

Do you have Symbolic Math Toolbox? Function int is in this toolbox.