function pp_plot_para(r1,K,n01,c,a,Th,d,n02,paratochange,inival,endval,stepval,animation)
%This function plots simulated dynamics and isolines together (predator-prey)

%Prey:
% 1 r1 - intrinsic growth rate 
% 2 K - prey carrying capacity
% 3 n01 - initial population sizes 

%Predator:
% 4 c - prey conversion rate
% 5 a - encounter rate
% 6 Th - handling time
% 7 d - predator mortality
% 8 n02 - initial population size

%further parameters:
% paratochange: number of the para to be changed (see above)
% inival: initial parameter value
% endval: final parameter value
% stepval: parameter change step

%animation - whether to animate

para(1) = r1; 
para(2) = K; 
para(4) = c; 
para(5) = a;
para(6) = Th;
para(7) = d;
paraname = ['prey growth rate      ';'prey carrying capacity';'                      ';'prey conversion rate  ';'encounter rate        ';'handling time         ';'predator mortality    ']';

xx=[];
yy=[];

for value = inival:stepval:endval
   para(paratochange) = value;
   [y,R0]=pp_lines(para(1),para(2),para(4),para(5),para(6),para(7));
   xx =[xx y(1,:) R0];
   yy =[yy y(2,:)];

   if (R0 == 0)
      return;
   end
end
maxx = max(xx);
maxy = max(yy);

hold on;  
   
for value = inival:stepval:endval
   para(paratochange) = value;
   [y,R0]=pp_lines(para(1),para(2),para(4),para(5),para(6),para(7));
   text = sprintf('%s %7.2g ', paraname(:,paratochange), value);
   plot(y(1,:),y(2,:));
   axis([0, maxx, 0, maxy])
   title(text);
   line([R0 R0 ],[0 maxy],'color','g');
   if (animation > 0)
      pause(1);
   end
   
end