% plotslnf.m 

%plotslnf.m 

function [tplot,ftxslnplot]=plotslnf(ftxsln,ax,tf,x0) 
% Used by SYS3PLOT.M to plot the user-derived solution to an ODE system. 
% User enters solutions via INITSYS.M 

axis(ax); 
plotN=50; 
t=ax(1):(ax(2)-ax(1))/(plotN-1):ax(2); 
ttemp=t; xmin=x0; xmax=x0; tf=ax(2); 
t=tf; 
if eval(ftxsln)>-inf 
	xmin=min([x0,eval(ftxsln)]); 
	if eval(ftxsln)<inf 
		xmax=max([x0,eval(ftxsln)]); 
		ftxslnplot=[]; 
		t=[]; 
		for j=1:length(ttemp) 
			t=ttemp(j); 
			ftxslnplot=[ftxslnplot,eval(ftxsln)]; 
			if min(ftxslnplot)>-inf 
				if xmin>min(ftxslnplot) 
					xmin=min(ftxslnplot); 
				end 
			end 
			if max(ftxslnplot)<inf 
				if xmax<max(ftxslnplot) 
					xmax=max(ftxslnplot); 
				end 
			end 
		end 
		t=ttemp; 
		tplot=ttemp; 
	else 
		disp('Your solution has a + inf result somewhere') 
		break 
	end 
else 
	disp('Your solution has a - infinity result somewhere') 
	break 
end