Include quadratic terms in experimental design

I'm designing a discrete choice experiment with four attributes: three topic-related attributes and one cost attribute. Each attribute has four levels. I want to estimate a quadratic utility function of the form:
V_i = beta_1 [Attribute1]_i + beta_2 [Attribute1]_i ^2 + beta_3 [Attribute2]_i + beta_4 [Attribute2 %]_i ^2 + beta_5 [Attribute3]_i + beta_6 [Attribute3]_i ^2 + beta_7[Cost ~ lognormal]
I saw in the Ngene manual (Section 5.6) that only simple linear functions are supported via fcn(...). Specifically, it states that "Currently only simple linear functions are supported in Ngene, therefore only plus ( + ), minus ( - ), attribute names and constants are allowed."
Does this mean I cannot include squared terms in the utility function (I tried it in Ngene, and it failed)? If so, are there any workarounds or supported methods for implementing quadratic terms in Ngene?
Alternatively, I can generate an efficient design using only the linear terms, under the assumption that the variation in the linear attribute will still carry the relevant information needed to estimate a quadratic term during model estimation, given the deterministic relationship between an attribute and its square. That is, I’d estimate the utility function with both linear and squared terms even if the design only explicitly includes the linear levels. However, if the quadratic terms are not explicitly included in the design, it isn't included in the Fisher information matrix and will affect the efficiency of the design.
Any insights or advice would be greatly appreciated!
V_i = beta_1 [Attribute1]_i + beta_2 [Attribute1]_i ^2 + beta_3 [Attribute2]_i + beta_4 [Attribute2 %]_i ^2 + beta_5 [Attribute3]_i + beta_6 [Attribute3]_i ^2 + beta_7[Cost ~ lognormal]
I saw in the Ngene manual (Section 5.6) that only simple linear functions are supported via fcn(...). Specifically, it states that "Currently only simple linear functions are supported in Ngene, therefore only plus ( + ), minus ( - ), attribute names and constants are allowed."
Does this mean I cannot include squared terms in the utility function (I tried it in Ngene, and it failed)? If so, are there any workarounds or supported methods for implementing quadratic terms in Ngene?
Alternatively, I can generate an efficient design using only the linear terms, under the assumption that the variation in the linear attribute will still carry the relevant information needed to estimate a quadratic term during model estimation, given the deterministic relationship between an attribute and its square. That is, I’d estimate the utility function with both linear and squared terms even if the design only explicitly includes the linear levels. However, if the quadratic terms are not explicitly included in the design, it isn't included in the Fisher information matrix and will affect the efficiency of the design.
Any insights or advice would be greatly appreciated!