function tutorial    % Create the base stations    par.cellradius=1000;    par.sps=1;    par.km=1;    par.lm=2;    par.ncluster=2;    [xyb,fib,rombvec]=crecells(par.cellradius,par.sps,par.km,par.lm,par.ncluster);    % Create new mobiles    par.vmean=10;    par.amean=5;    par.dt=10;       nmob=30;    xym=nans(nmob,1);    xyv=nans(nmob,1);    m=(1:nmob).';    [xym,xyv]=mobmove(xym,xyv,par.vmean,par.amean,par.dt,rombvec);    % Create the antenna lobe    lobevector=zeros(360,1);    % Calculate the gain between bases and mobiles    par.attconst=-21;    par.alpha=4;    par.sigma=8;    par.raa=0.5;    par.corrdist=110;    [lognmap, mapvec] = crelognmap(xyb, rombvec, par.corrdist);    g=pathgain(xym,xyb,fib,rombvec,par.attconst,par.alpha, par.sigma, par.raa, lobevector, lognmap, mapvec);    % Plot the system    figure(1);    plothex(xyb,fib);    hold on;    plot(xym,'b*');    hold off;