%
%	THE ASSUMED PRAMETER FOR HUMAN COCHLEA	
%
%	plots the mean radius of scala vestibuli and scala tympani of
%	a human cochlea, then the basilar membrane width, loading
%	data from BM_Width.MAT and xpandyp.MAT


xp=0;
yp=0;
xm=0;
ym=0;
load xpandyp
bmL=33.5; % BM length in humans (mm) = BETA unit

bmw_man=0;
load BM_Width.mat -ascii  % data are in BETA units 

x=bmL*BM_Width(:,1);	
y=bmL*BM_Width(:,2);

figure(1)

subplot(2,1,1)
axis([0,33.5,-2,2]);
plot(xp*bmL, yp*bmL, xm*bmL, -ym*bmL, (0:.05:1)*bmL, zeros(size(0:.05:1)),'w')

title('Radius of human cochlea (scala vestibuli & scala tympani) in mm')
xlabel('Distance from stapes (mm)') 
ylabel('Radius (mm)') 

subplot(212)
axis([0, 33.5, 0, 0.6]);
plot(x,y,'-red',x,y,'-g')
title('Basilar membrane width characteristic in humans');
ylabel('Witdh in (mm)');
xlabel('Length in (mm)') 
