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);    % Plot the system    figure(1);    plothex(xyb,fib);    hold on;    plot(xym,'b*');    hold off;