load m81vC.strip
t=m81vC(1:16,1)-2448620;
c11=m81vC(1:16,2);
ec11=m81vC(1:16,3);
c12=m81vC(1:16,4);
ec12=m81vC(1:16,5);
c13=m81vC(1:16,6);
ec13=m81vC(1:16,7);
c14=m81vC(1:16,8);
ec14=m81vC(1:16,9);
c15=m81vC(1:16,10);
ec15=m81vC(1:16,11);
c16=m81vC(1:16,12);
ec16=m81vC(1:16,13);
subplot(2,3,1)
plot(t,c11,'o',t,c11)
adderror(t,c11,ec11)
title('C11')
subplot(2,3,2)
plot(t,c12,'o',t,c12)
adderror(t,c12,ec12)
title('C12')
subplot(2,3,3)
plot(t,c13,'o',t,c13)
adderror(t,c13,ec13)
title('C13')
subplot(2,3,4)
plot(t,c14,'o',t,c14)
adderror(t,c14,ec14)
title('C14')
subplot(2,3,5)
plot(t,c15,'o',t,c15)
adderror(t,c15,ec15)
title('C15')
subplot(2,3,6)
plot(t,c16,'o',t,c16)
adderror(t,c16,ec16)
title('C16')