A numerical implementation of the unified Fokas transform for evolution problems on a finite interval