[lderesult,ri,rj,lr]=lde(ri,rj,lr);

function [lderesult,ri,rj,lr]=lde(ri,rj,lr);
lderesult=[];
persistent i ; 

;
ri_shape=size(ri);ri=reshape(ri,1,[]);
rj_shape=size(rj);rj=reshape(rj,1,[]);
if isempty(i), i=0; end;
lderesult = true;
for i = 1 : lr;
if( ri(i)~=rj(i) )
lderesult = false;
break;
end;
end;
ri_shape=zeros(ri_shape);ri_shape(:)=ri(1:numel(ri_shape));ri=ri_shape;
rj_shape=zeros(rj_shape);rj_shape(:)=rj(1:numel(rj_shape));rj=rj_shape;
csnil=dbstack(1); csnil=csnil(1).name(1)~='@';
if csnil&&~isempty(inputname(2)), assignin('caller','FUntemp',rj); evalin('caller',[inputname(2),'=FUntemp;']); end
if csnil&&~isempty(inputname(1)), assignin('caller','FUntemp',ri); evalin('caller',[inputname(1),'=FUntemp;']); end
if csnil&&~isempty(inputname(3)), assignin('caller','FUntemp',lr); evalin('caller',[inputname(3),'=FUntemp;']); end
end %function lde