% Parameters
%------------------------------------------------------------
m = 1;   % mass
k = 1;   % boltzmann constant
T = 500; % temperature

lim = 100;
v   = linspace(-lim,lim);
p   = 1/sqrt(2*pi*m*k*T)*exp(-v.^2/(2*m*k*T));

plot(v,p);
