#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Fri May 12 15:23:40 2023

@author: tex
"""

import numpy as np
import matplotlib.pyplot as plt 
import statistics as st

N=100
Umin=np.random.normal(12.6,0.2,100)
Umax=np.random.normal(33.8,0.2,100)
rho=Umin/Umax
r=(1-rho)/(1+rho)
print('r=',st.mean(r))
print('erreur=',st.pstdev(r))


N=100
Umin=np.random.normal(14.5,0.2,100)
Umax=np.random.normal(36.8,0.2,100)
rho=Umin/Umax
r=(1-rho)/(1+rho)
print('r=',st.mean(r))
print('erreur=',st.pstdev(r))