In [1]:

    

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


    

In [2]:

    

adaptive = np.loadtxt('data/ex18.out')
time = adaptive[:,0]
h = adaptive[:,1]
C1 = adaptive[:,2]
C2 = adaptive[:,3]
C3 = adaptive[:,4]
C4 = adaptive[:,5]
C5 = adaptive[:,6]


    

In [3]:

    

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


    

In [4]:

    

plt.plot(time, C2)
plt.ylabel(r'$C_2(t)$')
plt.xlabel(r'$t$')
plt.grid()


    

In [5]:

    

plt.plot(time, C3)
plt.ylabel(r'$C_3(t)$')
plt.yscale('log')
plt.xlabel(r'$t$')
plt.grid()


    

In [6]:

    

plt.plot(time, C4)
plt.ylabel(r'$C_4(t)$')
plt.xlabel(r'$t$')
plt.grid()


    

In [7]:

    

plt.plot(time, C5)
plt.ylabel(r'$C_5(t)$')
plt.yscale('log')
plt.xlabel(r'$t$')
plt.grid()


    

In [ ]: