load(interpol) $
load("opsubst") $
cspline([[0,0],[1,1],[2,4]]) $
f(x):=''% $
zeroize(x,y,z) := 0

opsubst(zeroize, diff, diff(f(x),x)), eval, eval;