This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
#1D-Turing pattern | |
#参考:三浦岳、発生の数理、京都大学学術出版会、2015、6.3 | |
#julia v0.6 | |
using Plots | |
gr() | |
dx=0.02 | |
dt=0.01 | |
N=convert(Int,1/dx) |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
#陰解法 | |
#参考:三浦岳、発生の数理、京都大学学術出版会、2015、4.5 | |
#julia v0.6 | |
using Plots | |
gr() | |
A=circshift(eye(10),(0,1))+circshift(eye(10),(0,-1))-2*eye(10) | |
dt=0.05 | |
initialDistribution=[0,0,0,1,1,1,0,0,0,0] |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
#SOD-model | |
#参考:三浦岳、発生の数理、京都大学学術出版会、2015、4.4 | |
#julia v0.6 | |
using Plots | |
gr() | |
dx=0.05 # | |
dt=0.1 #単位時間 | |
Du=0.01 | |
deg=0.5 |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
#Source-Sink model | |
#参考:三浦岳、発生の数理、京都大学学術出版会、2015、4.3 | |
#julia v0.6 | |
using Plots | |
gr() | |
dx=0.05 # | |
dt=0.1 #単位時間 | |
u0=zeros(1/dx) #初期濃度 |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
#Source-only model | |
#参考:三浦岳、発生の数理、京都大学学術出版会、2015、p36-43 | |
#julia v0.6 | |
using Plots | |
gr() | |
dx=0.05 | |
dt=0.1 | |
u0=zeros(1/dx) #初期濃度 | |
up=zeros(1/dx);up[1]=1.0 #産生濃度 |
NewerOlder