set(hp1,'xdata',bbb(1,:),'ydata',bbb(2,:),'zdata',bbb(3,:));
    M=size(bbb,2);
    name=strcat('laser_',num2str(i),'.txt');
    fid=fopen(name,'a+');
    for kk =1:M
        fprintf(fid,'%g %g %g\r\n',bbb(1,kk),bbb(2,kk),bbb(3,kk));
    end
    fclose(fid);

lineSeg=[1;1];
lineSeg=[lineSeg [1; 1]];

function  DrawBoundingBox(xBnd,yBnd,zBnd,lineWidth,color)
% DrawBoundingBox(xBnd,yBnd,zBnd,lineWidth,color)
%
% This function draws the wireframe box that as described by the limits in
% xBnd, yBnd, and zBnd
%
% INPUTS:
%   xBnd = [xLow, xUpp]
%   yBnd = [yLow, yUpp]
%   zBnd = [zLow, zUpp]
%

hold on;

% Draw the bottom:
plot3(...
    xBnd([1,1,2,2,1]),...
    yBnd([1,2,2,1,1]),...
    zBnd([1,1,1,1,1]),...
    'LineWidth',lineWidth','color',color);

% Draw the top:
plot3(...
    xBnd([1,1,2,2,1]),...
    yBnd([1,2,2,1,1]),...
    zBnd([2,2,2,2,2]),...
    'LineWidth',lineWidth','color',color);

% Draw the sides:
plot3(...
    xBnd([1,1]),...
    yBnd([1,1]),...
    zBnd([1,2]),...
    'LineWidth',lineWidth','color',color);
plot3(...
    xBnd([1,1]),...
    yBnd([2,2]),...
    zBnd([1,2]),...
    'LineWidth',lineWidth','color',color);
plot3(...
    xBnd([2,2]),...
    yBnd([2,2]),...
    zBnd([1,2]),...
    'LineWidth',lineWidth','color',color);
plot3(...
    xBnd([2,2]),...
    yBnd([1,1]),...
    zBnd([1,2]),...
    'LineWidth',lineWidth','color',color);

end

plot3(model(1,:),model(2,:),model(3,:),'r.',data(1,:),data(2,:),data(3,:),'b.','MarkerSize',1), hold on, axis equal