import numpy as np
from __future__ import division

def ErrorRule4(p1,A,dA,p2,B,dB,p3,C,dC,LF):
    dLF = LF * (np.sqrt(((p1*dA/A)**2)+((p2 *dB/B)**2)+((p3*dC/C)**2)))
    return dLF

def ErrorRule3(dA,dB):
    dC = np.sqrt(dA**2+dB**2)
    return dC

def f(lam,g,HM,Mu):
    freq = ((1/lam)(np.sqrt((HM*g)/Mu)))
    return freq

lam1 = .690
lam2 = .510
lam3 = .706
lam4 = .526

Unlam1 = .001
Unlam2 = .001
Unlam3 = .001
Unlam4 = .001

McordY = .001430
McordB = .000440
UnMcordY = .000002
UnMcordB = .000002

LcordY = 1.686
UnLcordY = .001
LcordB = 1.365
UnLcordB = .001

Mhang1 = .616
Mhang2 = .340
Mhang3 = .200
Mhang4 = .120
UnMhang1 = .002
UnMhang2 = .002
UnMhang3 = .002
UnMhang4 = .002

MuCordY = .0008
MuCordB = .0003


f1 = 125.895
f2 = 126.542
f3 = 114.588
f4 = 119.03
fAVG = 121.23


z = ErrorRule3(.000002,.00001)
print(z)
Uf1 = ErrorRule4(1,Unlam1,lam1,1,UnMhang1,Mhang1,1,.00001,.0008,f1)
print(Uf1)