phi=seq(-.95,.95,0.05)
factor=(1+phi)/(1-phi)
plot(y=factor,x=phi,type='l')
abline(h=1,col='blue')