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
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
|
clc
clear
for i=1:60
for j=1:60
A(i,j)=0
end
end
for i=2:9
A(i,i)=-1
A(i,10+i)=1
end
for i=1:10:51
A(i,i)=1
end
for i=10:10:60
A(i,i)=1
end
for i=12:19
A(i,i-10)=1
A(i,i-1)=1
A(i,i)=-4
A(i,i+1)=1
A(i,i+10)=1
end
for i=22:29
A(i,i-10)=1
A(i,i-1)=1
A(i,i)=-4
A(i,i+1)=1
A(i,i+10)=1
end
for i=32:39
A(i,i-10)=1
A(i,i-1)=1
A(i,i)=-4
A(i,i+1)=1
A(i,i+10)=1
end
for i=42:49
A(i,i-10)=1
A(i,i-1)=1
A(i,i)=-4
A(i,i+1)=1
A(i,i+10)=1
end
for i=52:59
A(i,i-10)=-1
A(i,i)=(1+(2*0.073)/0.16)
end
disp(A)
for i=2:9
B(i,1)=30
end
for i=1:10:51
B(i,1)=20
end
for i=10:10:60
B(i,1)=20
end
for i=12:19
B(i,1)=0.6/(20*0.16)
end
for i=22:29
B(i,1)=0.6/(20*0.16)
end
for i=32:39
B(i,1)=0.6/(20*0.16)
end
for i=42:49
B(i,1)=0.6/(20*0.16)
end
for i=52:59
B(i,1)=25*(2*0.073)/0.16
end
X=A\B
for i=1:10:60
T=[ X(i,1), X(i+1,1),X(i+2, 1), X(i+3,1), X(i+4,1),X(i+5,1),X(i+6,1),X(i+7,1),X(i+8,1),X(i+9,1)]
disp(T)
end
|
