%Vojt�ch Laitl 2016
%S.V.A.T. FJFI �VUT
epsilon = 8.8542*10^-12
e = 1.602*10^-19
k_B = 1.38*10^-23
m_el = 9.109*10^-31
c = 3*10^8
h = 6.626*10^-34
load loop_voltage_21072.txt
U = loop_voltage_21072(:,2)
R = 0.0097 %read from the data for discharge No.21072
P0 = 0.001
P = U.^2./R
L = 10.*log10(P./P0) %gain as a ratio of power disipated is solved
save -ascii gain.txt L
dn = -1/4.*L.^-4 %basic differential gear sovled
E_k = - (36*pi)^-2.*epsilon^-4/3.*n/e^-4
p = [0; [diff(E_k)]]
omega = p.*c./h %electron plasma frequency during peak events determined
save -ascii electron_plasma_frequency.txt omega
n = (omega.^2*m_el*epsilon)/e^2 %peak electron density approached
save -ascii peak_electron_density.txt n
n0 = 3.7*10^18 %approximation of electron density through equilibrium state
N = n+n0 %converges to the mean value
save -ascii total_electron_density.txt N
dT = -E_k./k_B
T0 = 21.9*11600
lambda_De = sqrt((epsilon.*k_B.*T0)./(10^17.*e^2)) %debye radius approach done within the basic model
T = dT.*(N./lambda_De)./11600 %model approximation
save -ascii temperature_eV.txt T
T_sum = T + T0 %solved within the equilibrium state determination
save -ascii total_temperature.txr T_sum
dN_D = 4/3*pi.*lambda_De.*n %plasma parameter growth in differential gear
save -ascii plasma_paramter_growth.txt dN_D
N_D = 4/3*pi.*lambda_De.*N %total plasma parameter
save -ascii total_plasma_parameter.txt N_D
W_k = 3/2.*k_B.*T %kinetic energy growth
save -ascii kinetic_energy_peak.txt W_k
W_k_sum = 3/2.*k_B.*T_smu %total kinetic energy
save -ascii total_kinetic_energy.txt W_k_sum
W_p = (2.*W_k)./(9.*N_D.^2/3) %potential energy peak
save -ascii potential_energy_peak.txt W_p
W_p_sum = (2.*W_k_sum)./(9.*N_D.^2/3) %total potential energy
save -ascii total_potential_energy.txt W_p_sum
v = (3*k_B.*T./(n.*m_el)).^1/2 .* lambda_De %velocity variations
save -ascii velocity.txt v %total right the basic state has the zero derivation of fluctuation scope
a = [0;[diff(v)]]
save -ascii acceleration.txt a
s = (2*k_B.*T.^3./(n.*m_el)).^1/2 .*lambda_De.^3 %total right if the basic-state 
save -ascii fluctuation_scope.txt s