function max_len_seq_sine
input Real weights[:] = {1, 1, 1, 1, 1, 1, 1};
input Integer harmonics[size(weights, 1)] = {1, 2, 4, 8, 16, 32, 64};
input Real bias = 0 "Bias to control middle value";
output Real y[integer(max(harmonics)^2)];
end max_len_seq_sine;