function u0 = pdex1ic(x)
u0 = sin(pi*x);