This work presents the development/assembly of a computational system for flow analysis in porous media, fractured and fractured porous media and in combination of both media, considering steady or transient states under saturated and unsaturated conditions. The system comprehends four computational programs, three of them of specific functions interconnected by Cplusplus programing routines and the last program is an output viewer. FracGen 3D program (Telles, 2006) generates fractures or fracture families in a determinist or probabilistic way. ICEM CFD v.14 program divides the interest domain in sub-domains by means of the element finite mesh generation. FTPF-3D program (Telles, 2006) uses the element finite method to discretize the governing equations in the space domain and the difference finite method for the time domain and for solving the nonlinearity is used the iterative Picard or BFGS method, so that, finally, Pos3D viewer program is answerable by visualization of the results. In the present dissertation five examples were developed, two of them for the validation of this procedure and the three others applied to a typical slope in Rio de Janeiro, which include vertical fractures and relief joints on their slopes. All those studied cases evaluate the influence of the fractures on porous media in terms of pressure and total heads and velocity fields for verifying of the hydraulic behavior of solid masses and eventual instabilities.