function Q_new = Tq_( X, D, Alpha , h )
%Tq is a discrete stocastic kernel used to calculate the discrete state
%evolution in DTSHS model.
%   Q_new: new discrete states of heater; On => 1, OFF=>0.
%   X: Continous state at time k (temp), 
%   D, Alpha: sigmoidal fun parameters.


Q_new = zeros(1, 5);

%   Calucate the Sigmoidal  fun Rho
%   Rho = Sigmodal_DTSHS( X, D, Alpha);

    Rho = (X .^ D) ./ ((Alpha .^ D) + (X .^ D));




%Sample new states 

Q_new = random('Binomial',1,(1-Rho));


end



