/* evaluates long-linear approximation */
a1 = (1-_bigsb)/_gam;
b1 = _bigsb*(_gam+1)/_gam - 1;
c1 = _gam/(1-_phi)-_bigsb;
disc = b1^2 - 4*a1*c1;
if disc < 0;
"discriminant less than zero, can't compute slope";
"a" a1;
"b" b1;
"c" c1;
"disc" disc;
disc = 0;
else;
disc = (disc)^0.5;
endif;
b = seqb(-10,10,0.01);
newb = ( _gam/(1-_phi) - (1-_bigsb)*b - _bigsb).*(1-_phi)./_gam.*(b-_gam);
graphjc;
xy(b,b~newb);
/*
graphjc;
px = _sg;
_plegstr = "p/d\000p/dx\000P/da\000ppa";
_plegctl = 1;
py =_lnpd~_lnpdx~_lnpda~_lnppa;
xy(px,py);
*/