Peak_Fe.m

n=20;
%Use uniform77 as axon distribution 1:0.25:20=77
b=zeros(77,n);
c=zeros(77,n);
for j=1:n
    neuroFEM('Apply','neuroFEM_ana')
    a=get(findobj('tag','ZengAnalysis_NeuroFEM'),'userdata');
    j
    for i=1:77
        temp=max(diff(a.Group.vestim(1:a.Group.nseg(i),i),2));
        b(i,j)=temp(1);
        temp=min(diff(a.Group.vestim(1:a.Group.nseg(i),i),2));
        c(i,j)=temp(1);
    end
end
figure
plot(b)
hold on;plot(c)
dlmwrite('s700-15g11v-Depo-Fe-E0.txt',[[1:0.25:20]' b],'\t')
dlmwrite('s700-15g11v-Hyper-Fe-E0.txt',[[1:0.25:20]' c],'\t')