function yp = whale(t,y,flag)
%                Parameters
        r1 = .05;
        r2 = .08;
        alf1 = .00000001;
        alf2 = .00000001;
        kk1 = 150000;
        kk2 = 400000;
        gam1 = kk2.*alf1;
        gam2 = kk1.*alf2;
%             
        yp1= r1.*y(1).*(1-y(1)) - gam1.*y(1).*y(2);  % equation for dy1/dt
        yp2= r2.*y(2).*(1-y(2)) - gam2.*y(1).*y(2);  % equation for dy2/dt
        yp = [yp1;yp2];