kop-3

import matplotlib.pyplot as plt
import numpy as np
m=1200
n=1200
x0=0.3
L=np.full(m*n,0.00000)
Y=np.full(m*n,0.00000)
l=-0.1326
I =[i for i in range(n)]
J =[j for j in range(m)]
for j in J:
    xi=x0
    for i in I:
        xi= -(31.3344 * l) + (11.24 * l * xi) - (l * xi * xi)
        L[m*(j-1)+i-1]=l
        Y[m*(j-1)+i-1]=xi
    x0=xi
    print(l)
    l += 0.00013
print(Y)
print(L)
plt.plot(L,Y,ls=" ", marker=",")
plt.show()