function yhat = and_1028(beta, w)
%Parameter Assignations
A = beta(1);
B = beta(2);
C = beta(3);
D = beta(4);
E = beta(5);

%independent variables
arab = w(:,1);
theo = w(:,2);


%Model for Pbad, (fn(arab, aTc))
C0_araC = 10^1 ;%uM
%n_pBad = 2.8;

 
%Kd_pBad = .09*10^3;%(uM)(?) source: openwetware supplemental materials
%Kd_tbs = 1.3*10^-3;%(uM) for tal8, for tal 1 the value is 240 nM, source: Meckler et al.
%Kd_RS =  .32;%source: Gallivan et al (uM)


%C_araC = (C0_araC) .* (arab.^n_pBad ./ (Kd_pBad.^n_pBad + arab.^n_pBad));

%Cf_araC = C0_araC - C_araC;

%K1_pBad = 0.009;  K2_pBad = 37.5./C0_araC;  K3_pBad = 3.4./C0_araC; 

%PBad = (K1_pBad + Kd_pBad.*C_araC)./ (1+K1_pBad+K2_pBad.*C_araC+ K3_pBad .*Cf_araC);

P = 1;%relative promoter strength

%Model 
yhat = A*P./(B*(1+ (theo.*(C./(1+D./arab))./E)));
