%
%  influent disturbance function
%
function y = inffunc( t, p )
%
y = p(1) + p(2) * ( sin( ( 2 * pi * t / 24 ) ...
                                 + ( pi * p(4) ) ) ...
         + p(3) * sin( ( 2 * pi * t / 12 ) ...
		                 + ( pi * ( p(4) + p(5) ) ) ) );
%