function plot_fjc_axes(fname)

A = load(fname);
points = A(:,1:3);
dirs = A(:,4:6);
startpoints = points + 1 * dirs;
endpoints = points - 1 * dirs;

lines.x = [startpoints(:,1)'; endpoints(:,1)'];
lines.y = [startpoints(:,2)'; endpoints(:,2)'];
lines.z = [startpoints(:,3)'; endpoints(:,3)'];

figure;
line(lines.x,lines.y,lines.z);
