from mpl_toolkits import mplot3d
import numpy as np
import matplotlib.pyplot as plt
from pathlib import Path
import os

os.chdir(Path(__file__).parent)

func=np.loadtxt("lorentz_prooton.dat")
x=func[:,0]
y=func[:,1]
z=func[:,2]

fig=plt.figure(figsize=(10,10))
p=plt.axes(projection="3d")
#p.set_aspect("equal")
p.set_xlim(-40000,40000)
p.set_ylim(-40000,40000)
p.set_zlim(-40000,40000)
p.plot(x,y,z)
plt.grid()
plt.show()