Here is one way to do the function, assuming you are looking for a numerical output, and that you pass it t as well as n:
n = ...
t = ...
output = fn(n,t);
function f = fn(n,t)
c = (-1)^(n-1)*t^(2-2*n)/(1-t^-4);
S = 0;
for k = 0:n-1
tm = t^(-4*n*k);
P = 1;
for i = 0:k
P = P*(1-t^(4*i-4*n));
end
S = S + tm*P;
end
f = c*S;
end
Best Answer