function [ c, ceq ] = confuneq_sus( x )
%UNTITLED3 Summary of this function goes here

%% konstanty

k = 1.0015;                 %koeficient Pitotovi trubice
D = 0.250;                  %průměr spalinovodo     [m]              (opravit/zkontrolovat)
S = (pi*D^2)/4;             %průřez (obsah) spalinovodu     [m2]

%%


%   Nonlinear inequality constraints
c = [];                %nejsou
%   Nonlinear equality constraints

% Výpočet rychlosti
rychlost = k*((574.2*x(2)+156842.77)./(x(8)*1000)).^(1/2).*(x(3)).^(1/2);

% Výpočet průtoku
prutok     = S*rychlost;  %prutok "spalin" [m3/s]
prutok_in = prutok.*(x(7)./x(2));
prutok_hod = prutok*3600; %prutok "spalin" [m3/h] - pro kontrolu, měl by být okolo 600m3/h

% Výpočet saturačního tlaku
P_sat_in = (1./((x(7)+273.15).^8.2)).*exp(77.345+0.0057.*(x(7)+273.15)-7235.*(x(7)+273.15).^(-1));       %[Pa]
P_sat_out = (1./((x(2)+273.15).^8.2)).*exp(77.345+0.0057.*(x(2)+273.15)-7235.*(x(2)+273.15).^(-1));      %[Pa]

% Výpočet absolutní vlhkosti
abs_vlhkost_in  = 0.62198.*(((x(6)./100).*P_sat_in)./((x(8).*1000)-((x(6)./100).*P_sat_in)));       %[kg_vody/kg_such.vzduchu]
abs_vlhkost_out = 0.62198.*(((x(1)./100).*P_sat_out)./((x(8).*1000)-((x(1)./100).*P_sat_out)));     %[kg_vody/kg_such.vzduchu]


% Výpočet hustoty
ro_in = (1.316e-3./(x(7)+273.15)).*(2.65.*(x(8).*1000)-(x(6)./100).*P_sat_in);   %[kg/m3] 
%ro_in_1 = (3.483e-3./(teplota_lab+273.15)).*(tlak_lab.*1000-0.379.*(rel_vlhkost_lab./100).*P_sat_in);   %[kg/m3]     
ro_out = (1.316e-3./(x(2)+273.15)).*(2.65.*(x(8).*1000)-(x(1)/100).*P_sat_out);   %[kg/m3]

% Výpočet průtoku v [kg/s]
prutok_kg_in = ro_in.*prutok_in;                       %[kg/s] - nesprávný
prutok_kg_out = ro_out.*prutok;                  %[kg/s]

% Výpočet odpařené vody
MW_stredni_such.vzduch = 28.964;              %[kg/kmol]
V_m = 22.414;                                 %[m3N/kmol]

m_voda_in  = (abs_vlhkost_in.*prutok.*MW_stredni_such.vzduch)./V_m;  %[kg/s]    (přepočet na normální m3 ??)
m_voda_out = (abs_vlhkost_out.*prutok.*MW_stredni_such.vzduch)./V_m; %[kg/s]
m_delta    = m_voda_out-m_voda_in;                                   %[kg/s]
%m_evap = sum(m_delta);                                                          %(opravit) nesmyslné číslo

% Výpočet entalpie
c_a = 1010;       %[J/kg.K]
c_v = 1840;       %[J/kg.K]
l   = 2500;       %[kJ/kg]
h_in = c_a.*x(7)+abs_vlhkost_in.*(l*1000+c_v.*x(7));      %[J/kg]   %(dosadit teplotu v Kelvinech? zkontrolovat)
h_out = c_a.*x(2)+abs_vlhkost_out.*(l*1000+c_v.*x(2));      %[J/kg]   %(dosadit teplotu v Kelvinech? zkontrolovat)

% Ohřev vzduchu
Q_h_vz = prutok_kg_out.*(h_out-h_in)/1000;       %[kW]

% Teplo potřebné na odpaření
Q_evap = m_delta*l*1000/1000;                    %[kW]

% Tepelné ztráty se neuvažují


%0 = Q_h - Q_o + Q_i - Q_evap
%ceq(1) = x(5) - Q_h_vz - Q_evap;
ceq(1) = x(5) - Q_h_vz - Q_evap;

%   Hmotnostní bilance
%ceq(2) = m_voda_in + hmotnost(i-1)-hmotnost(length(hmotnost)) - m_voda_out - x(4) -hmotnost(length(hmotnost));
ceq(2) = m_voda_in + x(4) - m_voda_out;
end

