I present numerical simulations of intensity and degree of polarization of light reflected by Earth-like exoplanets. The results are presented as a function of wavelength, and for a few different phase angles and a few different points on the planet. At this stage the aim is to show the working code and test a few different set ups of the star-planet system in order to find preferable configurations for observations. Not surprisingly, phase angle 90◦ shows the largest degree of polarization. For beneficial wavelength regions, visual light shows a larger overall degree of polarization, while NIR shows very clear absorption patterns in the degree of polarization, making detection of the atmospheric composition possible.