ua_1=function(x,t){ # # Function ua_1 computes the analytical solution for # the Fitzhugh-Nagumo equation # ua=1/(1+exp((1/2^0.5)*x/D^0.5+(a-1/2)*t)); return(c(ua)); }