Added solution from Kassam and Trefethen, ETD