%T7_B4_Isaac_Paha

load pellets

Mean_value1=mean(pellets)
S1=std(pellets)
x=1:30;

figure(1)
errorbar(x,Mean_value1,S1)  
xlabel('the number of days')
ylabel('the number Pellets')

Mean_value2=mean(pellets,2) 
S2=std(pellets,0,2) 
figure(2)
errorbar(x,Mean_value2,S2)
xlabel('the number of animals')
ylabel('the number of pellets')



SEM1=S1/sqrt(30)

SEM2=S2/sqrt(30)






