#!/usr/bin/gnuplot # maxwell_rozdeleni.plt # Maxwellovo rozdeleni rychlosti. # Miroslav Broz (miroslav.broz@email.cz), May 5th 2011 k = 1.38e-23 # J/K; Boltzmannova konstanta m_u = 1.66e-27 # kg T = 300. # K m_air = 29.*m_u # kg; "molekula vzduchu" m_hydrogen = 1.*m_u # kg; atom vodiku f(v,m) = 4./sqrt(pi) * (m/(2.*k*T))**(3./2.) * v**2 * exp(-m*v**2/(2.*k*T)) ######################################################################## set xl "v / m/s" set yl "f(v) dv" set xr [0:10000] #set logscale y set samples 1000 p f(x,m_air),\ f(x,m_hydrogen) pa -1 set term png small set out "maxwell_rozdeleni.png" rep