Write a Python program to differentiate the following equation:
y = 5*x^2 + 3*x + 1
Not applicable
def differentiate(x):
  return 10*x + 3;

# Test
x = 2 
print("The derivative of y = 5*x^2 + 3*x + 1 at x = 2 is:", differentiate(x)) # Output: 17