function y = u0tilde(x,a)
y = zeros(size(x));
indici = (a <= x);
y(indici) = sin(x(indici)).^50;
