/
porousfisher_odefn.m
65 lines (58 loc) · 2.58 KB
/
porousfisher_odefn.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
% Nondimensionalised Porous-Fisher model applied to tissue engineering pore infilling.
%
% Authors: Pascal R Buenzli, Matthew J Simpson, 2019-2020
function dudt = odefn(t,u,N,alpha2,dx)
u=reshape(u,N,N);
dudt1=zeros(N,N);
for i=2:N-1
for j=2:N-1
dudt1(i,j)=alpha2*((u(i+1,j)+u(i,j))*(u(i+1,j)-u(i,j))-(u(i,j)+u(i-1,j))*(u(i,j)-u(i-1,j)))/(2.*dx^2) ...
+alpha2*((u(i,j+1)+u(i,j))*(u(i,j+1)-u(i,j))-(u(i,j)+u(i,j-1))*(u(i,j)-u(i,j-1)))/(2.*dx^2) ...
+u(i,j)*(1.0-u(i,j));
end
end
i=1;
for j=2:N-1
dudt1(i,j)=alpha2*((u(i,j)+u(i+1,j))*(u(i+1,j)-u(i,j))-(u(i,j)+u(i+1,j))*(u(i+1,j)-u(i,j)))/(2.*dx^2) ...
+alpha2*((u(i,j+1)+u(i,j))*(u(i,j+1)-u(i,j))-(u(i,j)+u(i,j-1))*(u(i,j)-u(i,j-1)))/(2.*dx^2) ...
+u(i,j)*(1.0-u(i,j));
end
i=N;
for j=2:N-1
dudt1(i,j)=alpha2*((u(i,j)+u(i-1,j))*(u(i,j)-u(i-1,j))-(u(i,j)+u(i-1,j))*(u(i,j)-u(i-1,j)))/(2.*dx^2) ...
+alpha2*((u(i,j+1)+u(i,j))*(u(i,j+1)-u(i,j))-(u(i,j)+u(i,j-1))*(u(i,j)-u(i,j-1)))/(2.*dx^2) ...
+u(i,j)*(1.0-u(i,j));
end
j=1;
for i=2:N-1
dudt1(i,j)=alpha2*((u(i+1,j)+u(i,j))*(u(i+1,j)-u(i,j))-(u(i,j)+u(i-1,j))*(u(i,j)-u(i-1,j)))/(2.*dx^2) ...
+alpha2*((u(i,j)+u(i,j+1))*(u(i,j+1)-u(i,j))-(u(i,j)+u(i,j+1))*(u(i,j+1)-u(i,j)))/(2.*dx^2) ...
+u(i,j)*(1.0-u(i,j));
end
j=N;
for i=2:N-1
dudt1(i,j)=alpha2*((u(i+1,j)+u(i,j))*(u(i+1,j)-u(i,j))-(u(i,j)+u(i-1,j))*(u(i,j)-u(i-1,j)))/(2.*dx^2) ...
+alpha2*((u(i,j)+u(i,j-1))*(u(i,j)-u(i,j-1))-(u(i,j)+u(i,j-1))*(u(i,j)-u(i,j-1)))/(2.*dx^2) ...
+u(i,j)*(1.0-u(i,j));
end
i=1;
j=1;
dudt1(i,j)=alpha2*((u(i,j)+u(i+1,j))*(u(i+1,j)-u(i,j))-(u(i,j)+u(i+1,j))*(u(i+1,j)-u(i,j)))/(2.*dx^2) ...
+alpha2*((u(i,j)+u(i,j+1))*(u(i,j+1)-u(i,j))-(u(i,j)+u(i,j+1))*(u(i,j+1)-u(i,j)))/(2.*dx^2) ...
+u(i,j)*(1.0-u(i,j));
i=N;
j=1;
dudt1(i,j)=alpha2*((u(i,j)+u(i-1,j))*(u(i,j)-u(i-1,j))-(u(i,j)+u(i-1,j))*(u(i,j)-u(i-1,j)))/(2.*dx^2) ...
+alpha2*((u(i,j)+u(i,j+1))*(u(i,j+1)-u(i,j))-(u(i,j)+u(i,j+1))*(u(i,j+1)-u(i,j)))/(2.*dx^2) ...
+u(i,j)*(1.0-u(i,j));
i=1;
j=N;
dudt1(i,j)=alpha2*((u(i,j)+u(i+1,j))*(u(i+1,j)-u(i,j))-(u(i,j)+u(i+1,j))*(u(i+1,j)-u(i,j)))/(2.*dx^2) ...
+alpha2*((u(i,j)+u(i,j-1))*(u(i,j)-u(i,j-1))-(u(i,j)+u(i,j-1))*(u(i,j)-u(i,j-1)))/(2.*dx^2) ...
+u(i,j)*(1.0-u(i,j));
i=N;
j=N;
dudt1(i,j)=alpha2*((u(i,j)+u(i-1,j))*(u(i,j)-u(i-1,j))-(u(i,j)+u(i-1,j))*(u(i,j)-u(i-1,j)))/(2.*dx^2) ...
+alpha2*((u(i,j)+u(i,j-1))*(u(i,j)-u(i,j-1))-(u(i,j)+u(i,j-1))*(u(i,j)-u(i,j-1)))/(2.*dx^2) ...
+u(i,j)*(1.0-u(i,j));
dudt=reshape(dudt1,[],1);