From c6b68daf98277e9e14695e0b776daa457883af3e Mon Sep 17 00:00:00 2001
From: Indigo5684 <159226326+Indigo5684@users.noreply.github.com>
Date: Tue, 30 Sep 2025 13:19:35 -0500
Subject: [PATCH] Modern Phys Lab Planck's Constant Update

---
 statics/physics/modern-phys-lab/data.tar.xz | Bin 2292 -> 0 bytes
 statics/physics/modern-phys-lab/run.py      |  38 ++++++++++++--------
 2 files changed, 24 insertions(+), 14 deletions(-)
 delete mode 100644 statics/physics/modern-phys-lab/data.tar.xz

diff --git a/statics/physics/modern-phys-lab/data.tar.xz b/statics/physics/modern-phys-lab/data.tar.xz
deleted file mode 100644
index fc9c7b2ce85ad07fd663b8ac843c2f9a0d37ecfe..0000000000000000000000000000000000000000
GIT binary patch
literal 0
HcmV?d00001

literal 2292
zcmV<Q2n+Z9H+ooF000E$*0e?f03iVu0001VFXf})EB^?wT>vo~3lQQ+UR{Z0h<KsS
z7W|`N*Y#+5*>RtJ!CT;JX;<m1dfdFmeZJV@V72f(<3_29u$3z8Ms8hN%|wlk6bNX)
zN;{<5+GA(<=!TI1l`rWN`UrW0kWSb+D(N>>h~w~G93eEJ4gl*QTd=33cx$4)|2e;b
zsHKp+(>?hBMzIi9M{1^Pi?Y!2$;nX!xq+VU#dBh-HxRHFPiwzlSYy{9F-%DF9Y#uM
z&Bw?bxYS$|gTp1kKO+m!LCvW+V(|zZm!R#@dOdgTJ^0vwP9;Lk7yilr=sJL{#nKha
z|L1|_*EdNIUOa_#h_E87wmpLTqdeD2#ci%aIP270sFagUFgmEYWzEEYW2mPUnNdiE
z;@>8~u%kPQ%ugzQuD{H>tGN}xJPbYkDn<?kr{CH|MQJ^LMt6J}^;w6GIOI%9Kc$ov
z9g*AU<#idC#NqKE*IAW%Sbl#2-V31h4rH<V{Ni9y9caR%8kBgMm$W?`b-|9NG6D26
z!jbyZ?yJ`pmdaEC$F-6-3&bHR_vREQ+PqPU56CsY+3NfXzCHU-8I-(wrNKWtZDt>b
zd24s1uYzZ0HV5ylqiQiGZIt|Cz>JM%nRMo!0g;IEI2R_&B%Y@BB9-okwI^u1g`Q3c
z7^(_gxFsCkfgk|Btn9lj^Z9W1PQA-WF>4dR*FA%Ob-_P3=w_HbA4@NJeFh1m)TF!;
zCVB(^DhwOVkQk|tR>qmPr?e=8%H#-vtI$_Cu8ZbHG5T?ycIZ+voMqvSNX08Qs^_JF
zpGA|P`&UL%08x(^uqtj7^AO1`MGA|LEG>X&!C#`Ah+^CeMA=YM)jWd;<E<`qLU>j%
zqn3?h+o5VrKd<TKhopi8b&QM@!>Az$af^&!TdJ5Ib_n+F*w4Tr>Nd_Vogv)b4$z!A
z4yk&&Z@UpIe5vBrcz$(QqTptOwJZ%bVyJL~j2+o+r?2|Z?8!7zG)3!P4PGL;gNB);
z9Ae3IY05rq?AbouX;Y+477PqCjOFhEXRsny{*srKWg)mH{!*@@#^uP|HC}omiG2}j
zS^f?i@dOLAF2nqpD({2hl{@guudHBWj;`CDX*@CcI99@Cm6#GQw4mfqC8_=>m4qcg
z=w1XGCmuuo@}!7KER5tt!yt&vI8+$QHlYO$fzSnuMh5u|YkaOWX4BIzX*UYT4chZ^
zEBm*+$-bK;*QApHQ6wdB-ggI2FZ-2i4%vB&^m3TcbX6KCJVtV2?wUsG9U6oibnY<i
zFZj?fhqByAE%}IUPNOBPN>`x>gRG&do(?1X0kHBeoTd@`LSz5EgS8^L<|&unT|3$^
zhU29Trdwb<hf-INX<KnHnN}?DdLW@|+!W_zup))-c`c{5BiO4!xP6qz+HG4B-U;Da
zKw1MRjDQ=PV0op^(5G(`Kspn<H52nrsZwf$j)QiCm|5RpRW{Anjl%v45rbF~RNek$
zrOUynq61-RlKCSmx1+b%f9<cw&)SO4)e!Tsw^p5#7NCgC=Gz+NUGoymoC)C2wlAy`
zUTJkJesb=$#{`jY_6~wSO}xu{y937FD_!Q~ZpVK@u>|3xEN6ni3`x0^g3dINxfreV
z2XR)V@(rCHc0h<y9Bj{$U(e;l_#N~6+j8nTzcBK7ro#G+%A4I##&Li~ZQbtTC;_J@
zh03B!<h83NMFm6RJ7MN(oiIda+Hn}{C<;6#@5y?|B3;&yLp1jv?m&NNmH}pvnE<=8
zIzAmNWMOFJ9UL3y$(|%_LFmM%Nu=VNs=%-0;FA*3+w$QR6!8Tf;^^h(WONKW8D$<9
zsU~~65~ru?Uj2v{0qEwN5)4jDDP95m_HR?&%nI~+GJ981k4n3S)@=`H6!|c`2Ai$6
z1-@gR-WLqMPVuryI|%P_KaxWLq&`9~{00fSGCo!5RNAMVOkpnzH8Bi4DGywx|7dp^
zN_z16(LFr17UbF=kkNuzu7jUNm+CwgL$hgY*?q<HZ_P|j3YEBqQKMHrx>`N2ZBi1M
zkpE^@gM{oKH8!@!4|dJ8r{Szm%jp74l=o|=jL8;v7|mG_Uy@Rx>@htFHTWN8VL5eo
zu?R{Yl-=j0ZtyDq$A<3Ke4ot>fjHw0A1;3a?>c;;P2N*)JEUAid$95%V&RokIM!CB
zUk}YYGHETaG?z(BUK3IxR5?Zrprt)r^Db`arV`7#AL`F)-v`N{{%bN^qf>c4K7JT7
z29R)lXRTJX*LZWhEAgsQ3p`m#GxU`W8%d<zi7_*9?WJ_i%X5X0p`3D?UE7Xe#D{hg
zpJgtz6!VD;15}3sea(k^T($qWTbGEtM?;mzdMFF|NMp@+H6{vmf_E^~Lf~neiR;AR
zPXty7Ctj>K0I>db@shTjyqqJS7;rnj!1f@QW7qwli4_ZD4*C%e*5+$nlU^cCAWtE#
zmVulB-H}=2r7A@9Mz&#XQ@(=?`izrXSGAFT*@pE!Y)tUui==0NY%r9tKaR2L8EF3s
zHtkaEYtRhk3p9<?T!ebm7dcnWu)UnFvfWP;4!v8|^D5OG<g-Gy{y129`A%}=ciXCa
zBKjQN${Q4mv3j$*DErj~jG@*Q>SNU)UB-aYZ|8}1L(O>CP$*)V=1hka{8B6XXA3F!
zH~DcKlfYVM;DQb-BLYoi+VHWqdj;vBsbxyexDp>?sl7o<kKeknQcqU}RIStzYe1@A
z{o|XlgGuv>JEtt}dtoIIlrYQ@r}<LUt6C>3MBl=tmIZ0P6G!g2X9=6++u)g30qzF9
zY=|Oi0L0`UF+*>jZmgPub3dnOPL<FZlL-|Z^P<|lMbNx+>J7Wu(Q$J#sN|ViCT2?0
zb_Le5-e7I@-{HP!GT6FGT}$wHcu9C+KSInyC5Dlj$<D~jPmJrEmDhyjwOf!hqJ~`q
z6E#+-Lr%4^o=uVH2p15$`4dotbx<XX7Ks1=00012i&x8R&M<@k0nHJBSO5UFu2sCT
O#Ao{g000001X)@qNn%p~

diff --git a/statics/physics/modern-phys-lab/run.py b/statics/physics/modern-phys-lab/run.py
index e2cb1dc..1831a83 100644
--- a/statics/physics/modern-phys-lab/run.py
+++ b/statics/physics/modern-phys-lab/run.py
@@ -5,11 +5,11 @@ import scipy
 
 colors = ['blue', 'green', 'orange', 'purple', 'red']
 wavelengths = {
-  'blue': 465,
-  'green': 520,
-  'orange': 589,
-  'purple': 390,
-  'red': 622
+  'blue': 465 * 10**(-9),
+  'green': 520 * 10**(-9),
+  'orange': 589 * 10**(-9),
+  'purple': 390 * 10**(-9),
+  'red': 622 * 10**(-9)
 }
 
 u_0 = np.empty((len(colors)))
@@ -28,7 +28,7 @@ for color_index, color in enumerate(colors):
     for index, line in enumerate(lines):
       parts = line.split('\t')
       voltage[index] = float(parts[0])
-      current[index] = float(parts[1])
+      current[index] = float(parts[1]) / 1000
 
     ## Line 1
     line_1_bound = len(lines) // 3
@@ -53,18 +53,28 @@ for color_index, color in enumerate(colors):
     ax.legend()
     ax.set_title(f"Voltage v Current ({color})")
     ax.set_xlabel("Voltage (V)")
-    ax.set_ylabel("Current (mA)")
+    ax.set_ylabel("Current (A)")
 
-    u_0[color_index] = intersect_y[0]
+    u_0[color_index] = intersect_x[0]
     wavelength_inv[color_index] = 1 / wavelengths[color]
 
 plt.show()
 
-line = polynomial.Polynomial.fit(u_0, wavelength_inv, 1)
-space = np.linspace(min(u_0), max(u_0))
-plt.plot(space, line(space), color='orange', label='Trendline 2')
+line = polynomial.Polynomial.fit(wavelength_inv, u_0, 1).convert()
+line = line.convert(domain = [min(wavelength_inv), max(wavelength_inv)]).convert()
 
-plt.scatter(u_0, wavelength_inv)
-plt.xlabel("U_0")
-plt.ylabel("1/Wavelength")
+s_2 = (1 / (len(wavelength_inv) - 2)) * sum((u_0 - line.coef[1] * wavelength_inv - line.coef[0]) ** 2)
+delta_p = len(wavelength_inv) * sum(wavelength_inv ** 2) - (sum(wavelength_inv) ** 2)
+slope_error = np.sqrt(len(wavelength_inv) / delta_p) * np.sqrt(s_2)
+print(delta_p)
+print(slope_error)
+print(s_2)
+
+space = np.linspace(min(wavelength_inv), max(wavelength_inv))
+plt.plot(space, line(space), color='orange', label=f'{line.convert()}')
+
+plt.scatter(wavelength_inv,u_0)
+plt.ylabel("U_0")
+plt.xlabel("1/Wavelength")
+plt.legend()
 plt.show()