TestData = randn(1000,3); 
for k=1:size(TestData,2);
    
mean_new(k,:) = mean(TestData(:,k));    
std_big(k,:)  = mean(TestData(:,k))  +  std(TestData(:,k));
std_low(k,:)  = mean(TestData(:,k))  -  std(TestData(:,k));

datavec = TestData(:,k);
figure(k)
x = 1:length(datavec);

hold on
plot(datavec (std_low(k,1) < datavec  & datavec < std_big(k,1)) ,'b.')

 
line(xlim,[mean_new(k,1),mean_new(k,1)],'Color','green')
line(xlim,[std_low(k,1),std_low(k,1)],'Color','red')
line(xlim,[std_big(k,1),std_big(k,1)],'Color','red')


 
plot(datavec (datavec  > std_big(k,1)),'r.')
plot(datavec (datavec  < std_low(k,1)),'r.')

lgd= legend (['Data Points'],['Ideal = ' num2str(mean_new(k,1))],['UpperLimit = ' num2str(std_big(k,1))], ['LowerLimit = ' num2str(std_low(k,1))]  );
set(lgd,'Location','northEastoutside','FontSize',12);
hold off
end