ส่งงาน week10

week10/RichExtra.py

+from math import e
+def Extrapolation(h1,p):
+    g1=(e**(-(1+h1))-2*e**(-1)+e**(-(1-h1)))/(h1**2)
+    h2=h1/2
+    g2=(e**(-(1+h2))-2*e**(-1)+e**(-(1-h2)))/(h2**2)
+    return (((2**p)*g2)-g1)/(2**p-1)
+Extrapolation(0.64,2)
+
+def G(h):
+    cda = {0.64:0.380610,
+           0.32:0.371035,
+           0.16:0.368711,
+           0.08:0.368281,
+            0.04:0.36875,
+            0.02:0.37,       
+            0.01:0.38,      
+            0.005:0.40,       
+            0.0025:0.48,       
+            0.00125:1.28}
+    return cda[h]
+def RichExtra(h,p):
+    return (((2**p)*G(h/2))-G(h))/(2**p-1)
+RichExtra(0.64,2)
+
+####Output 0.36784333333333336