Return to Answer

Use besselj --- the Bessel function of first kind --- to generate J1. I suppose you have to vary a and r to generate the "bubble".

I generated the following by varying x and y from -1:0.01:1 and plotting meshing points (x,y,f), I don't know if this is what you want.

###Code a = 2pi; [X Y] = meshgrid(-1:0.01:1,-1:0.01:1); R = sqrt(X.^2+Y.^2); f = (2besselj(1,a*R(:))./R(:)).^2; mesh(X,Y,reshape(f,size(X))); axis vis3d;

Log plot

Doresdoom suggestion, I replaced axis vis3d; with set(gca,'Zscale','Log').

alt text

###Mesh alt text

Use besselj --- the Bessel function of first kind --- to generate J1. I suppose you have to vary a and r to generate the "bubble".

I generated the following by varying x and y from -1:0.01:1 and plotting meshing points (x,y,f), I don't know if this is what you want.

Code

a = 2*pi;
[X Y] = meshgrid(-1:0.01:1,-1:0.01:1);
R = sqrt(X.^2+Y.^2);
f = (2*besselj(1,a*R(:))./R(:)).^2;
mesh(X,Y,reshape(f,size(X)));
axis vis3d;

Log plot

Doresdoom suggestion, I replaced axis vis3d; with set(gca,'Zscale','Log').

alt text

Mesh

alt text

Use besselj --- the Bessel function of first kind --- to generate J1. I suppose you have to vary a and r to generate the "bubble".

I generated the following by varying x and y from -1:0.01:1 and plotting meshing points (x,y,f), I don't know if this is what you want.

###Code a = 2pi; [X Y] = meshgrid(-1:0.01:1,-1:0.01:1); R = sqrt(X.^2+Y.^2); f = (2besselj(1,a*R(:))./R(:)).^2; mesh(X,Y,reshape(f,size(X))); axis vis3d;

Log plot

Doresdoom suggestion, I replaced axis vis3d; with set(gca,'Zscale','Log').

alt text http://lh6.ggpht.com/_4VSKzbK-Ieo/S-GTem8qv9I/AAAAAAAAAD8/PbgKKZuvoog/s1024/untitled.jpg

###Mesh alt text

Use besselj --- the Bessel function of first kind --- to generate J1. I suppose you have to vary a and r to generate the "bubble".

I generated the following by varying x and y from -1:0.01:1 and plotting meshing points (x,y,f), I don't know if this is what you want.

###Code a = 2pi; [X Y] = meshgrid(-1:0.01:1,-1:0.01:1); R = sqrt(X.^2+Y.^2); f = (2besselj(1,a*R(:))./R(:)).^2; mesh(X,Y,reshape(f,size(X))); axis vis3d;

Log plot

Doresdoom suggestion, I replaced axis vis3d; with set(gca,'Zscale','Log').

alt text

###Mesh alt text

Use besselj --- the Bessel function of first kind --- to generate J1. I suppose you have to vary a and r to generate the "bubble".

I generated the following by varying x and y from -1:0.01:1 and plotting meshing points (x,y,f), I don't know if this is what you want.

###Code [X Y] = meshgrid(-1:0.01:1,-1:0.01:1); R = sqrt(X.^2+Y.^2); f = (2besselj(1,aR(:))./R(:)).^2; mesh(X,Y,reshape(f,size(X))); axis vis3d;

###Mesh alt text

Use besselj --- the Bessel function of first kind --- to generate J1. I suppose you have to vary a and r to generate the "bubble".

I generated the following by varying x and y from -1:0.01:1 and plotting meshing points (x,y,f), I don't know if this is what you want.

###Code a = 2pi; [X Y] = meshgrid(-1:0.01:1,-1:0.01:1); R = sqrt(X.^2+Y.^2); f = (2besselj(1,a*R(:))./R(:)).^2; mesh(X,Y,reshape(f,size(X))); axis vis3d;

Log plot

Doresdoom suggestion, I replaced axis vis3d; with set(gca,'Zscale','Log').

alt text http://lh6.ggpht.com/_4VSKzbK-Ieo/S-GTem8qv9I/AAAAAAAAAD8/PbgKKZuvoog/s1024/untitled.jpg

###Mesh alt text