function fv = c_4_p_7_fn(x)
global r;
denom = ( 1 + x.^2 ) .* ( 1 + r^2 * x.^2 ) .* ( 1 + r^4 * x.^2 );
fv = 1 ./ denom;