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

/* [wxMaxima: input   start ] */
kill(all);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
definisco funzioni spostamento genericamente lineari
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
u(x,y) := a1 + a2*x + a3*y;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
v(x,y) := a4 + a5*x + a6*y;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]

   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
eqns : [
  u(x_i,y_i) = u_i , v(x_i,y_i)=v_i ,
  u(x_j,y_j) = u_j , v(x_j,y_j)=v_j ,
  u(x_k,y_k) = u_k , v(x_k,y_k)=v_k
];
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
concatenazione: 
eqns : append(  eqns , [ ANCH_IO , PURE_IO] );
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
unkn : [ a1 , a2 , a3 , a4 , a5 , a6 ];
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
ho definito un sistema lineare di equazioni, che provvedo
a risolvere 
   [wxMaxima: comment end   ] */

/* [wxMaxima: comment start ]
linsolve( eqns , unkn ) , globalsolve=true;
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
sols : linsolve( eqns , unkn ) , globalsolve=true;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
assegnazione "manuale" senza globalsolve=true
a1 : ev (a1 , sols);
   [wxMaxima: comment end   ] */

/* [wxMaxima: comment start ]

   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
epsx : diff ( u(x,y) , x , 1 );
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
epsy : diff ( v(x,y) , y , 1 );
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
gammaxy : diff ( v(x,y) , x , 1 ) + diff ( u(x,y) , y , 1 );
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
creo un vettore (lista) eps
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
eps : [ epsx , epsy , gammaxy ] ;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
controllo la natura di costante in x,y del di cui sopra
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
diff ( eps , x , 1);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
diff ( eps , y , 1);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
caso soluzioni multiple (sistema di equazioni di secondo grado)
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
sols:solve( [x^2 + x*y + 1 = 0, y^2 + 2*x -4=0],[x,y]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
assegno soluzione 1
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
[x,y] : ev ([x,y] , sols[1]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
x;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
y;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
dimentico soluzione 1
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
kill(x,y);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
sostituisco soluzioni
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
[x,y] : ev ([x,y] , sols[2]);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
dimentico soluzione 2
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
kill(x,y);
/* [wxMaxima: input   end   ] */

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