Navigation Menu

Skip to content

Instantly share code, notes, and snippets.

View SoTaInverSpinel's full-sized avatar

STInverSpinel SoTaInverSpinel

View GitHub Profile
#1D-Turing pattern
#参考:三浦岳、発生の数理、京都大学学術出版会、2015、6.3
#julia v0.6
using Plots
gr()
dx=0.02
dt=0.01
N=convert(Int,1/dx)
#陰解法
#参考:三浦岳、発生の数理、京都大学学術出版会、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]
#SOD-model
#参考:三浦岳、発生の数理、京都大学学術出版会、2015、4.4
#julia v0.6
using Plots
gr()
dx=0.05 #
dt=0.1 #単位時間
Du=0.01
deg=0.5
#Source-Sink model
#参考:三浦岳、発生の数理、京都大学学術出版会、2015、4.3
#julia v0.6
using Plots
gr()
dx=0.05 #
dt=0.1 #単位時間
u0=zeros(1/dx) #初期濃度
#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 #産生濃度