f=freq_el(:,10);
min(f)
max(f)
sr=std(f);
mr=mean(f);
edges=0.0026:0.001:1.75*10^-2;
x_axis=3.1*10^-3:0.001:1.8*10^-2;
n=histc(f,edges);
bar(x_axis,n,1)
hh=bar(x_axis,n,1)
np = get(gca,'NextPlot');
set(gca,'NextPlot','add')
x=(-3*sr+mr:0.1*sr:3*sr+mr)';
y = normpdf(x,mr,sr);
y=197*(y*0.001);
hh1 = plot(x,y,'r-','LineWidth',2);
edges=0.9:0.5:7;
n=histc(v_cross,edges);
x_axis=0.925:0.5:7.25;
hh=bar(x_axis,n,1)
sr=std(v_cross);
mr=mean(v_cross);
x=(-3*sr+mr:0.1*sr:3*sr+mr)';
y = normpdf(x,mr,sr);
y=197*(y*0.5);
np = get(gca,'NextPlot');
set(gca,'NextPlot','add')
hh1 = plot(x,y,'r-','LineWidth',2);