function []=vector_field(fun,x1ray,x2ray,tailend);
% Carol Lucas
% vector_field.m
% plots initial 2 seconds for a variety of initial conditions

clear t;
mess='Initial point indicated with O'
t=[0:.01:tailend];
lt=length(t);

for n=1:length(x1ray);
	xinit(1,1)=x1ray(n);
for m=1:length(x2ray);
	xinit(2,1)=x2ray(m);
	
	[t,x]=ode45(fun,t,xinit);
	plot(x(1,1),x(1,2),'co',x(:,1),x(:,2),'c');pause(.5);
	end
	end