function y = f(x,t)

global beta c

 y= exp(-t)* ( -sin(pi*x) + beta*pi*pi*sin(pi*x) +c*sin(pi*x) );