import time
from machine import Pin
import math

print(" Serie de Fourier - Onda Cuadrada ")

def graficador_cuadrada(nmax):  # numero de armonicos
    xmin = 0
    xmax = 60
    fx = 0.0
    delta = 0.5
    T = 60
    a0 = 0
    x = xmin
    while x < xmax:
        sum_val = 0
        w0 = 2 * math.pi / T
        n = 1
        while n <= nmax:
            an = 0.0
            bn = (2 / (n * math.pi)) * (1 - math.pow(-1, n))
            sum_val += (an * math.cos(n * w0 * x) + bn * math.sin(n * w0 * x))
            n += 1
        fx = a0 / 2 + sum_val
        print(fx)
        x += delta

while True:
    for narmonicos in range(0, 25, 1):
        graficador_cuadrada(narmonicos)