%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
        %
        %    alf.m
        %
        %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
        function x1 = alf(t)
        global cglob
        x1(1)=cglob(1)*(t) - cglob(4) - cglob(3)*sin( cglob(2)*(t+cglob(5)));