def partial_difference_quotient(f, v, i, h):
	"""
	Compute the ith partial difference quotient of f at v
	"""
	w = [v_j + (h if j== i else 0) 		# add h to just the ith element of h
			for j, v_j in enumerate(v)] 

	return (f(w) - f(v)) / h