In this paper heat conduction equation for an eccentric spherical annulus with the inner surface kept at a constant temperature and the outer surface subjected to convection is solved analytically. Eccentric problem domain is first transformed into a concentric domain via formulating the problem in bispherical co-ordinate system. Since an analytical Green's function for the heat conduction equation in bispherical co-ordinate for an eccentric sphere subject to boundary condition of third type can not be found, an analytical Green's function obtained for Dirichlet boundary condition is employed in the solution. Utilizing this Green's function yields a boundary integral equation for the unknown normal derivative of the surface temperature distribution. The resulting boundary integral equation is solved analytically using method of moments. The method has been applied to heat generating eccentric spherical annuli and results are compared to the simulation results of FLUENT CFD code. A very good agreement was observed in temperature distribution computations for various geometrical configurations and a wide range of Biot number. Variation of heat dissipation with radii and eccentricity ratios are studied and a very good agreement with FLUENT has been observed.