Modified Bessel function of the zero and first order

These are my codes

% Figure 1 % Motionless breather figure clear; clc; % hbar(Naturalpunits)=1; hbar=6.5821220*10e-16; Wp=10e12; T=300; K=8.617385*10e-5; W=0.7071; m=9.1093897*10e-31; beta=1/(K*T); d=(100:22:300)*10e-10; delta=0.01:0.05:0.5; za=(delta.*beta); %% B1=besseli(0,za); B2=besseli(1,za); [D,DELTA]=meshgrid(d,delta); [b2,b1]=meshgrid(B2,B1); %Wo=((Wp^2*n*d.^2.*delta.*B2)/(hbar^2*B1)); Q=Wp^2*D.^2.*DELTA.*b2.*m; R=(b1.*hbar^2); Wo=(Q./R).^2; mesh(Wo) %% %surfl(Wo) %colormap(jet) % change color m surf(D,DELTA,abs(Wo)); surf(d,delta,Wo); %surf(d,delta,abs(Wo.^2)) xlabel('d') ylabel('delta') zlabel('Wo^2') hold on grid on

[Torsten]

1. Please format your code (i.e. use the {} Code - button).

2. What is your question ?

Best wishes

Torsten.

[Samuel Suakye]

Not getting the expected graph As

[Samuel Suakye]

figure clear; clc; {hbar=6.5821220e-16}; {Wp=1.0e12;} {T=300;} {K=8.617385e-5;} {W=0.7071;} {m=9.1093897e-31;} {beta=1/(K*T);} {d=(100e-10:22e-10:300e-10);} {delta=0.01:0.05:0.5;} {za=(delta.*beta);} %% {I_0=besseli(0,za); } {I_1=besseli(1,za);}

{[I1,I0,] = meshgrid(I_1, I_0);} {[D,Delta] = meshgrid(d,delta);}

{Wo = {(((Wp.^2).*m.*(D.^2).*Delta.*I1)./((hbar.^2).*I0)).^0.5};}

%% {surf(d,delta,Wo)} %surf(D, Delta, Wo) {xlabel('d')} {ylabel('delta')} {zlabel('Wo^2')} {hold on} {grid on}

the problem is am suppose to get a contour like this but am not getting it