function ggprime=pderhs(t,gg,opt,k,dx,A) ggprime=(k/dx^2)*A*gg;