def surf_eos(V):
return feos.vinet(V, 1.2, 1.2, 1.2, 1.2) # , convert_to_sympy=sympy_)
def f4(V):
a = 3 * V ** 3 - 8 * V - 7 * V ** 2 + 2 * V ** 4
print "aaaaaaaaaaaaaaaa|", a, "||", type(a)
return a
def f4der(V):
return 8 * V ** 3 + 9 * V ** 2 - 14 * V - 8
f4_ = 3 * V ** 3 - 8 * V - 7 * V ** 2 + 2 * V ** 4
f4der_ = 8 * V ** 3 + 9 * V ** 2 - 14 * V - 8
print feos.vinet(V, 1.2, 1.2, 1.2, 1.2, convert_to_sympy=True)
print minimize_scalar(surf_eos).x
print minimize_scalar(f2).x
print "<<"
print minimize_scalar(f3, bounds=(-4, -1), method="bounded").x
print "vvv bei -1.79 und 1.15 vvvvvvvvv"
print minimize_scalar(f4, bounds=(-3, -1), method="bounded").x
print minimize_scalar(f4, bounds=(0, 1000), method="bounded").x
print "DER bei 0.47 und xxx vvvvvvvvv"
print minimize_scalar(f4der, bounds=(-3, 10), method="bounded").x
print "bfgs:vvv"
from scipy import optimize
print optimize.fmin_bfgs(f4, 7)
