In this work we performed a numerical analysis of the the outer region of the Pluto-Charon system by the insertion of a sample of test particles initially in P-type prograde and retrograde orbits, in the system formed by Pluto, Charon, Nix and Hydra. These numerical integrations generated diagrams of semi-major axis versus eccentricity which define regions of particles in stable orbits and regions of collision and escape. In the vicinity of the satellites Nix and Hydra were identified chaotic regions, where particles located in this region have their eccentricities and semi-major axis increased provoking an ejection or collision with a massive body of the system. A set of particles remained in regions near the orbits of Nix and Hydra, possibly coorbitais with them. For both cases, prograde and retrograde, the "stable" region is larger in the outer region of the system, after Hydra's orbit, depending on the value of eccentricity. Numerical simulations were also performed by inserting some massive hypothetical satellites beyond the Charon's orbit and the effects on the orbits of Nix and Hydra were analyzed. A preliminary numerical study of the effects of the solar radiation force on a sample of particles with radii of 1μm, 3μm, 5μm e 10μm was performed. This study showed that particles under the effects of the Poynting-Robertson drag decay on a time scale between 1.45 × 10 6 years (particles of 1μm in radius) and 1.45×10 7 years (particles of 10μm in radius), while the radiation pressure caused variations of the eccentricities of the particles causing in some cases collisions with the planet.