* Reset for initial conditions

aa1 = a1;
aa2 = a2;
aa3 = a3;

PBACK =  1.17 ;
PARTFRACT1 = 1;
PARTFRACT2 = 1;
PARTFRACT21 = 1;
partfract(t) = partfract21;
cost1(T) = (PBACK*SIGMA(T)/EXPCOST2)* ( (BACKRAT-1+ EXP (-gback* (ORD(T)-1) ) )/BACKRAT);
PARTFRACT(T)$(ord(T)<25) = Partfract21 + (PARTFRACT2-Partfract21)*exp(-DPARTFRACT*(ORD(T)-2));
partfract("1")= PARTFRACT1;

TATM.up(t)  = 10 ;
mat.up(T)= 4000;
miu.up(t)= 1;
miu.lo(t)= .001;
k.lo(t) = 1;
k.up(t) = 1000000;
K0 = 137;
miu.fx("1")=.005;
partfract("1")=	0.25372	;