In this thesis, we study bounds for zeros of harmonic functions $p(z)-overline{z}$ and $r(z)-overline{z}$, where $p$ is a polynomial and $r$ is a rational function. A surprising applications of results obtained in this study is that we are able to bound the number of gravitational lensed images in some astronomical models.