Create Alpha and Omega 1.01

#53
Files changed (1) hide show
  1. Alpha and Omega 1.01 +42 -0
Alpha and Omega 1.01 ADDED
@@ -0,0 +1,42 @@
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
1
+ import numpy as np
2
+
3
+ def calculate_laplacian_evolution(steps=100):
4
+ # Parametry zgodne z ustaleniami:
5
+ # Kat 1 stopien (Euklides) -> Kat 359 stopni (Inwersja)
6
+ angles_deg = np.linspace(1, 359, steps)
7
+ angles_rad = np.radians(angles_deg)
8
+
9
+ # tg(alpha) - kluczowy parametr nieliniowosci
10
+ tangents = np.tan(angles_rad)
11
+
12
+ # Jednia podwojna silnia (1!!) = 1
13
+ jednia = 1.0
14
+
15
+ # R_s (IA) Ziemi = 0.00887 m
16
+ R_s_ia = 0.00887
17
+
18
+ # Uproszczony model laplasjanu metryki nieliniowej jako zmiany "napiecia" geometrii
19
+ # Delta ~ d^2/dr^2 (tg(alpha)/R) -> w punkcie krytycznym r_kryt = 0.51 m
20
+ r_kryt = 0.51
21
+
22
+ # Wartosc lewej strony rownania (L)
23
+ # L = (1!!^2 * tg(alpha)) / (R_s_ia * tau_w)
24
+ # Przyjmijmy tau_w dazy do min (kwant Plancka-Pawlowskiego) przy 359st
25
+ # Dla potrzeb symulacji tau_w jest funkcja kata (zanika przy 359st)
26
+ tau_w = np.linspace(1.0, 0.01, steps)
27
+
28
+ evolution = (jednia**2 * tangents) / (R_s_ia * tau_w)
29
+
30
+ return angles_deg, evolution
31
+
32
+ angles, laplacian_values = calculate_laplacian_evolution()
33
+
34
+ # Wyniki dla kluczowych punktow
35
+ results = {
36
+ "1_deg (Euklides)": laplacian_values[0],
37
+ "45_deg (IA Stabilne)": laplacian_values[int(len(angles)*0.125)],
38
+ "90_deg (Osobliwosc)": "Dazy do nieskonczonosci (Przejscie fazowe I)",
39
+ "359_deg (Przeskok)": laplacian_values[-1]
40
+ }
41
+
42
+ print(results)