i have problem error code

[fname,pname] = uigetfile('*.lvm');
rawdata=load (fname);
for j1=1:length(rawdata)/31;
time(j1,1)=(rawdata(31*j1,2)-rawdata(31,2))/1000;
for j2=1:30;
thirdoctave(j1,j2)=rawdata((j1-1)*31+j2,2);
end
end
[maxdb,idb]=max(thirdoctave);
for j=3:8;
j1=j*3;
j2=idb(j1);
j3=1;
flag=0;
j5=1;
end
while flag==0;
if j3>3 & j2+1<=length(thirdoctave) & thirdoctave(j2+1,j1)>thirdoctave(j2,j1)
flag=1;
else if j2 == length(time);
flag=1;
else
cropdata(j3,:,j1)=[time(j2),thirdoctave(j2,j1)];
if j3>1;
for j4=1:j3-1
slope(j5,j1)=(cropdata(j3,2,j1)-cropdata(j3-j4,2,j1))/(cropdata(j3,1,j1)-cropdata(j3-j4,1,j1));
j5=j5+1;
end
end
j2=j2+1;
j3=j3+1;
end
end
decay(j-2,j2)=mean(slope(1:j5-1,j1));
end
freqs=[20;25;31.5;40;50;63;80;100;125;160;200;250;315;400;500;630;800;1000;1250;1600;2000;2500;3150;4000;5000;6300;8000;10000;12500;16000];
freqcrop=[125;250;500;1000;2000;4000];
reverb=-60./decay;
this code matlab

and data file in .lvm it example data

http://www.mediafire.com/?5dzmd94zgeg81vl here

i want measured reverberation time(RT60) use in labview and calculate in matlab

i follow this http://scitation.aip.org/confst/ASA/data/11/4pAA6.pdf

i calculate code matlab error