Numerical method of slope failure probability based on Bishop model