function AlphaP=FormAlphaP(Init_L,deltL,Init_U,deltU)
tdeltLinx=find(deltL<0);
tdeltUinx=find(deltU<0);
t1=-Init_L(tdeltLinx)./deltL(tdeltLinx);
t2=-Init_U(tdeltUinx)./deltU(tdeltUinx);
t3=[t1;t2];
t4=min(t3);
AlphaP=0.9995*min([t4 1]);
end