# encoding: utf-8

from scipy.optimize import fsolve

from numpy import *

beita = 1

gama = 1

#以下为Mohr-Columb准则计算

# def f(x):

# OH = (0.33/(1-0.33)+beita)*(14.12-0.73*4.68)+0.73*4.68

# Oh = (0.33/(1-0.33)+gama)*(14.12-0.73*4.68)+0.73*4.68

# Or = 0.946*x+0.25272

# Osita = 3*OH-Oh-x*0.683552239-1.481

# Oz = 12.639+0.66*(OH-Oh)+0.316447761*x

# x = float(x)

# return Osita-3.4164-(Or-3.4164)*4.4835+10.29

# result1 = fsolve(f, 7)

# print '摩尔-库伦准则计算结果：', result1*10/6/9.8, 'g/cm^3'

# 以下为Druker-Prager准则计算

def g(y):

OH = (0.33/(1-0.33)+beita)*(14.12-0.73*4.68)+