Add file phi_no_interpolation.m that computes phi from scratch each time. This is slower than interpolation but more accurate.