function [dy] = ecdiflaser(t,y)
%Ecuaciones del láser para el modelo de Lorenz-Haken
%   la función devuelve tres valores: F',P' y d'. 
global s b r

F = y(1);
P = y(2);
d = y(3);

dF=s*(P-F);
dP=-P+d*F;
dd=b*(r-d-F*P);

dy = [dF; dP; dd];
end