Chemical etching is employed as micromachining manufacturing process to produce micron-size components. As a semiconductor wafer is extremely expensive due to many processing steps involved in the making thereof, the need to critically control the etching end point in an etching process is highly desirable. It was found that not only the etchant and temperature determine the exact anisotropy of etched silicon. The angle between the silicon surface and the mask was also shown to play an important role. In this paper, angular dependence of the etching rate is calculated on the base of the silicon symmetry properties, by means of the interpolation technique using experimentally obtained values of the principal <100>, <110>, <111> directions in KOH solutions. The calculations are performed using an extension of the sparse field method for solving three dimensional (3D) level set equations that describe the morphological surface evolution during etching process. The analysis of the obtained results confirm that regardless of the initial shape the profile evolution ends with the crystal form composed of the fastest etching planes, {110} in our model.