You should use "exp(k*p)" instead of plain old "e".  Also, you have several "D" and that should be made local.  I got a plot after fixing those syntax errors but I can't tell you if it is right!

@mmcdara Thank you all for your great insight!