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]:

    
impl_eul1 = np.loadtxt('data/ex14_impl_eul_1.out')
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]:

    
impl_eul_e1 = np.loadtxt('data/ex14_impl_eul_e-1.out')
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]:

    
impl_eul_e2 = np.loadtxt('data/ex14_impl_eul_e-2.out')
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]:

    
impl_eul_e3 = np.loadtxt('data/ex14_impl_eul_e-3.out')
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]:

    
impl_mp_1 = np.loadtxt('data/ex14_impl_mp_1.out')
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]:

    
impl_mp_e1 = np.loadtxt('data/ex14_impl_mp_e-1.out')
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]:

    
impl_mp_e2 = np.loadtxt('data/ex14_impl_mp_e-2.out')
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]:

    
impl_mp_e3 = np.loadtxt('data/ex14_impl_mp_e-3.out')
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 [ ]: