% Hemi-cylindrical shield and environment obstructed by cylinder
function Fij = HemiCylPastFullCylView(ri,rj)
  R=ri./rj;
  Fij = R.*(2/pi*(sqrt(R.^-2-1)+asin(R))-1);
end
