Today, computational acoustics is one of major branch of underwater acoustics that play important role in the advanced laboratories in the world. It is important due to high financial cost of operational tests in the water. There exist many numerical methods in this field that among them only, ray tracing, normal modes, fast field programming and parabolic equation method have gained more attention. PE method is a suitable method in low frequencies and shallow water as well as valid for all environments. In this paper, Numerical implementation of PE method based on c++ and MATLAB routines is given over some of well known standard test problems. In the implementation, according to both mathematically and physically complexities of problem, we chose general manner for some of important parameters such as depth and frequency of source, range of receiver, and the whole depth of simulation. Some results show that we have surface and bottom reflections from 4500 meters in range. TL error has a value about 18.4 percent. As increasing source strength from 200 dB, we can eliminate noise effects. TL in stations 1 to 4 in Caspian sea is shown as well as ray simulation in graphs 2, 6, 8, 10.