% cstr_cusp.m
% generates the cusp diagram, based on case 2 parameters
%
  case2;
%
clear Tempr Tjs Tfs concs;
   Tempr = 280:5:400;
% vary UA or F
   npts = length(Tempr);
   CSTR_PAR(8) = 0.1;
   F_plt = ones(npts)*CSTR_PAR(8);
   [Tjs,Tfs,concs] = cstr_io(Tempr);
   plot3(F_plt,Tjs,Tempr,'w-')
   hold on
 for kindex = 1:9;
   CSTR_PAR(8) =  0.1+ 0.1*kindex;
   F_plt = ones(npts)*CSTR_PAR(8);
   [Tjs,Tfs,concs] = cstr_io(Tempr);
   plot3(F_plt,Tjs,Tempr,'w-')
 end
   hold off
%   axis([0.25 1.375 260 360 300 380])
   axis([0.1 1.0 220 400 280 400])
   view(60,20)
   xlabel('F/V')
   ylabel('Temp_j')
   zlabel('Temp_r')