#ms.ple("torque", M, dM)
#ms.pve("deflection force", Dt, dDt)
#ms.ps("linreg vs graph", 0.0217, 0.0004, Dt, dDt)
# 1.2 Deflecting force determination by period time measurement
Jd = md / 2.0 * rd**2
dJd = rd * ms.sqrt((0.5 * rd * dmd)**2 + (md * drd)**2)
Dp = 4.0 * ms.pi**2 * Jd / (T2**2 - T1**2)
dDp = 4.0 * ms.pi**2 / (T2**2 - T1**2) * ms.sqrt(dJd**2 +
((2.0 * Jd * T1 * dT1)**2 +
(2.0 * Jd * T2 * dT2)**2) /
(T2**2 - T1**2)**2)
print()
ms.pve("T1", T1, dT1, False)
ms.pve("T2", T2, dT2)
ms.pve("Jd", Jd, dJd)
ms.pve("Dp", Dp, dDp)
# Deviation between 1.1 and 1.2
ms.ps("Deviation", Dt, Dp, dDt, dDp, False)
D = ms.mean_value([Dt, Dp])
dD = 1 / 2 * ms.sqrt(dDt**2 + dDp**2)
ms.pve("Mean Value", D, dD)
# 2 Verification of the theorem of Steiner
a2 = []
da2 = []
J = []
dJ = []
return deg + mins / 60.0
lamHg = [690.7, 623.4, 579.1, 577.0, 546.1, 499.2, 491.6, 435.8, 407.8, 404.7]
phiHg = [mtd(181, 56), mtd(180, 34), mtd(180, 18), mtd(180, 16), mtd(179, 59), mtd(179, 30), mtd(179, 26), mtd(178, 31), mtd(177, 55), mtd(177, 52)]
dphiHg = [mtd(0, 5)] * len(phiHg)
phiHe = [mtd(180, 49), mtd(180, 22), mtd(179, 34), mtd(179, 28), mtd(179, 9), mtd(178, 48)]
dphiHe = [mtd(0, 5)] * len(phiHe)
phiH = [mtd(180, 45), mtd(179, 21), mtd(178, 29), mtd(178, 0)]
dphiH = [mtd(0, 5)] * len(phiH)
for i in range(len(phiHg)):
ms.pv("lambda[" + str(i) + "]", lamHg[i])
print()
for i in range(len(phiHg)):
ms.pve("deltaHg(lambda)[" + str(i) + "]", phiHg[i], dphiHg[i])
print()
for i in range(len(phiHe)):
ms.pve("deltaHe(lambda)[" + str(i) + "]", phiHe[i], dphiHe[i])
print()
for i in range(len(phiH)):
ms.pve("deltaH(lambda)[" + str(i) + "]", phiH[i], dphiH[i])
for i in range(len(phiHg)):
phiHg[i] *= m.pi / 180.0
for i in range(len(phiHe)):
phiHe[i] *= m.pi / 180.0
for i in range(len(phiH)):
phiH[i] *= m.pi / 180.0
import math as m
import measure as ms
# Gasthermometer
p0 = 1013.25
pl = 1020.2
dpl = 0.2
p = [916, 945, 991, 1027, 1061, 1092, 1125, 1163, 198, 1226, 1245]
dp = [1, 5, 6, 5, 5, 5, 3, 2, 5, 4, 2]
pc = 653
dpc = 2
pn = 255
dpn = 1
pnb = p[len(p) - 1] * p0 / pl
dpnb = p0 / pl * m.sqrt(dp[len(p) - 1]**2 + p[len(p) - 1]**2 / pl**2 * dpl**2)
ms.pve("pnb", pnb, dpnb)
mw = 497.61e-3
dmw = 0.13e-3
mw -= mk
dmw = m.sqrt(dmw**2 + dmk**2)
Tw = 55.8
dTw = 0.2
Tl = 24.1
dTl = 0.1
Tm = 53.25 # 53.19
dTm = 0.27 # 0.4
W = mw * cw * (Tw - Tm) / (Tm - Tl)
dW = m.sqrt((Tw - Tm)**2 * ((cw * dmw)**2 + (mw * dcw)**2) + (mw * cw)**2 *
(dTw**2 + (((Tl - Tw) * dTm)**2 + ((Tw - Tm) * dTl)**2) /
(Tm - Tl)**2)) / (Tm - Tl)
ms.pve("mw", mw, dmw)
ms.pve("Wasserwert", W, dW)
# T=100°C
R = 8.3144598
Mp = [207.2e-3, 26.9815385e-3, 12.011e-3]
p = 1024.6
dp = 2.0
mp = [501.773e-3, 131.08e-3, 125.05e-3]
dmp = [0.06e-3, 0.04e-3, 0.02e-3]
mw = [612.38e-3, 600.03e-3, 615.13e-3]
dmw = [0.04e-3, 0.09e-3, 0.05e-3]
Tw = [27.4, 24.7, 23.4]
dTw = [0.2, 0.2, 0.2]
Tp = [100 + 0.0276 * (p - p0) for i in range(len(mp))]
#print(f)
#print(d_f)
#print(sigma)
# Auflösungsvermögen Mikroskop
s = 41.5
ds = 1.0
b = 250
f = 40
lam = 550e-6
beta = b / f - 1.0
a = [0.2, 0.5, 0.4]
da = ms.std_dev_m(a)
a = ms.mean_value(a)
alpha = m.atan(a / (2 * s))
dalpha = 1.0 / m.sqrt(4 * s**2 + a**2) * m.sqrt(da**2 + (a / s)**2 * ds**2)
g = 1 / (2 * beta)
gMin = 0.61 * lam / m.sin(alpha)
dgMin = abs(0.61 * lam * (m.cos(alpha) / (m.sin(alpha)**2)) * dalpha)
ms.pve("a", a, da)
ms.pv("beta", beta)
ms.pve("alpha", alpha * 180 / m.pi, dalpha * 180 / m.pi)
ms.pv("g", g)
ms.pve("gMin", gMin, dgMin)
ms.pv("sigma", abs(gMin - g) / dgMin)