{VERSION 6 0 "IBM INTEL NT" "6.0" } {USTYLETAB {CSTYLE "Maple Input" -1 0 "Courier" 0 1 255 0 0 1 0 1 0 0 1 0 0 0 0 1 }{PSTYLE "Normal" -1 0 1 {CSTYLE "" -1 -1 "Times" 1 12 0 0 0 1 2 2 2 2 2 2 1 1 1 1 }1 1 0 0 0 0 1 0 1 0 2 2 0 1 }} {SECT 0 {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 20 "restart:with(plots); " }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 48 "Exact traveling wave solution of Fisher equation" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "u:=x->1/(1+e xp(a*x))^2;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 48 "d*diff(u(x), x,x)+c*diff(u(x),x)+k*u(x)*(1-u(x));" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "simplify(%);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "v:=(x,t)->1/(1+exp(-x-5*t))^2;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 52 "diff(v(x,t),t)-diff(v(x,t),x,x)-6*v(x,t)*(1-v(x,t)); " }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 12 "simplify(%);" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 47 "animate(v(x,t), x=-20..20, t =-2..2,frames=500);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "w:=( x,t)->1/(1+exp(x+5*t))^2;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 47 "animate(w(x,t), x=-20..20, t=-2..2,frames=500);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 64 "diff(w(x,t),t)-diff(w(x,t),x,x)-6*w(x,t)* (1-w(x,t));simplify(%);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 64 " diff(w(x,t),t)+diff(w(x,t),x,x)+6*w(x,t)*(1-w(x,t));simplify(%);" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 20 "restart:with(plots);" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 18 "Try another power?" }}{PARA 0 "> " 0 "" {MPLTEXT 1 0 23 " u:=x->1/(1+exp(a*x))^3;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 60 " d*diff(u(x),x,x)+c*diff(u(x),x)+k*u(x)*(1-u(x));simplify(%);" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 23 "v:=x->1/(1+exp(a*x^2));" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 60 "d*diff(v(x),x,x)+c*diff(v(x) ,x)+k*v(x)*(1-v(x));simplify(%);" }}}{EXCHG {PARA 0 "" 0 "" {TEXT -1 45 "Convergence of solution to the traveling wave" }}{PARA 0 "> " 0 " " {MPLTEXT 1 0 23 "restart:with(PDEtools);" }}}{EXCHG {PARA 0 "> " 0 " " {MPLTEXT 1 0 56 "eq:=diff(u(t,x),t)=diff(u(t,x),x,x)+6*u(t,x)*(1-u(t ,x));" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 52 "IBC:=\{u(0,x)=Heav iside(x-4), u(t,-20)=0, u(t,20)=1\};" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 47 "pds:=pdsolve(eq, IBC, numeric, spacestep=1/40);" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 29 "pds:-animate(t=4,frames=200) ;" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 79 "IBC1:=\{u(0,x)=3*Heavi side(x-2)*sin(x)-2*Heaviside(x-4), u(t,-20)=0, u(t,20)=1\};" }}} {EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 49 "pds1:=pdsolve(eq, IBC1, nume ric, spacestep=1/40);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 30 "pd s1:-animate(t=4,frames=200);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 79 "IBC2:=\{u(0,x)=3*Heaviside(x-2)*sin(x)-2*Heaviside(x-4), u(t,-20 )=0, u(t,20)=1\};" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 46 "pds2:= pdsolve(eq, IBC2, numeric, spacestep=1);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 32 "pds2:-animate(t=0.1,frames=100);" }}}{EXCHG {PARA 0 " > " 0 "" {MPLTEXT 1 0 48 "IBC3:=\{u(0,x)=(x+20)/40, u(t,-20)=0, u(t,20 )=1\};" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 46 "pds3:=pdsolve(eq, IBC3, numeric, spacestep=1);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 31 "pds3:-animate(t=0.01,frames=3);" }}}{EXCHG {PARA 0 "> " 0 "" {MPLTEXT 1 0 0 "" }}}}{MARK "33" 0 }{VIEWOPTS 1 1 0 1 1 1803 1 1 1 1 } {PAGENUMBERS 0 1 2 33 1 1 }