loes=NSolve[2 a^3/(a^2-1)==15,a][[2]]; a=a/.loes; r=2 a^3/(1+a^2)^2; q=2 a^3/(a^2-1); g[u_]=r(1-u/q); h[u_]=u/(1+u^2); abb=Plot[{g[u],h[u]},{u,0,q},ImageSize->{400,300}];