#Las variables son n: número total de pagos e i: tasa de interés por periodo de conversión
n,i=var("n,i")
#El monto total S de la anualidad está dado por S=R*s(n,i)
s(n,i)=((1+i)**n-1)/i
show("s(n,i)=", s(n,i))

#El valor presente A de la anualidad está dado por A=R*a(n,i)
a(n,i)=s(n,i)/(1+i)**n
show("a(n,i)=", a(n,i))

#precaucion: sólo modifique estos parametros
#Inicializamos nuestros valores 
#r=tasa de interés anual, N=frecuencia de conversión, T:plazo en años 
A=2939
r=0.9017
N=6
T=1
############
i=r/N; n=N*T
print "i=",i
print "n=",n

#También lo podemos hacer directamente usando el valor presente de la anualidad
R = A/a(n,i)
print a(n,i)
print "R=", R
print A*(1.025)**6

#Construimos una tabla de amortización
#Primero creamos una tabla vacía
amortizacion = []
#Ahora, creamos la tabla; vea las notas de clase para la explicación completa
Ak=A
for k in range(1,n+1):
    #esta linea crea los valores para cada periodo y los anexa a la tabla
    amortizacion.append([k, round(Ak,2), round(Ak*i,2), round(R,2), round(R-Ak*i,2)])
    #nuevo capital insulto al inicio del siguiente periodo
    Ak = Ak - (R-Ak*i)
#Calculamos los totales
#TIV = total de intereses vencidos
TIV = sum([amortizacion[k][2] for k in range(n)])
#TP = total de pagos; debe igualar el monto de la anualidad
TP = sum([amortizacion[k][3] for k in range(n)])
#TCP = total del capital pagado al final de cada periodo; verifique que TCP+TIV=TP
TCP = sum([amortizacion[k][4] for k in range(n)])
#finalmente, anexamos estos valores a la tabla de amortizacion
amortizacion.append(["Totales", " " ,TIV, TP, TCP])

#imprimimos nuestra tabla de amortización en el formato "table"
print table(amortizacion)

#Calculamos el monto del pago periódico R usando el monto de la anualidad
S = A*(1+i)**n; print S
R = S/s(n,i)
print R