u1_anal=function(x,t,k,D,c){ # # Function u1_anal computes the analytical solution for u1(x,t) # z=x-c*t; u1a=1/(1+exp(-c*z/D)); # # Return solution return(c(u1a)); }