%

*A program to plot the dispersion relation for electromagnetic waves in*

%

*plasma with B_o = 0*

clear all

close all

%

*assign the value of constants*

no = 1e18; %

*neutral density in m^-3*

e = 1.6021e-19; %

*charge of electron in C*

Eo = 8.8541e-12; %

*permittivity of free space in Fm^-1*

me = 9.1095e-31; %

*mass of electron in kg*

w_p = sqrt(no*e^2/Eo/me);%

*plasma frequency*

c = 3e8; %

*speed of light in vacuum*

k = linspace(-2*pi/2000,2*pi/2000,1000);

%

*mathematical relation of w with k*

w = sqrt(w_p^2 + c^2*k.^2); %

*frequency of electromagnetic wave*

plot(k,w)

xlabel('wave vector(k)'), ylabel('ferquency of electromagnetic wave (w)'),title('Dispersion relation for electromagnetic waves in plasma with B_o = 0')

grid on

## No comments:

Post a Comment