area-metric-gravity/field-equations/scalar_stationary.maple at master · nilsalex/area-metric-gravity · GitHub

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
with(LinearAlgebra);

s2 := -s1/4;
s5 := -s4/3;
s7 := -s6/3;
s8 := s1 + 4*s4;
s9 := s1/4 + s4;
s10 := (3*s1)/4 + 3*s4;
s12 := -s1/3 - s11 - (4*s4)/3;
s15 := -s14;
s17 := 4*s6;
s18 := s6;
s19 := 3*s6;
s20 := s14 - s6;
s21 := -s14 - s6/3;
s22 := s16;
s23 := -s1/4 - s11 - s4;
s24 := s1/4 + s11 + s4;
s25 := -s13;
s26 := (-2*s1)/3;
s27 := -s1/2;
s28 := s1/6;
s29 := (-3*s1)/4 + s3;
s30 := (-2*s3)/3;
s31 := s1/3;
s32 := (2*s4)/9;
s33 := (2*s6)/9;
s34 := -s1 + (4*s3)/3;
s35 := (-3*s1)/4 + s3;
s36 := (-2*s3)/3;
s38 := (-3*s1)/2 + 2*s3 - s37;
s40 := s1/2 - (2*s3)/3;
s41 := s1/6 + (2*s3)/9 + (2*s4)/3;
s42 := (2*s6)/3;
s43 := -2*s1;
s44 := -3*s1 + 4*s3;
s45 := (4*s1)/3 + (8*s4)/3;
s46 := (8*s6)/3;

mat := Matrix([

[s1, s2, s3, - s5 * b^2, - s7 * b^2],
[s8, s9, s10, - s12 * b^2 + s13, - s15 * b^2 + s16],
[s17, s18, s19, - s21 * b^2 + s22, - s24 * b^2 + s25],
[-s26 * b^2, - s28 * b^2, - s30 * b^2, s32 * b^4, s33 * b^4],
[-s34 * b^2 , - s36 * b^2, - s38 * b^2 + s39, s41 * b^4, s42 * b^4],
[0, -s43 * b^2, -s44 * b^2, s45 * b^4, s46 * b^4]

]);

vec := Vector([0,0,0,0,0,M]);

sol := LinearSolve(mat, vec);