Python解非线性方程

# 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)+0.73*4.68
    Oh = (0.33/(1-0.33)+gama)*(14.12-0.73*4.68)+0.73*4.68
    Or = 0.946*y+0.25272
    Osita = 3*OH-Oh-y*0.683552239-1.481
    Oz = 12.639+0.66*(OH-Oh)+0.316447761*y
    y = float(y)
    return (math.sqrt(1/6*((Or-Osita)**2&