function Yex = Yexfct(t,x)

global alpha
l
for 
if (x > (2*t))
    Yex= cos (2*pi*x-2*t); %= Y0(x-2t)
else
    Yex = alpha;
end