% plotsln.m 
% Used to plot an explicit solution to an ODE. Solution is defined via 
% INIT.M. (Can actually be used to plot any function defined in terms of 
% the independent variable t 

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; 
		plot(t,ftxslnplot) 
		xlabel(['x = ',ftxsln]) 
		grid on 
	else 
		disp('Your solution has a + inf result somewhere') 
		break 
	end 
else 
	disp('Your solution has a - infinity result somewhere') 
	break 
end