# Centrale 2 MP : Panneaux solaires anti reflets


import matplotlib.pyplot as plt
import numpy as np

## Valeurs numériques

N = 1.5
n = np.sqrt(N) # à justifier
c = 3.0E8
omega0 = 3.14E15
deltaomega = 2.35E15

## À adapter ou à justifier

R=((N-1)/(N+1))**2

def delta(e):
    return(2*n*e)
    
def RAR(e):
    return(0.5*R*(1+(np.sin(deltaomega*delta(e)/(2*c)))*(np.cos(omega0*delta(e)/c))/(deltaomega*delta(e)/(2*c))))
    
plt.figure()
e = np.linspace(0,500*10**(-9),100)
plt.plot(e,RAR(e))
plt.xlabel("Epaisseur")
plt.ylabel("Coefficient de réflexion")
plt.grid()
plt.show()