function F= FO_d(x, alfa, gama, csi, phi, mu)
%% defining vector x of variables
%% x(1)=dt
%% x(2)=d
%% x(3)=dstar

%%%% giving initial values for the variables

%%% define the function

F=[x(1)*(1/alfa-(phi+gama)*x(1)-csi)+x(2)*(1/alfa-2*(phi+gama)*x(1)-csi);
    x(1)*(1/alfa-(phi+gama)*x(1)-csi-mu)+(1/alfa-2*(phi+gama)*x(1)-csi-mu)*x(3);
    (((x(3)-x(2)))/x(3))*(mu/(phi*alfa))-x(1)*(1/alfa-(phi+gama)*x(1)-csi)];
end