function y = cubic ( x )
  y = 816*x.^3 - 3835*x.^2 + 6000*x - 3125;
  return
end
