ciie.m

This software gave volume MoM solution by the CG-FFT method

function ciie=ciie(x,y,z,a,b,c,ts,te,ps,pe,x0,y0,z0);

ciie=0;
r=(x0-x)^2/a^2+(y0-y)^2/b^2+(z0-z)^2/c^2; 
theta=180/pi*acos((z0-z)/sqrt((x0-x)^2+(y0-y)^2+(z0-z)^2));
phi=180/pi*atan2((y0-y),(x0-x));
if (phi<0), phi=phi+360;  end;

if (r<=1),
   if (theta>=ts),
      if (theta<=te),
         if (phi>=ps),
             if (phi<=pe),
                ciie=1;
             end;
          end;
       end;
   end;
end;

end;