Yes, your worksheet shows the same plot as I got when I corrected exp(k*p) issue and the "local D" issue.

Just as a aside, you include the package "linalg" which has been "deprecated" and it is not recommended by

the "Big" guns; you should use "LinearAlgebra" instead.

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!