Matlab PDE boundry condition confusion.

i have a pde where i want to be able to manipulate variables to see how the function acts. The problem is that i don't know how to set boundry conditions or g=graph the square of the function. The function should be: (k is a constant)

DuDt = k*(DuDx)/Dx + (x^2)*u

with an initial condition:(b and s are constants)

u0=e^[(x-b)^2/s]

and boundary conditions

u(-5,t)=0 u(5,t)=0

any help would be very much appreciated