By computer simulation the growing process of porous silicon under p(-)-Si and p(+)-Si anodization in HF solution is studied. The model of electrochemical etching of p(+)-Si includes the relief selective mechanism, which allows one to establish the relationship between anodization conditions (current density, HF concentration, temperature and doping level) and the topological characteristics of porous silicon (PS). The simulation of p(-)-Si dissolution is based on the model of diffusion limited aggregation (DLA), taking into account the thermal generation of holes and the quantum confinement effect. The various morphology of simulated PS structures exhibits a close resemblance to that of experimental ones formed in p(+)-Si and p(-)-Si wafers. For simulated p(-)-Si-based PS layers the porosity profiles and fractal dimension are calculated. It is shown that PS in p(-)-Si is multifractal with fractal dimension varying monotonously from 0.1 to 3 with size increase.