function y = f2(x,t)

global K C

 y= (C+K*pi*pi)*sin(pi*x);