pwd

ls *

files = !ls *

files

import bisect

bisect?

bisect??

bisect.bisect_right?

import math
def f(a, b):
    return a*math.cos(2*math.pi/b)

f(1, 2)

f(1, 0)

import time
for i in range(10):
    print i,
    time.sleep(0.5)

from IPython.core.display import Image
Image(filename='mirabellug.jpeg')  # fichier local

Image('http://ipython.org/_static/IPy_header.png')

from IPython.display import HTML

s = """<table>
<tr>
<th>Header 1</th>
<th>Header 2</th>
</tr>
<tr>
<td>row 1, cell 1</td>
<td>row 1, cell 2</td>
</tr>
<tr>
<td>row 2, cell 1</td>
<td>row 2, cell 2</td>
</tr>
</table>"""

h = HTML(s)
h

HTML("<iframe src=http://www.mirabellug.org width=1000 height=500>")

from IPython.lib.display import YouTubeVideo
YouTubeVideo("t4gjl-uwUHc")

%pylab inline

f, (ax1, ax2) = subplots(ncols=2)    # renvoie une figure et deux axes associés sur deux colonnes
x = linspace(0, 2*pi, 50)       # tableau de 50 points uniformément répartis entre 0 et 2pi
line1, = ax1.plot(x, sin(x), 'b-')   # courbe de la fonction sinus, en trait continu bleu
line2, = ax1.plot(x, cos(x), '-', color='#333333')   # fonction cosinus, en trait continu et couleur définie par un code rgb
line3, = ax2.plot(x, x**2, 'g-')   # parabole, en vert 
line4, = ax2.plot(x, 9*x - 4*x**2 + 0.5*x**3, 'k-')    # polynôme plus complexe, en noir

ax1.set_title("Sin et cos", fontsize=20)
ax2.set_title(u"Polynômes", fontsize=20)
ax1.set_xlabel("x", fontsize=16)
ax2.set_xlabel("x", fontsize=16)
ax1.set_ylabel("y", fontsize=16)
ax1.grid(True)
ax2.grid(True)
line2.set_color('#00adef')
line3.set_marker('*')

f

Devol_sans = genfromtxt("evolution_diametre_sans_eclaircies.csv", skiprows=1, names=True, delimiter=";") 

Devol_sans['Dmoy']

plot(Devol_sans['temps'], Devol_sans['Ddom'], 'k--')

plot(Devol_sans['temps'], Devol_sans['Ddom'], 'k--', label=u'Diamètre dominant (cm)')
plot(Devol_sans['temps'], Devol_sans['Dmoy'], 'b--', label=u'Diamètre moyen (cm)')
legend(loc='upper left')
grid(True)
title(u"Évolution des diamètres dominant et moyen\nen l'absence d'éclaircie", fontsize=14)
xlabel(u"Âge du peuplement (années)", fontsize=12)

Devol_avec = genfromtxt("evolution_diametre_avec_eclaircies.csv", skiprows=1, names=True, delimiter=";") 

plot(Devol_avec['temps'], Devol_avec['Ddom'], 'k--', label=u'Diamètre dominant (cm)')
plot(Devol_avec['temps'], Devol_avec['Dmoy'], 'b--', label=u'Diamètre moyen (cm)')
grid(True)
legend(loc='upper left')
title(u"Évolution des diamètres dominant et moyen\navec des éclaircies", fontsize=14)
xlabel(u"Âge du peuplement (années)", fontsize=12)

f, (ax1, ax2) = subplots(ncols=2, figsize=(10,6))
ax1.plot(Devol_sans['temps'], Devol_sans['Ddom'], 'k--', label=u'Diamètre dominant (cm)')
ax1.plot(Devol_sans['temps'], Devol_sans['Dmoy'], 'b--', label=u'Diamètre moyen (cm)')
ax1.grid(True)
ax1.legend(loc='upper left')   # légende en haut à gauche de l'axe
ax1.set_title(u"Sans éclaircies")
ax2.plot(Devol_avec['temps'], Devol_avec['Ddom'], 'k--', label=u'Diamètre dominant (cm)')
ax2.plot(Devol_avec['temps'], Devol_avec['Dmoy'], 'b--', label=u'Diamètre moyen (cm)')
ax2.grid(True)
ax2.legend(loc='upper left')
ax2.set_title(u"Avec éclaircies")

Ddistrib_sans = genfromtxt("classes_diametre_sans_eclaircies.csv", skiprows=1, names=True, dtype=['S5', 'i'], delimiter=";")

nb_classes = Ddistrib_sans.size
fig, ax = subplots(figsize=(7,4))
ax.set_xlabel(u"Diamètre (cm)", fontsize=14)
ax.set_ylabel("Nombre d'arbres", fontsize=14)
ax.set_title(u"Répartition par classes de diamètres\nen l'absence d'éclaircie", fontsize=16)
width=0.5
ax.set_xticks(arange(nb_classes) + width/2)
ax.set_xticklabels(Ddistrib_sans['diametre'], rotation=70, fontsize=11)
barres = ax.bar(arange(nb_classes), Ddistrib_sans['N'], width=width)

Ddistrib_avec = genfromtxt("classes_diametre_avec_eclaircies.csv", skiprows=1, names=True, dtype=['S5', 'i'], delimiter=";")

nb_classes = Ddistrib_avec.size
fig, ax = subplots()
ax.set_xlabel(u"Diamètre (cm)", fontsize=14)
ax.set_ylabel("Nombre d'arbres", fontsize=14)
ax.set_title(u"Répartition par classes de diamètres\navec des éclaircies", fontsize=16)
width=0.5
ax.set_xticks(arange(nb_classes) + width/2)
ax.set_xticklabels(Ddistrib_avec['diametre'], rotation=70, fontsize=11)
barres = ax.bar(arange(nb_classes), Ddistrib_avec['N'], width=width)

%loadpy http://matplotlib.org/mpl_examples/pylab_examples/animation_demo.py

%%script perl
$_="krJhruaesrltre c a cnp,ohet";$_.=$1,print$2while s/(..)(.)//;

%%script ruby
puts "Hello from Ruby #{RUBY_VERSION}"

import this