Integrating sin(x) using Simpson Rule

Simpson.py

#Integracion Simpson en funcion seno
#
from math import *
def f(x):
  return sin(x)

def simpson_rule(a,b):
  c=(a+b)/2.0
  h=abs(b-a)/2.0
  return h*(f(a)+4.0*f(c)+f(b))/3.0

val_a = input("Ingrese el valor a: ")
val_b = input("Ingrese el valor b: ")
res = simpson_rule(val_a, val_b)
print "Valor aproximado de la integracion: "+str(res)