syms y;
syms x;
syms b;

syms o;

o = (1/pi)*(sin(n*w*x))/n;

o = serieFourier(100, 1, 0.5, o);

fplot(o, [0, 0.25]);


function y = serieFourier(n, f, b, s)
    y = 0;
    w = 2*pi*f;
    for i=1:n
       y = y + s;
    end
    y = b + y;
end