ADDED mttroot/mtt/lib/examples/Chemical/ReactorTF/MakeFigure.m Index: mttroot/mtt/lib/examples/Chemical/ReactorTF/MakeFigure.m ================================================================== --- /dev/null +++ mttroot/mtt/lib/examples/Chemical/ReactorTF/MakeFigure.m @@ -0,0 +1,22 @@ +## Makes the figures + +par = ReactorTF_numpar; # Parameters +sym = Reactor_sympar; # Parameter indices + +F_s= [90:10:500]; # Range of flows + +Z = []; +for f_s=F_s + par(sym.f_s) = f_s; + z = sort(eig(ReactorTF_sm(par))); + Z = [Z z]; +endfor + + +grid; xlabel("f_s"); ylabel("z_1"); +plot(F_s,Z(1,:)); +psfig("ReactorTF_zero_1"); + +grid; xlabel("f_s"); ylabel("z_2"); +plot(F_s,Z(2,:)); +psfig("ReactorTF_zero_2"); ADDED mttroot/mtt/lib/examples/Chemical/ReactorTF/Makefile Index: mttroot/mtt/lib/examples/Chemical/ReactorTF/Makefile ================================================================== --- /dev/null +++ mttroot/mtt/lib/examples/Chemical/ReactorTF/Makefile @@ -0,0 +1,14 @@ +## Makes the schematic diagram and the zero figure +all: ReactorTF_pic.ps ReactorTF_zero_1.ps ReactorTF_zero_2.ps + +ReactorTF_pic.ps: ReactorTF_pic.fig + fig2dev -Lps ReactorTF_pic.fig> ReactorTF_pic.ps + +ReactorTF_zero_1.ps: ReactorTF_abg.fig + mtt -q ReactorTF sm m; mtt -q ReactorTF numpar m; + mtt -q Reactor sympar m; + octave MakeFigure.m + +ReactorTF_zero_2.ps: ReactorTF_zero_1.ps + touch ReactorTF_zero_2.ps +