typedef int[-2147483648,2147483647] int32_t;
const int N = 210;//210;//10-50-150-200
broadcast chan agads;
clock y, x, z;
bool PriE1a=false, NoPriE1a=true;
int32_t hrPTXEa = 0;
int flag1=1;
// ESTADOS ////////////////////////////////////////
int32_t Na3E1D = 0,
Na2E2D = 0,
K2E2D = 0,
E1D = N,
E2D = 0,
//Nuevos estados///////////////////////////////////////////
E1PTXEf = 0,
E1PTXEa = 0,
PTXEfa = 0,
PTXEff = 0, caPTXEast = 0,
// PROTÓMEROS /////////////////////////////////////////////////////////
Na2EPTXNa2E = 0,
K2E2PTXK2E = 0,
//New Ext. AM - 16-11-15
Na2E2PTXEa = 0,
PTXNa2ED = 0,
PTXK2ED = 0;
const double KdK = 12.0;
const double KdNa = 4.0;
double NaI = 150.0,
NaO = 160.0,
KI = 0.0,
KO = 0.0,
PTXo = 0.0;
const double alphaPTXEaf = 1.33356;
const double LalphaPTXEaf = 1.66e-4;
const double betaPTXEfa = 5.085;
const double alpha1 = 1.0e2;
const double beta1 = 1.0e-2;
const double alpha2 = 2.0e2;
const double beta2 = 8.0e2;
const double alpha3 = 1.0e-1;
const double beta3 = 3.0e2;
const double alpha4 = 50.0;
const double beta4 = 1.0e-1;
const double alphaP1 = 3.23e3;
const double betaP1New = 1.66e-4;
const double betaP1 = 1.66e-4;
const double alphaP2 = 4.3e2;
const double betaP2 = 2.1e-9;
const double alphaP3 = 5.15e2;
const double betaP3 = 6.6e-7;
const double alphaP4 = 1.93e-4;
const double betaP4 = 2.5e-4;
const double alphaP5= 1.7e-4;
const double betaP5 = 1.31e-4;
const double alphaP6 = 1.06e-1*7;
const double betaP6 = 2.0e-2;
const double alpha61New = 1.3353499339*10;
const double beta61New = 1.33356*0.1;
const double alphaP7 = 3.2e-2;
const double betaP7 = 1.07e-1;
const double alphaP8 = 8.8e-5;
const double betaP8 = 2.65e1;
const double alphaP9 = 2.15e-2;
const double betaP9 = 7.4e-3;
bool flag = false;
void Primeras(bool pe1a, bool nope1a)
{
if (not pe1a && E1PTXEa == 0 && not flag)
{
PriE1a = true;
NoPriE1a = false;
flag = true;
z=0;
}
else
if (PriE1a && not NoPriE1a)
{
PriE1a = false;
NoPriE1a = true;
}
}
ArtGadsFig2B
NewProt
NewProtOne
NewProtTwo
Prototipo3
CIArtGadsFig2B
LavadoPTXFig2B
ArtGadsFig5A
Str_Fig2
Str_Fig5
r1
double rbp, aux;
r_1
double rbp;
r2
double rbp, sr; //state rate: que representa un % de descuento en una parte de un total 20% de 10 es el 2% de 100, luego 20%*10/100.
r_2
double rbp, sr;
r3
double rbp;
r_3
double rbp;
r4
double rbp;
r_4
double rbp;
Template0
rp1
double rbp, sr;
r_p1
rp2
double rbp;
r_p2
rp3
double rbp;
r_p3
double rbp;
rp4
double rbp;
r_p4
double rbp;
rp6
double rbp, sr;
r_p6
double rbp, sr;
rp61
r_p61
rp8
r_p8
double rbp;
Template1
CIp1
CIp1One
CIp1KxNa
p1Popu
p_1
p2
p_2
p3
p_3
p4
p_4
p5
p_5
p51
p_51
p6New
p_6New
p60New
p_60New
p61New
p_61New
p62New
p_62New
Oldp6
Oldp_6
Oldp61
Oldp_61
p7
p_7
p71
p_71
p8
p_8
p81
p_81
p9
p_9
p91
p_91
p10
p_10
p101
p_101
p11
p_11
p12
p_12
p13
p_13
p14
p_14
p141
p_141
system
//Str_Fig5,//va con ptx=0
//ArtGadsFig5A,
//ArtGadsFig2B, //Colocar acá condiciones iniciales de parametros.. PTX=2
//CIArtGadsFig2B, //Saco para ver como va con la populación
// en las 40ut con las N bombas... luego reducir para ver si la probab aumenta
// ArtGadsFig2BtoNa3E1,//Para chequear cdo sale de a e f hacia E1 o Na3E1
//NewProt,
//NewProtOne,
//NewProtTwo,
Prototipo3,
r1, r2, r3, r4,
r_1, r_2, r_3, r_4,
rp1, r_p1,
rp2, r_p2, rp3, r_p3,
//rp4, r_p4, = p14!!
rp8, r_p8,
rp6, r_p6,
rp61, r_p61,
CIp1One,
//CIp1,//saca PTX cuando liga //E1PTXEf2E1PTXEa, Popu->Populación - CI -> Canais Isolados
p1Popu,// 1 o mas bombas... ArtGadsb2B saca la ptx a las 60 ut
//CIp1KxNa,
p_1,//E1PTXEa2E1PTXEf,
p2,//E1PTXEa2PTXEf_a,
p_2,//PTXEf_a2E1PTXEa,
p3,//PTXEf_a2PTXEf_f,
p_3,//PTXEf_f2PTXEf_a,
p4,//PTXEf_f2E1PTXEf,
p_4,//E1PTXEf2PTXEf_f,
p5,//E1PTXa2K2E2PTXK2E,//1 Saco AM
p_5,//K2E2PTXK2E2E1PTXa,
p51,
p_51,
//New estentios of AM - 16-11-15
p6New, p_6New,
p60New, p_60New,
p61New, p_61New,
p62New, p_62New,
/* Olds
p6, //Na2EPTXNa2E2E1PTXa,//2 Saco AM
p_6, //E1PTXEa2Na2EPTXNa2E,
p61,
p_61, */
p7,//PTXNa2E_2toPTXEf_a,//3 Saco AM
p_7,//PTXEf_a2PTXNa2E_2,
p71,
p_71,
p8,//PTXNa2E_2toPTXEf_f, //4
p_8,//PTXEf_f2PTXNa2E_2,
p81,
p_81,
p9,//PTXEf_a2PTXK2E_2, //5 Saco AM
p_9,//PTXK2E2PTXEf_a,
p91,
p_91,
p10,//PTXEf_f2PTXK2E_2, //6
p_10,//PTXK2EtoPTXEf_f,
p101,
p_101,
p11,//K2E2PTXK2EtoK2E2_2, //7
p_11,//K2E2_2toK2E2PTXK2E,
p12,//Na2EPTXNa2EtoNa2E2_2, //8
p_12,//Na2E2toNa2EPTXNa2E
p13,
p_13,
p14, //Na2E:PTXNa2E a E1:PTXEf ....... por p4, por NaO
p_14,
p141, //Na2E:PTXNa2E a E1:PTXEf ....... por p5, por NaI
p_141
;
Pr [<=800] (<>(z>50 && z<70 && PTXo== 0 && E1PTXEf ==0 && E1PTXEa ==0 && PTXEff==0 && PTXEfa==0
&& PTXNa2ED==0 && Na2EPTXNa2E==0
&& PTXK2ED==0 && K2E2PTXK2E==0))
Pregunta B - resetear el PTXo apenas liga, sino que vamos a tirar el sodio intra y extra y vamos a colocar
potasio intra y extra en igual cantidad que el sodio.
K2E2PTXK2E
Pr [<=300] (<>(z>50 && z<70 && PTXo== 0 && E1PTXEf ==0 && E1PTXEa ==0 && PTXEff==0 && PTXEfa==0 && PTXNa2ED==0 && Na2EPTXNa2E==0))
Pregunta A - Se resetea, cerar, apenas ligue el canal a abierto dejando Sodio intra y extra en los valores que traían.
Pr[<=800](<> (x<=50 && E1PTXEf>=1))
Pr [<=1200] (<>(z>150 && z<450 && PTXo== 0 && E1PTXEf ==0 && E1PTXEa ==0 && PTXEff==0 && PTXEfa==0))
Pr [<=600] (<>(z>150 && PTXo== 0 && E1PTXEf ==0 && E1PTXEa ==0 && PTXEff==0 && PTXEfa==0))
Pr [<=1200] (<>(z>1100 && PTXo== 0 && E1PTXEf ==0 && E1PTXEa ==0 && PTXEff==0 && PTXEfa==0))
Pr [<=1200] (<>(z>50 && z<450 && PTXo== 0 && E1PTXEf ==0 && E1PTXEa ==0 && PTXEff==0 && PTXEfa==0))
Pr [<=135] (<>(z>50 && z<70 && PTXo== 0 && E1PTXEf ==0 && E1PTXEa ==0 && PTXEff==0 && PTXEfa==0))
Pr [<=135] (<>(z>50 && z<70 && PTXo== 0 && E1PTXEf ==0 && E1PTXEa ==0 && PTXEff==0 && PTXEfa==0 && Na2EPTXNa2E>0))
En 210ut se tiene que colocar PTX entre [0, 60), al los 60 se saca la PTX y tiene que ir y venir de "f" a "a" hasta 90 ut,
esto da 150...
Pr[<=120] (A[](( E1PTXEf>0 && E1PTXEa==0) imply A<> E1PTXEa==1)
Pr [<=120 ] (<>(E1PTXEa >=1 && E1PTXEa <2))
Pr [<=120 ] (<>(E1PTXEf + E1PTXEa ==1))
Pr (<>[0,10] y>5&& y<=10 && (PTXEff + PTXEfa) ==N)
Pr [<=40 ] (<>(PTXEff + PTXEfa >=5))
Pr [<=60 ] (<>(y>15 && y <= 55 && PTXEff + PTXEfa >=N)) >= 0.15
Pregunta 2 –
¿Cuál es la probabilidad de que todas las bombas permanezcan entre PTXEff y PTXEfa durante 40 segundos?
Pr[<=120] ((E1PTXEf>0 && E1PTXEa==0 && y>=15) U (E1PTXEf + E1PTXEa==1 && y<=120))
Pr[<=120] ((E1PTXEf>0 && E1PTXEa==0 && y>=15) U (E1PTXEa==1 && PriE1a && y<=120))
Esta esta mal.... no tengo comosaber si f es igual a f-1!
Pr[<=120] ((E1PTXEf>0 && E1PTXEa==0 && y>=15) U (E1PTXEa==1 && E1PTXEf == E1PTXEf - 1 && y<=120))
Pr[<=120] ((E1PTXEf>0 && E1PTXEa==0 && y>=15) U (E1PTXEa==1 && y<=120))
Pregunta 1 - ¿Cuál es la probabilidad de que se ligue una bomba, esto es, de transitar de E1PTXEf para E1PTXEa,
cuando la primera se liga, en una populación, solo una bomba? No en muchas.
(E1PTXEf>0 && E1PTXEa==0) --> E1PTXEa==1
A[] (p imply A<> q)
Pr [<=1200 ] (<>(PTXNa2ED >=1))
Pr [ y<=950] (<>(x>=0 && x<=900 && (Na2E2PTXEa==1 || E1PTXEa ==1 || PTXEfa ==1)))
simulate 50 [ y<=200] {Na2EPTXNa2E, Na2E2PTXEa, E1D, E1PTXEf, E1PTXEa, PTXEfa, PTXEff, Na3E1D, Na2E2D, E2D, PTXNa2ED}
simulate 50 [ y<=1200] {Na2EPTXNa2E, Na2E2PTXEa, E1D, E1PTXEf, E1PTXEa, PTXEfa, PTXEff, Na3E1D, Na2E2D, E2D, PTXNa2ED}
simulate 5 [ y<=1200] {Na2EPTXNa2E, Na2E2PTXEa, E1D, E1PTXEf, E1PTXEa, PTXEfa, PTXEff, Na3E1D, Na2E2D, E2D, PTXNa2ED}
, Na3E1D, Na2E2D, E2D, PTXNa2ED
Pr [<=300] (<>(z<60 && Na2E2PTXEa==0) imply (z>=60 && z<=120 && Na2E2PTXEa==1))
Pr [<=300] (<>(z<60 && Na2E2PTXEa==0) && (z>=60 && z<=120 && Na2E2PTXEa==1))
simulate 1 [ y<=5000] {Na2EPTXNa2E, Na2E2PTXEa, E1D, E1PTXEf, E1PTXEa, PTXEfa, PTXEff}
simulate 1 [ y<=150] {Na2EPTXNa2E, Na2E2PTXEa, E1D, E1PTXEf, E1PTXEa, PTXEfa, PTXEff}
simulate 1 [ y<=120] {Na2E2PTXEa, E1D, E1PTXEf, E1PTXEa, PTXEfa, PTXEff}
simulate 1 [ y<=120] {Na2E2PTXEa, E1D, E1PTXEf, E1PTXEa, PTXEfa, PTXEff}
simulate 1 [ y<=95] {Na2E2PTXEa, E1D, E1PTXEf, E1PTXEa, PTXEfa, PTXEff}
simulate 1 [ y<=95] {E1D, E1PTXEf, E1PTXEa, PTXEfa, PTXEff}
simulate 10 [ y<=1000] {E1D, E1PTXEf, E1PTXEa, PTXEfa, PTXEff}
simulate 1 [ y<=1200] {Na3E1D, Na2E2D, K2E2D, E1D, E2D, E1PTXEf, E1PTXEa, PTXEfa, PTXEff, K2E2PTXK2E, Na2EPTXNa2E, PTXNa2ED, PTXK2ED}
Str_Fig5, PTXo=0
simulate 1 [ y<=850] {Na3E1D, Na2E2D, E1D, E2D, E1PTXEf, E1PTXEa, PTXEfa, PTXEff}
Pr [<=150 ] (<>(z>15 && PTXo== 0 && PTXEfa==0 && PTXEff==0))
Pr [<=150 ] (<>(z>15 && PTXo==0.0 && PTXEfa==0 && PTXEff==0))
simulate 1 [ y<=500] {Na3E1D, Na2E2D, K2E2D, E1D, E2D, E1PTXEf, E1PTXEa, PTXEfa, PTXEff, K2E2PTXK2E, Na2EPTXNa2E, PTXNa2ED, PTXK2ED}