% define input arguments fs, wallrefl, r, src, mic, n, r

[h]=my_ISM(fs,mic,n,wallrefl,rm,src);  % generates RIR using my_ISM.m function

[decay_curve,EDT,T60]=IR_Analyze(h,fs); % plots Schroeder decay curve of RIR
                                        % and generates EDT & T60 values
                                        
% to implement the RIR 
[audio,fs]=audioread('anechoic_voice.wav');
[y]=fconv(audio,h);