% myrandom.m  
% Units are mm.

% domain size
bx = 160;
by = 40;
bz = 40;

lf = 12; % rod length
nf = 500; % Number of rods

rns = rand(nf,3);      % Start
rne = rand(nf,3)-0.5;  % End

% Start Points
for i = 1:nf
    rns(i,1) = rns(i,1)*bx;
    rns(i,2) = rns(i,2)*by;
    rns(i,3) = rns(i,3)*bz;
end

% Unit Deltas
delta = zeros(nf,1);
for i = 1:nf
   temp = rne(i,:);
   delta(i) = norm(temp);
end

% Length Deltas
rne = lf*rne./delta;

% End Points
rne = rns + rne;


fileID = fopen('scfibers.txt','w');
for i = 1:nf
    fprintf(fileID,'%12.8f %12.8f %12.8frn',rns(i,1),rns(i,2),rns(i,3));
    fprintf(fileID,'%12.8f %12.8f %12.8frnrn',rne(i,1),rne(i,2),rne(i,3));
end
fclose(fileID);