function y = f(x)

% for step function initial data:
y  = x<=(pi-1);