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

/* [wxMaxima: comment start ]
pulisco memoria
   [wxMaxima: comment end   ] */

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

/* [wxMaxima: comment start ]
programma che calcola i coefficienti della forma lineare degli spostamenti
in un elemento tria3
   [wxMaxima: comment end   ] */

/* [wxMaxima: comment start ]
definisco generica forma lineare per gli spostamenti u e v
   [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 ]
lista di equazioni che legano i valori di u,v in corrispondenza
dei nodi agli spostamenti nodali
   [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: input   start ] */
unkn : [ a1,a2,a3,a4,a5,a6 ];
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
ho un sistema di equazioni, procedo a risolverlo
   [wxMaxima: comment end   ] */

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

/* [wxMaxima: comment start ]
campo spostamenti
   [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 );
/* [wxMaxima: input   end   ] */

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

/* [wxMaxima: comment start ]
verifico natura costante in x,y delle espressioni di cui sopra
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
diff( [ epsx,epsy,gammaxy ] , x);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
diff( [ epsx,epsy,gammaxy ] , y);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
vettore (lista) contenente le componenti di deformazione
   [wxMaxima: comment end   ] */

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

/* [wxMaxima: comment start ]
utilizzo di coefmatrix
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
delta : [u_i,v_i,u_j,v_j,u_k,v_k];
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
coefmatrix(eps,delta);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
ratsubst ( duearea  ,  (x_j-x_i)*y_k+(x_i-x_k)*y_j+(x_k-x_j)*y_i  , %  );
/* [wxMaxima: input   end   ] */

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