The finite element numerical simulation in classical response surface method was substituted by iterative calculation of implicit stability coefficient of slope in Sarma model.