My case is 3 level 3 factor experiment, and I would like to fit the curved surface with non-linear fitting. The reason why I use 3 level design is that I know the curve between them should be the exponential, quadratic, or even higher relation. It seems to me that Response Surface Design(CCD or Box-Behnken) only support 2 level factors. Do you have any idea how to do that? Thanks.