%%timeit -n1 -r1
from numpy.linalg import eig
from scipy import random

# Losowa macierz 1000x1000 - 10^6 elementów
M=random.random((1000,1000))
# Symetryzacja
M=(M+M.T)
# Diagonalizacja
V, L=eig(M)
print '\nSuma wartości własnych:', sum(V), '\n'

%%timeit -n1 -r1
r=random.normal(size=1e7)
hist(r,bins=50,normed=True);
x=linspace(-5,5,100)
plot(x, exp(-x*x/2)/sqrt(2*pi),'r-');

%%timeit -n1 -r1
nx = 50 ; ny = 50;               ## Rozmiar siatki
epsilon = 1e-3                   ## Próg zbieżnośći
imax = 10000                     ## Maksymalna liczba iteracji

T = zeros((nx+1, ny+1)); TN = zeros((nx+1, ny+1))   ## Pole temperatury
T[[0,nx],:] = 100                                   ## Warunki brzegowe

k=1
while k<imax:
    TN[1:-1,1:-1]=(T[:-2,1:-1]+T[2:,1:-1]+T[1:-1,:-2]+T[1:-1,2:])/4.0
    errmax=np.max(TN-T)
    T[1:-1,1:-1]=TN[1:-1,1:-1]
    k+=1
    if errmax < epsilon:
        print "Zbieżność po ", k, " iteracjach."
        break

X, Y = meshgrid(linspace(0, nx, nx+1), linspace(0, ny, ny+1))
contourf(X, Y, T.T, 64, rstride=1, cstride=1) ; colorbar(); xlabel("X"); ylabel("Y"); show();

%lsmagic

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

%%bash
echo "Jestem BASH ($BASH)" 

%%perl
@numbers = (1, 2, 3, 4, 5, 6);
print $#numbers

%%latex
\begin{equation}
i\hbar\frac{\partial}{\partial t} \Psi(\mathbf{r},t) = \left [ \frac{-\hbar^2}{2m}\nabla^2 + V(\mathbf{r},t)\right ] \Psi(\mathbf{r},t)
\end{equation}