Problem implementing KL expansion with square-exponential kernel: output just looks like Gaussian
I've had issue with this -- please see the errata of the reference. The eigenfunctions are not scaled in the original version of the book. This errata should solve it.