This is an old revision of the document!
25. 11. Zdrojaky
12. 11.
5. 11.
#!/usr/bin/env python from pylab import * N = 500 xmax = 5 x = linspace(-xmax,xmax,N) dx = x[1]-x[0] T = -1/2.*(eye(N,k=-1) - 2*eye(N) + eye(N,k=1))/dx/dx V = diag(x*x/2.) plot(x, diag(V), "k-") H = T + V E, psi = linalg.eigh(H) o = 0 twinx() for p,e in zip(psi.T, E)[:10]: prob = abs(p)**2 plot(x, o+prob, label="E={:.4}".format(e)) o += 1.1*amax(prob) legend() savefig("osc.png") show()
29. 10.
data = rand(2,10) # data.shape == (2,10) # takto NE ! #x,y = [],[] #for i in range(10): # x.append(data[0,i]) # y.append(data[1,i]) # takto ANO :) x,y = data
22. 10.