a0 = 1 / 4;
x = 0:0.01:4;
f10 = [];
f50 = [];
f100 = [];

function [a]=An(n)
    //a = (%pi * n * (t - 2) * sin(%pi * n * t / 2) + 2 * cos(%pi * n * t / 2)) / (%pi ^ 2 * n ^ 2) + 2 * sin(%pi * n * t / 2) / (%pi * n);
    a = (2 * cos(%pi * n) / ((%pi * n) ^ 2) - 2 / ((%pi * n) ^ 2) + 2 * sin(2 * %pi * n) / (%pi * n) - 2 * sin(%pi * n) / (%pi * n)) / 2;
endfunction

function [b]=Bn(n)
    b = (2 * sin(%pi * n) / ((%pi * n) ^ 2) - 2 / (%pi * n) - 2 * cos(2 * %pi * n) / (%pi * n) + 2 * cos(%pi * n) / (%pi * n)) / 2;
endfunction

for i = 0:0.01:4
    t = i
    f10($ + 1) = a0
    f50($ + 1) = a0
    f100($ + 1) = a0

    for j = 1:100
        if (j <= 10)
            f10($) = f10($) + An(j) * cos(%pi * t * j / 2) + Bn(j) * sin(%pi * t * j / 2)
        end
        if (j <= 25)
            f50($) = f50($) + An(j) * cos(%pi * t * j / 2) + Bn(j) * sin(%pi * t * j / 2)
        end
        if (j <= 100)
            f100($) = f100($) + An(j) * cos(%pi * t * j / 2) + Bn(j) * sin(%pi * t * j / 2)
        end
    end
end

plot2d(x,f10)
plot2d(x,f100)