In [32]:

    

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


    

In [33]:

    

expl_omeg1 = np.loadtxt('ex08_expl_omeg1_1e-3.out')
time = expl_omeg1[:,0]
V1 = expl_omeg1[:,1]
V2 = expl_omeg1[:,2]


    

In [34]:

    

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


    

In [35]:

    

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


    

In [14]:

    

impl_omeg1 = np.loadtxt('ex08_impl_omeg1_1e-3.out')
time = expl_omeg1[:,0]
V1 = expl_omeg1[:,1]
V2 = expl_omeg1[:,2]


    

In [36]:

    

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


    

In [38]:

    

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


    

In [39]:

    

expl_omeg1000 = np.loadtxt('ex08_expl_omeg1000_1e-4.out')
time = expl_omeg1[:,0]
V1 = expl_omeg1[:,1]
V2 = expl_omeg1[:,2]


    

In [41]:

    

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


    

In [42]:

    

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


    

In [43]:

    

expl_omeg1000 = np.loadtxt('ex08_impl_omeg1000_1e-3.out')
time = expl_omeg1[:,0]
V1 = expl_omeg1[:,1]
V2 = expl_omeg1[:,2]


    

In [30]:

    

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


    

In [31]:

    

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