               #Definizioni Globali
set term latex;
unset label;
unset border;
unset title;
               #Livelli di Fermi
unset xtics;
g(x,kT,a) = 1/(exp((x-a)/kT)+1); 
set data style points;
set function style lines;
set style line 1 lt 1 lw 1 pt 1 ps 2;
set style line 2 lt 1 lw 1 pt 1 ps 2;
set style line 3 lt 1 lw 1 pt 1 ps 2;
set style line 4 lt 1 lw 1 pt 1 ps 2;
set style line 5 lt 1 lw 1 pt 1 ps 2;
set style line 6 lt 1 lw 1 pt 1 ps 2;
set size .6,.8;
set label 1 '$\E_F$' at 1.9,-.05;
set label 2 '$\T_1$' at 2.05,0.05;
set label 3 '$\T_2$' at 2.4,.2;
set label 4 '$\T\simeq0$' at 1.4,.1;
set xlabel '$\E$';
set ylabel '$\mez{n}$';
set output 'FermiaOkupado.tex';
plot [0:3][0:1.3]g(x,0.01,2) title '',g(x,0.1,2) title '', g(x,0.2,2) title '';
unset label
                     #Boltzmann Ekvacio sur Fermi Distribuo
set output 'BoltzmannTrasporto.tex'
set label 1 '$f_0$' at 1.8,.5;
set label 2 '$f$' at 2.2,.5;
set size .4,.4;
unset ytics;
plot [0:3][0:1.3]g(x,0.1,2) title '',g(x,0.1,2.5) title '';
f(x)=g(x,0.1,2.5)-g(x,0.1,2);
unset label;
set output 'DefDeT.tex';
set size .4,.4;
set label 1 '$\De\T f$' at 1.5,.5;
plot [0:3][0:1.3] f(x) title '';
unset xlabel;
unset ylabel;
unset ytics;
unset label;
                   #Polarizajxo
set size 1,1;
unset border;
set style line 1 lt 1 lw 1 pt 1 ps 2;
f(x)=e1+a*(e2-e1)/(a**2-x**2);
g(x)=b/((a**2-x**2)**2-g);
e1=.1;e2=.1;a=1.;b=1.;g=1.;
set output 'Polarizajxo.tex';
plot  [0:100][-3:3] f(x) title '',g(x) title '';
                   #Funzione Sinc
g(x)=sin(x)/x;
set arrow 1 from -20,0 to 20,0 nohead;
set term texdraw;
set output 'Sinc.tex';
set size .5,.5;
plot [-20:20] g(x) title'';
unset arrow;
                     #Paketo
g(x) = x<-4 ? exp(x*.5+2)*sin(x*3.5) : x>4? exp(-x*.5+2)*sin(3.5*x) : sin(3.5*x);
set output 'Paketo.tex';
set size 1,.4;
plot [-10:10] g(x) title '';
                      #Batado
g(x) = sin(x)+sin(x*.8);
f(x) = 2*cos(x*.1);
set style line 2 lt 4 lw 4 pt 1 ps 1;
set style line 3 lt 4 lw 4 pt 1 ps 1;
set output 'Batado.tex';
set size 1,.4;
plot [-50:50] g(x) title '',f(x) title '',-f(x) title '';
                      #Bose-Einstein-a kondenso
f(x) = x + x**2/(2**(1.5)) - x/(x-1);
set xlabel '$z$'
set ylabel '$\Lambda^3\no$';
set output 'BEKondenso.tex';
set size .4,.4;
set arrow 1 from 1,0 to 1,10 nohead;
plot [0:1][0:10] f(x) title'';
unset arrow
                    #Faza transado
f(x) = 1 - (x/273)**(1.5);
set xlabel '$\T$';
set ylabel '$\frac{\mez{n_0}}{\mez N}$';
set output 'FazaTransado.tex';
plot [0:300][0:1.2] f(x) title '';
unset xlabel;
unset ylabel;

