% 1-simplex
syms x1 t1 t2 positive;
F = x1*t1 / (x1*t1 + (1-x1)*t2);
F= sqrt(F);
J = jacobian(F,x1);
g = J*J + (1-sum(J))^2;
 
