From 747903f4e363cdc34076527ddd19d1563cf82cb6 Mon Sep 17 00:00:00 2001
From: Andrej <andrej@baza>
Date: Wed, 9 Jun 2021 11:09:26 +0200
Subject: [PATCH] zacel pripravljati primere za seminar

---
 __pycache__/boris.cpython-37.pyc | Bin 2886 -> 3462 bytes
 boris.py                         |  14 ++++++++++++++
 report.py                        |  23 +++++++++++++++++++++++
 3 files changed, 37 insertions(+)
 create mode 100644 report.py

diff --git a/__pycache__/boris.cpython-37.pyc b/__pycache__/boris.cpython-37.pyc
index 434fcdb130acf84221517ccc00c9c87e7d2e435f..e9bc86a8a96377e40b05a2c7420dc70d56d1a773 100644
GIT binary patch
delta 900
zcmZWn&rcIU6rP#gZMVBn+EPj@s1y(=7-J;yhBc@O9^~RdSw$?oVL=Ne)5;IorV-Ga
zho&b@OxnM|7*F1e2cz-mZBHJ&aMi03-`muv@lEEN`M&vP=e^l^pL)_4o$l<k7)ERJ
z?cKUvkEZ$8`sfHhRiB-xjX)Q~4!C1O91;i2>40uX;@AnJ&<kl?i9!!Bh`ItY%NiVm
zKGeEo4Epa2^#zW@33SKdWH2TqZ+61~WFhNnE$$=(o`u0S?|~sFC2g6MJu)s6GL?qm
z{m^aU^g3x7mA(6%GvL=*ZtT^MPdaxRAB^*Sym4gg@rl=aq0<7Tmu8_hgVp$AduK_f
zx-d-aHug){%ij<xtjSxfs(X4<;AnV8Q}Y-W!X67x0MD>xkzs!+$VdULrE!LCy~S{3
zp#S4G0#*yGX;uxj4eiay$k1lE9`=lT#(y87PyXH!^1={;FqjYl4Xme;;P4RG&9J9Y
znjUY77G4XcWVTgN`y*qW5y})8-de}3;#j+)+fX@#!y}Btp9t>rXB2Pp?a9ZyI~m01
z{2@J#oZ(w5vH%NPt58|-?O<M%>n;+rv?<j}Vck_{(dipYwNkO99!7?DtHT4gcB}I_
zOP!-5Cu!L%KUmx-mMUd6MO1gyD)>f0snV9ZK(0K&MS^LBoT)N2r$48zBKV?IR+oug
z2^ehq0(MlElmG_VIAx*ekZY)xQwWTk*v?=dz&6G^cnX`k(TGPL+4H#molI)3@gkzE
zIh@`h>NGOe49_$ct+6Spz_->475Z**bx{ctxODWjg$+e0i#}gcL$nbyU}}PZo+2Le
Z&bo;Ywa@O4Cz4&1iYa0|CUlX~e*qOC*ct!;

delta 320
zcmZpZJ|@QJ#LLUY00a{j?M^tuIgw9-@xnxPHAbb0e%kfQDf}q{EsRmzDS|0NEeugS
zDZ(itK$<s2B}F_%0w}|m%AX>c%~W(IMJh!aC?k*}1EvL285anp$fn4p$fYwcWQ-C9
z^W{?Hp?r}Pg(%Tf?o{DakyL?H!Bo){#b&k`rYNx}@l?K4v1Ud_h7?AYU<OUq%>s<4
zSr~6kzQbz9s6APjt%*@*^E5UcCZNz^jzmVU$zI%z-2Olhd<OZV$Y=63?jkN9ATt2O
z^_?8V<H0B~c{Y!Rg)vaz7FTX!NkLA2Nls=`5fey&G4&QpZhS#e6p&2Kxy6`Vqzuxi
i4I*?vgx4hA2u6*`k-Ree!W@i1$Rxlhz{J8N#sUBnc}@}l

diff --git a/boris.py b/boris.py
index 028f6b9..a478979 100644
--- a/boris.py
+++ b/boris.py
@@ -74,6 +74,20 @@ def plot3(X,enote='m',zemlja=False):
 
     plt.show()
 
+def rplot(X, enote='m', xy=[0,1],os='xy'):
+    fig = plt.figure()
+    ax = fig.subplots()
+    ax.plot(X[:,xy[0]],X[:,xy[1]], 'k')
+    ax.grid()
+    dx = X[-1,xy[0]]-X[-2,xy[0]]
+    dy = X[-1,xy[1]]-X[-2,xy[1]]
+    ax.arrow(X[-1,xy[0]],X[-1,xy[1]],dx,dy,head_width=0.1,head_length=0.2,
+            overhang=0.3, zorder=10,color='k')
+    ax.set_xlabel(r'${}$ [{}]'.format(os[0],enote))
+    ax.set_ylabel(r'${}$ [{}]'.format(os[1],enote))
+    plt.show()
+
+
 
 if __name__ == "__main__": # testni del kode
     E = lambda x: np.array([0.0,0.2,0.3])
diff --git a/report.py b/report.py
new file mode 100644
index 0000000..d4f97f1
--- /dev/null
+++ b/report.py
@@ -0,0 +1,23 @@
+#!/usr/bin/python3 
+
+import numpy as np
+import matplotlib.pyplot as plt
+from matplotlib import rc
+rc('font',**{'family':'serif','serif':['Computer Modern']})
+rc('text', usetex=True)
+
+from boris import *
+
+B = lambda x: np.array([0.,0.,x[0]**2+x[1]**2])
+E = lambda x: np.array([0.,0.,0.])
+
+x0 = np.array([0.,-1.,0.])
+v0 = np.array([0.,1.,0.])
+
+B = lambda x: np.array([0.,0.,1.])
+E = lambda x: np.array([0.,-1.,0.])
+E = lambda x: np.array([np.cos(x[0]),0.3,0.])
+X,V = boris(x0,v0,E,B,1e-3, 100)
+
+plot3(X)
+rplot(X)