function y = kernel(tau)
global TOUT
y = min(TOUT,tau);