Function f

code/griddata/spiral_interp.py:32–34

(x,y)

import numpy as np
from scipy.interpolate import griddata
def f(x,y):
   return x + 2*(np.cos(7*x)+
      np.sin(5*y)) + 3*y*np.sin(6*x)
K = 0.005
# subsample f(x,y) at 114 points
t = np.linspace(0,800,num=1100)