Skip to content

Instantly share code, notes, and snippets.

@AgustinParmisano
Last active February 11, 2022 23:34
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save AgustinParmisano/239404c911c99afa35eacf57a11781a4 to your computer and use it in GitHub Desktop.
Save AgustinParmisano/239404c911c99afa35eacf57a11781a4 to your computer and use it in GitHub Desktop.
lorenz cohvec con python
from ast import literal_eval
import numpy as np
from matplotlib import pyplot as plt
data_file = open("CohVec.txt", "r")
data = data_file.read()
list_data1 = data.replace("{","[")
list_data2 = list_data1.replace("}","]")
data_matrix = literal_eval(list_data2)
i=0
for data_array in data_matrix:
i+=1
#print(data_array)
data_np_array= np.array(data_array)
#print(data_np_matrix)
lorenz = data_np_array.cumsum() / data_np_array.sum()
#print(lorenz)
plt.plot(np.linspace(0.0, 1.0, lorenz.size), lorenz)
if i % 100000 == 0:
print("Va {} lineas de 133334".format(i))
plt.plot([0,1], [0,1])
plt.show()
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment