/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/
/* [ Created with wxMaxima version 15.04.0 ] */

/* [wxMaxima: input   start ] */
kill(all)$
N:6$
t:makelist(concat(x_,i),i,0,2)$
p:makelist(concat(x_,i),i,3,5)$
x:makelist(concat(x_,i),i,0,N-1)$
f[0]:(x_0-x_3)^2+(x_1-x_4)^2+(x_2-x_5)^2$
f[1]:apply(F[0],t)$
f[2]:apply(F[1],t)$
f[3]:apply(F[2],p)$
f[4]:apply(F[3],p)$
g:[]$
for i:0 thru N-2 do
g:endcons(f[i],g)$
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
Jac:jacobian(g,x);Jac:transpose(Jac);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
A:submatrix(N,Jac);b:-Jac[N];
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
sol[N-1]:determinant(A)$
for i:1 thru N-1 do
(temp:A[i], A[i]:b, sol[i-1]:(determinant(A)), A[i]:temp)$
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
f1:('diff(F[0](x_0,x_1,x_2),x_1,1))*('diff(F[1](x_0,x_1,x_2),x_2,1))-('diff(F[0](x_0,x_1,x_2),x_2,1))*
('diff(F[1](x_0,x_1,x_2),x_1,1))=D1;/*10*(x_1-x_0)*/
f2:('diff(F[0](x_0,x_1,x_2),x_2,1))*('diff(F[1](x_0,x_1,x_2),x_0,1))-('diff(F[0](x_0,x_1,x_2),x_0,1))*
('diff(F[1](x_0,x_1,x_2),x_2,1))=D2;/*(x_0*(28-x_2)-x_1)*/
f3:('diff(F[0](x_0,x_1,x_2),x_0,1))*('diff(F[1](x_0,x_1,x_2),x_1,1))-('diff(F[0](x_0,x_1,x_2),x_1,1))*
('diff(F[1](x_0,x_1,x_2),x_0,1))=D3;/*(x_0*x_1-8*x_2/3)*/
f4:('diff(F[2](x_3,x_4,x_5),x_4,1))*('diff(F[3](x_3,x_4,x_5),x_5,1))-('diff(F[2](x_3,x_4,x_5),x_5,1))*
('diff(F[3](x_3,x_4,x_5),x_4,1))=D4;
f5:('diff(F[2](x_3,x_4,x_5),x_5,1))*('diff(F[3](x_3,x_4,x_5),x_3,1))-('diff(F[2](x_3,x_4,x_5),x_3,1))*
('diff(F[3](x_3,x_4,x_5),x_5,1))=D5;
f6:('diff(F[2](x_3,x_4,x_5),x_3,1))*('diff(F[3](x_3,x_4,x_5),x_4,1))-('diff(F[2](x_3,x_4,x_5),x_4,1))*
('diff(F[3](x_3,x_4,x_5),x_3,1))=D6;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
for i:0 thru N-1 do
(print(factorsum(scsimp(ratsimp(sol[i]),f1,f2,f3,f4,f5,f6))));
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
for i:0 thru N-1 do
(printf(true, "f[~a]=~a;~%",i,factorsum(scsimp(ratsimp(sol[i]),f1,f2,f3,f4,f5,f6))));
/* [wxMaxima: input   end   ] */

/* Maxima can't load/batch files which end with a comment! */
"Created with wxMaxima"$