function praktika(shot,tlak)

	objemKomory = 0.079; % [m^3]
	objem = 0.057; % [m^3]
	k_B = 1.38e-23; % [J/K]
	
	k_Rogowski = 1.1e7;
	
	load([int2str(shot) '/CH1.txt']);
	load([int2str(shot) '/CH2.txt']);
	load([int2str(shot) '/CH3.txt']);
	load([int2str(shot) '/CH4.txt']);
	
	casy = CH1(:,1); % [s]
	dt = casy(2)-casy(1); % [s]
	
	derivaceProudu = k_Rogowski * CH1(:,2); % [A/s]
	proud = cumsum(derivaceProudu)*dt; % [s]
%	keyboard
	
	B_t = cumsum(CH4(:,2));
	H_alfa = CH2(:,2);
	
	H_alfa_c = conv(H_alfa,ones(101,1)/101,'same');
	H_alfa_sorted = sort(H_alfa_c(1:50:end));
	H_alfa_min = H_alfa_sorted(round(length(H_alfa_sorted)*0.1));
	H_alfa_max = H_alfa_sorted(round(length(H_alfa_sorted)*0.95));
	plazmaInd = H_alfa_c>H_alfa_min+(H_alfa_max-H_alfa_min)*0.1;
	delkaVyboje = sum(plazmaInd)*dt;
	fprintf('delka vyboje: %f ms\n',delkaVyboje*1000);
	
	U_loop = CH3(:,2);
	
	I_p = proud - U_loop/0.01029;
	I_p = I_p * 1.14;
	fprintf('maximalni I_p: %f kA\n',max(conv(I_p,ones(101,1)/101,'same'))/1000);
	
	P_ohm = U_loop .* I_p;
	fprintf('maximalni P_ohm: %f kW\n',max(conv(P_ohm,ones(101,1)/101,'same'))/1000);
	fprintf('prumerny P_ohm: %f kW\n',mean(P_ohm(plazmaInd))/1000);
	
	T_e = (0.7*(I_p./U_loop)).^0.6666666667;
	T_e(imag(T_e)~=0) = 0;
	T_e(T_e > 50) = 0;
	[T_e_max,T_e_max_ind] = max(conv(T_e,ones(101,1)/101,'same'));
	fprintf('maximalni T_e: %f eV\n',T_e_max);
	
	hustota = tlak/1000/k_B/300;
	E_th_max = 3*k_B*(T_e_max*11600)*objem*hustota;
	fprintf('maximalni E_th: %f J\n',E_th_max);
	
	I_p_conv = conv(I_p,ones(101,1)/101,'same');
	U_loop_conv = conv(U_loop,ones(101,1)/101,'same');
	fprintf('tau_E: %f ms\n',E_th_max/I_p_conv(T_e_max_ind)/U_loop_conv(T_e_max_ind)*1000);
	
	
	close all

	figure;
	plot(casy*1e3,U_loop,'g');
	xlabel('t [ms]');
	ylabel('U_loop [V]');

	figure;
	plot(casy*1e3,B_t,'b');
	xlabel('t [ms]');
	ylabel('B_t [??]');

	figure;
	plot(casy*1e3,I_p,'r');
	xlabel('t [ms]');
	ylabel('I_p [A]');
	
end