     plot2d(['(f(x)),[discrete,p]],[x,0,10],
          [gnuplot_curve_styles,
               [&quot;with lines&quot;,&quot;with points pointsize 3&quot;]])$
(%i7) /* New call, but giving values at the derivatives */
      cspline(p,d1=0,dn=0);
              3          2
        1949 x    11437 x    17027 x   1247
(%o7)  (------- - -------- + ------- + ----) charfun2(x, minf, 3)
         2256       2256      2256     752
            3          2
      1547 x    35581 x    68068 x   173546
 + (- ------- + ---