In this book the solution representations obtained from a recently developed Fokas integral method for solving boundary value problems for linear evolution PDEs are evaluate numerically. In particular, the case of the linear KdV equation is considered. The Fokas method is quite general and it is therefore of wider interest to assess its competitiveness for numerical purposes. Until now pseudospectral methods have been know to be the most accurate numerical scheme for smooth functions. To compare the two methods, the linear KdV equation is computed numerically using both, a pseudospectral method and the direct evaluation of the integral representation. The two methods are compared for accuracy and speed of the numerical computation, showing that for linear evolutionary PDEs the numerical implementation of Fokas method is much faster and more accurate than a pseudospectral method. The nonlinear KdV equation is also looked at using pseudospectral methods and a motivation for a possible hybrid method which would use both the Fokas and pseudospectral methods together is given.