import numpy as np
import random as rd
import matplotlib.pyplot as plt

##

def NS(a):
    return('{:.2E}'.format(a))

########################################"
## TP 30 câble coaxial

def c(Z):
    return 1/Z/1e-10


def li(Z):
    return Z**2*1e-10

# On connait les valeurs et écart-type :

Zc=46
u_Zc=2


## tirages aléatoires


tirages=[]
for i in range(100000):
    tirages.append(c(rd.normalvariate(Zc,u_Zc)))

C,u_C=np.average(tirages),np.std(tirages)
print('c = ',NS(C),'u_c = ',NS(u_C))

##

tirages=[]
for i in range(100000):
    tirages.append(li(rd.normalvariate(Zc,u_Zc)))

l,u_l=np.average(tirages),np.std(tirages)
print('l = ',NS(l),'u_l = ',NS(u_l))


def F(x,y):
    return x*y/2

tau=1.04e-6
u_tau=1e-8

tirages=[]
for i in range(100000):
    tirages.append(F(rd.normalvariate(C,u_C),rd.normalvariate(tau,u_tau)))