# -*- coding: utf-8 -*-
"""
Created on Wed Sep  5 17:18:57 2018

@author: Benjamin
"""


import numpy as np               # Pour la fonction linspace, zeros et la trigo
import matplotlib.pyplot as plt  # Boîte à outils graphiques

# Definition de la fonction a tracer
def C(x):
    return(2*np.sqrt(x) / (1+x))

# Liste des abscisses
L_x = [k/10 for k in range(100)] # echelle de 0 a 9,9 par 0,1
# Liste des ordonnees
L_C = [C(x) for x in L_x]

plt.figure() # Creer nouvelle figure
plt.subplot(1,2,1) # 1er sous-graphe
plt.plot(L_x,L_C)  # Tracer.
plt.xlabel("x=I1/I2")   # Legende abscisses
plt.ylabel("Contraste") # Legende ordonnees
plt.title('Evolution du contraste  : 2*np.sqrt(x) / (1+x)') # Titre
plt.grid() # Grille sur le graphe



# Liste des abscisses
L_x_log = [10**(k/10) for k in range(-100,101)] # de 10**-10 a 10**10
L_C_log = [C(x) for x in L_x_log]

plt.subplot(1,2,2) # 2e sous-graphe
plt.plot(L_x_log,L_C_log)
plt.xlabel("x=I1/I2")
plt.ylabel("Contraste")
plt.title('Evolution du contraste en échelle log')
plt.yscale('log')
plt.xscale('log')
plt.grid()
plt.show() # Afficher le graphe