function [ N ] = funcao_nivel( PAR, t )
% [ N ] = funcao_nivel( PAR, t );

Nm = PAR(1);
AA = PAR(2);
fA = PAR(3);
AS = PAR(4);
fS = PAR(5);

N = Nm + ...
    AA*cos(2*pi*t/365 + fA ) + ...
    AS*cos(2*pi*t/(11*365) + fS );