Week 9: implemented cosine taper for wave function, need to work on direct method without lumping, wave propagation velocity is still not correct