Spline1

{

This example shows a 1D spline table, and compares derivatives with

analytic values.

See "Splinetable.pde" for a 2D table construction.

}

title '1D Spline table'

select

   regrid=off

definitions

    u = splinetable("gauss1.tbl")

    gu = exp(-16*x^2)

boundaries

    Region 1

       start(-1,-1)

       line to (1,-1) to (1,1) to (-1,1) to close

plots

           contour(u)

           contour(dx(u))

           contour(dy(u))

       elevation(u, gu) from(-1,0) to (1,0)

       elevation(dx(u), dx(gu)) from(-1,0) to (1,0)

       elevation(dxx(u), dxx(gu)) from(-1,0) to (1,0)

end