Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- % 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);
Add Comment
Please, Sign In to add comment