y=a0*sqrt((X1^((1+a1)/a1))/(X2^(1/a2)+a3^X3))
sqrt(X1^((a1+1)/a1)/(a3^X3+X2^(1/a2)))
(a0*(1/a1-(a1+1)/a1^2)*X1^((a1+1)/a1)*log(X1))/(2*sqrt(X1^((a1+1)/a1)/(a3^X3+X2^(1/a2)))*(a3^X3+X2^(1/a2)))
(a0*X1^((a1+1)/a1)*X2^(1/a2)*log(X2))/(2*a2^2*(a3^X3+X2^(1/a2))^2*sqrt(X1^((a1+1)/a1)/(a3^X3+X2^(1/a2))))
-(a0*a3^(X3-1)*X1^((a1+1)/a1)*X3)/(2*sqrt(X1^((a1+1)/a1)/(a3^X3+X2^(1/a2)))*(a3^X3+X2^(1/a2))^2)
genr a0 = 1
genr a1 = 6
genr a2 = 2
genr a3 = 2
y=a0*sqrt((X1^((1+a1)/a1))/(X2^(1/a2)+a3^X3))
params a0 a1 a2 a3