/* [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 ]
istruzione ratsubst 

ratsubst(
  sostituto,
  da_sostituirsi,
  espressione_entro_cui_operare
)
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
expr: (a*x^4 + b*x^2 + c);/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
ratsubst(
y,
x^2,
expr
);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
istruzioni assume e forget
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
sqrt(a^2*b^2);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
assume(a>0,b>0);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
sqrt(a^2*b^2);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
forget(a>0);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
sqrt(a^2*b^2);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
istruzione makelist
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
x:makelist(  pippo^2  ,  pippo  ,1,20,2);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
makelist(  %[i],      i, 1 , length(%) , 2 );
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
y:makelist ( ele/(ele+1) , ele , x );
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
grafici per punti
   [wxMaxima: comment end   ] */

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

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

/* [wxMaxima: comment start ]
manipolazioni trigonometriche
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
sin(t)^2+cos(t)^2;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
trigsimp(%);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
formula in coordinate polari r,t, con parametro alpha 
(es. fase)
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
sin(3*t+alpha);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]

   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
trigexpand(%);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
trigexpand(%);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
expand(%);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
ne ottengo una forma in coordinate cartesiane
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
ratsubst(x/r,cos(t),%)$
ratsubst(y/r,sin(t),%);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]

   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
x+sin(3*x)/sin(x);
trigexpand(%),trigexpandtimes=true,expand;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
equazioni non lineari e numeri complessi
   [wxMaxima: comment end   ] */

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

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

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

/* [wxMaxima: input   start ] */
realpart(y);
imagpart(y);
/* [wxMaxima: input   end   ] */

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

/* [wxMaxima: comment start ]
da forme trigonometriche a esponenziali complessi
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
a*cos(t)+b*%i*sin(t);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
exponentialize(%);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
demoivre(%);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]

   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
a*cos(y)+b*%i*sin(x);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
exponentialize(%);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
expand(%);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
demoivre(%);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
fullratsimp(%);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
blocchi e variabili locali ai blocchi
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
tmp:0;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
a:pippo;
b:pluto;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
block([tmp], tmp:a , a:b, b:tmp );
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
a;b;
/* [wxMaxima: input   end   ] */

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

/* [wxMaxima: input   start ] */
kill(a,b,tmp);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
integrazione numerica
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
quad_qag( sin(x) , x,0, %pi ,1  );
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
operazioni si matrici
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
A:zeromatrix(3,3);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
A[1,2]:a $
A[1,3]:b $
A[2,3]:c $
/* [wxMaxima: input   end   ] */

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

/* [wxMaxima: input   start ] */
A : A - transpose(A);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
v:  matrix( [x], [y], [z] );
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
transpose(v).A.v;
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
fullratsimp(%);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
genero matrice con 1 su parte sopradiagonale,
0 sulla diagonale, -1 nella parte sottodiagonale
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
A:genmatrix(
lambda( [i,j] ,   

if j>i then 
    1   
    elseif j<i then 
        -1 
else 
    a

)
,3
,3
);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
determinant(A);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
solve(%,a);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
da funzione senza nome a funzione con nome mediante
assegnazione a nome
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
fun : lambda( [i,j] , i + %i*j  );
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
fun(3,4);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]

   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
B:augcoefmatrix(
[
    a*x1+b*y1+c=f1,
    a*x2+b*y2+c=f2,
    a*x3+b*y3+c=f3
],
[ a,b,c]
);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
coefmatrix(
[
    a*x1+b*y1+c=f1,
    a*x2+b*y2+c=f2,
    a*x3+b*y3+c=f3
],
[ a,b,c]
);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: input   start ] */
determinant(%);
/* [wxMaxima: input   end   ] */

/* [wxMaxima: comment start ]
cambio di segno del termine noto nella matrice aumentata
(da stackoverflow)
   [wxMaxima: comment end   ] */

/* [wxMaxima: input   start ] */
transpose(
    append(
        transpose(
            submatrix(B,4)
        )
        ,
        -transpose(
            submatrix(B,1,2,3)
        )
    )
);
/* [wxMaxima: input   end   ] */

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