function out=plotLorenz(T,dt,rho,Ics)
%x0=plotLorenz(100,0.01,5,[0.1,0.1,3])
%x0=plotLorenz(100,0.01,24.1,x0)
x1=Ics(1); x2=Ics(2); x3=Ics(3);
sig=10.0;
bet=8/3;
% integrate forwards T units
[t,x] = ode45(@(t,x) lorenz_par(t,x,rho),[0:dt:T],[x1;x2;x3]);
n=length(t);
% show trajectory being analyzed
plot(x(:,1),x(:,3),'-','MarkerSize',2);
%figure
%plot3(x(:,1),x(:,2),x(:,3),'-','MarkerSize',2);
out=x(end,:);
