% Aligned coaxial disks
function Fij = DiskView(ri,rj,h)
  Ri=ri./h; Rj=rj./h; S=1+(1+Rj.^2)./Ri.^2;
  Fij = (S-sqrt(S.^2-4*(Rj./Ri).^2))/2;
end
