loes=NSolve[2 b^3/(b^2-1)==15,b]; a=b/.loes[[1]]; 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}];