[t,bz4,bz3,bz2,bz1,bz0] = textread('vsf.dat','%f %f %f %f %f %f ','headerlines',1)
bz4bz3bz2bz1bz0Figure = figure
axes('FontSize',14,'position',[0.15 0.15 0.75 0.75])
plot(t,bz4,'--.g',t,bz3,'.r',t,bz2,'.-b',t,bz1,'.--c',t,bz0,'.-.m','MarkerSize',10.0)
axis([-0.5 2 -25 30])
bz4XLabel=(t(length(t)))/2;
bz4YLabel=bz4((length(t))/2);
text(bz4XLabel,bz4YLabel,'Bz_4 \rightarrow ','HorizontalAlignment','right','FontSize',14)
bz3XLabel=t(1);
bz3YLabel=bz3(1);
text(bz3XLabel,bz3YLabel,'Bz_3 \rightarrow ','HorizontalAlignment','right','FontSize',14)
bz2XLabel=t(1);
bz2YLabel=bz2(1);
text(bz2XLabel,bz2YLabel,'Bz_2 \rightarrow ','HorizontalAlignment','right','FontSize',14)
bz1XLabel=(t(length(t)))/4;
bz1YLabel=bz1((length(t))/4);
text(bz1XLabel,bz1YLabel,'Bz_1 \rightarrow ','HorizontalAlignment','right','FontSize',14)
bz0XLabel=t(length(t));
bz0YLabel=bz0(length(t));
text(bz0XLabel,bz0YLabel,' \leftarrow Bz_0 ','HorizontalAlignment','left','FontSize',14)
TbzXLabel=1;
TbzYLabel=-22;
text(TbzXLabel,TbzYLabel,'t_b_r=1.0 s ','HorizontalAlignment','center','FontSize',14)
grid on
xlabel('Time, s','FontSize',14)
ylabel('Bz, mT','FontSize',14)
title('Vertical stray field','FontSize',14)
print -f1 -djpeg 'vsf.jpeg'
