Links zu Dokumentationen/Tutorials für IPython/Python/numpy/matplotlib/git sowie die Sourcodes findet ihr im GitHub Repo.

``````

In [1]:

import numpy as np
import matplotlib.pyplot as plt
%matplotlib inline

``````

## Implizit Euler

``````

In [2]:

time = impl_eul1[:,0]
V1 = impl_eul1[:,1]
V2 = impl_eul1[:,2]

``````
``````

In [3]:

plt.plot(time, V1)
plt.ylabel(r'\$V_{1}(t)\$')
plt.xlabel(r'\$t\$')
plt.grid()

``````
``````

``````
``````

In [4]:

plt.plot(time[:100], V1[:100])
plt.ylabel(r'\$V_{1}(t)\$')
plt.xlabel(r'\$t\$')
plt.grid()

``````
``````

``````
``````

In [5]:

time = impl_eul_e1[:,0]
V1 = impl_eul_e1[:,1]
V2 = impl_eul_e1[:,2]

``````
``````

In [6]:

plt.plot(time, V1)
plt.ylabel(r'\$V_{1}(t)\$')
plt.xlabel(r'\$t\$')
plt.grid()

``````
``````

``````
``````

In [7]:

plt.plot(time[:1000], V1[:1000])
plt.ylabel(r'\$V_{1}(t)\$')
plt.xlabel(r'\$t\$')
plt.grid()

``````
``````

``````
``````

In [8]:

time = impl_eul_e2[:,0]
V1 = impl_eul_e2[:,1]
V2 = impl_eul_e2[:,2]

``````
``````

In [9]:

plt.plot(time, V1)
plt.ylabel(r'\$V_{1}(t)\$')
plt.xlabel(r'\$t\$')
plt.grid()

``````
``````

``````
``````

In [10]:

plt.plot(time[:10000], V1[:10000])
plt.ylabel(r'\$V_{1}(t)\$')
plt.xlabel(r'\$t\$')
plt.grid()

``````
``````

``````
``````

In [11]:

time = impl_eul_e3[:,0]
V1 = impl_eul_e3[:,1]
V2 = impl_eul_e3[:,2]

``````
``````

In [12]:

plt.plot(time, V1)
plt.ylabel(r'\$V_{1}(t)\$')
plt.xlabel(r'\$t\$')
plt.grid()

``````
``````

``````
``````

In [13]:

plt.plot(time[:100000], V1[:100000])
plt.ylabel(r'\$V_{1}(t)\$')
plt.xlabel(r'\$t\$')
plt.grid()

``````
``````

``````

## Implizit MP

``````

In [14]:

time = impl_mp_1[:,0]
V1 = impl_mp_1[:,1]

``````
``````

In [15]:

plt.plot(time, V1)
plt.ylabel(r'\$V_{1}(t)\$')
plt.xlabel(r'\$t\$')
plt.grid()

``````
``````

``````
``````

In [16]:

time = impl_mp_e1[:,0]
V1 = impl_mp_e1[:,1]

``````
``````

In [17]:

plt.plot(time, V1)
plt.ylabel(r'\$V_{1}(t)\$')
plt.xlabel(r'\$t\$')
plt.grid()

``````
``````

``````
``````

In [18]:

time = impl_mp_e2[:,0]
V1 = impl_mp_e2[:,1]
V2 = impl_mp_e2[:,2]

``````
``````

In [19]:

plt.plot(time, V1)
plt.ylabel(r'\$V_{1}(t)\$')
plt.xlabel(r'\$t\$')
plt.grid()

``````
``````

``````
``````

In [20]:

time = impl_mp_e3[:,0]
V1 = impl_mp_e3[:,1]
V2 = impl_mp_e3[:,2]

``````
``````

In [21]:

plt.plot(time, V1)
plt.ylabel(r'\$V_{1}(t)\$')
plt.xlabel(r'\$t\$')
plt.grid()

``````
``````

``````
``````

In [22]:

plt.plot(time[:100000], V1[:100000])
plt.ylabel(r'\$V_{1}(t)\$')
plt.xlabel(r'\$t\$')
plt.grid()

``````
``````

``````
``````

In [ ]:

``````