Python - ¿Por qué marca error de sintaxis

#====================================

#MacLaurin.py

#====================================

#El siguiente código es un modulo perteneciente al pruebademclaurin.py

#Este modulo se encarga de reflejar los valores de la secuencia de McLaurin de un valor dado.

#Este modulo hace uso de funciones

import math #Importamos la biblioteca de matematicas, esto con el fín de utilizar los modulos de la biblioteca

def sen(x,n): #Mediante esta función es posible obtener el valor de seno de la serie de McLaurin del valor dado por el usuario.

	xrad= math.radians(x) #Convertimos el valor ingresado de x en radianes

	senx = 0

	for k in range(1, n):

		term = ((-1)**k)*(xrad**((2*k)+1))/math.factorial((2*k)+1)

		senx += term

	return senx

def exp(x,n):

	expx = 1

	for k in range(1, n):

		term = (x**k)/math.factorial(k)

		expx += term

	return expx

def tabulafuncionsen(x,nMin,nMax,incN):

	print("{:}			sen({:f})".format("n",x)) #Imprimimos el encabezado de la tabla

	print("_____________________________________________________________") #Imprimimos la división de la tabla

	for n in range(nMin,nMax+1,incN):

			print("{:d}			sen({:f})".format(n,sen(x,n))

def tabulafuncionexp(x,nMin,nMax,incN):

	print()

	print("{:}			exp({:f})".format("n",x)) #Imprimimos los titulos

	print("_____________________________________________________________") #Imprimimos la división de la tabla

	for n in range(nMin,nMax+1,incN):

		print("{:d}			{:f}".format(n,exp(x,n))

def tabulafuncion(x,nMin,nMax,incN,operacion):

	#Determinamos primero el tipo de función

	if op= 'Seno' or operacion == 's' or operacion == 'S' or operacion == 'sin' or operacion == 'Sin': #En caso de que la función sea Sen(x)

		tabulafuncionsen(x,nMin,nMax,incN)

	elif operacion == 'E' or operacion == 'Exp' or operacion == 'e' or operacion == 'exp' or operacion == 'Exponencial':

	#En caso de que sea Exp(x)

		tabulafuncionexp(x,nMin,nMax,incN)

#Prueba del modulo

if __name__ == "__main__":

	print("Tabula seno de 100 grados, con un limite minimo de 0 y máximo de 10, los culaes van incrementando de 1 en 1")

	tabulafuncion(100,0,10,1,"s")

print("{:d}			sen({:f})".format(n,sen(x,n))

print("{:d}			sen({:f})".format(n,sen(x,n)))