python 实现差商

差商计算中没有加入O(h),数学还在继续加强中，理解的可能有点不到位，希望理解有误的地方大家能给出宝贵建议

#coding=utf-8
from sympy import *
#向前差商的函数
def forward():
    f1.subs(x,0.5)
    g1 = (f.subs(x,0.5 + h) - f.subs(x,0.5))/ h
    return g1

#向后差商的函数
def backword():
    g1 = (f.subs(x,0.5) - f.subs(x,0.5 - h)) / h
    return g1

#中心差商的函数
def center():
    g1 = (f.subs(x,0.5 + h) - f.subs(x,0.5-h)) / (2 * h)
    return g1


h = input("请输入h的值:")
#print type(h)
x=Symbol("x")
f = -0.1 * x**4 + -0.15 * x**3 - 0.5 * x**2 - 0.25 * x +1.2
f1 = diff(f,x)
print f1
print "向前差商结果为："
print forward()
print "向后差商结果为："
print backword()
print "中心差商的结果为："
print center()