SoTaInverSpinel

#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 #産生濃度

SoTaInverSpinel

#Source-Sink model
#参考：三浦岳、発生の数理、京都大学学術出版会、2015、4.3
#julia v0.6
using Plots
gr()

dx=0.05 #
dt=0.1 #単位時間
u0=zeros(1/dx) #初期濃度

SoTaInverSpinel

#SOD-model
#参考：三浦岳、発生の数理、京都大学学術出版会、2015、4.4
#julia v0.6
using Plots
gr()

dx=0.05 #
dt=0.1 #単位時間
Du=0.01
deg=0.5

SoTaInverSpinel

#陰解法
#参考：三浦岳、発生の数理、京都大学学術出版会、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]

SoTaInverSpinel

#1D-Turing pattern
#参考：三浦岳、発生の数理、京都大学学術出版会、2015、6.3
#julia v0.6
using Plots
gr()

dx=0.02
dt=0.01
N=convert(Int,1/dx)

SoTaInverSpinel

#2D-Turing pattern陽解法
#参考：三浦岳、発生の数理、京都大学学術出版会、2015、6.6
#julia v0.6
using Plots
gr()

dx=0.04
dt=0.02
N=convert(Int,1/dx)

SoTaInverSpinel

#2D-Turing pattern陰解法
#参考：三浦岳、発生の数理、京都大学学術出版会、2015、6.6
#julia v0.6

#using FFTW
using Plots
gr(size=(320,300),titlefont=font("sans-serif",9))

dx=0.04
N=convert(Int,1/dx) #gridSize=25

SoTaInverSpinel

#2D-Turing pattern陰解法＋水玉模様
#参考：三浦岳、発生の数理、京都大学学術出版会、2015、6.8
#julia v0.6

#using FFTW
using Plots
gr()

dx=0.04
N=convert(Int,1/dx) #gridSize=25

SoTaInverSpinel

#BZ-reaction
#reference:http://discovery.ucl.ac.uk/17241/1/17241.pdf
#julia v0.6

using Plots
gr(size=(320,300))

h=400
w=400

SoTaInverSpinel

#Swarm Oscillators
#三浦岳、発生の数理、京都大学学術出版会、2015、あとがき
#Julia v1.0

using LinearAlgebra
using Plots
gr(size=(340,300),titlefont=font("sans-serif",9))

const n=50
const L=25