clear; %%% Plot data from newtn.f %%%
load xplot.dat
xplot = xplot.'; % Transpose matrix
temp = size(xplot);
subplot(121)
plot(xplot(1,:),xplot(2,:),'o',...
xplot(1,:),xplot(2,:),'-',xplot(1,temp(2)),xplot(2,temp(2)),'*');
xlabel('x'); ylabel('y');
subplot(122)
plot(xplot(1,:),xplot(3,:),'o',...
xplot(1,:),xplot(3,:),'-',xplot(1,temp(2)),xplot(3,temp(2)),'*');
xlabel('x'); ylabel('z');
subplot(111)
