We can now plot the violation amount scaled to a quiver plot along each of the sampled directions.
% Normalize the violation amount for plotting
nullVecMap = nullVecMap./max(abs(nullVecMap),[],'all');
% Define the spacetime point to sample
samplePoint = [1,5,5,5];
figure()
hold on
for i = 1:size(vectorField,2)
quiver3(0,0,0, ...
vectorField(2,i), vectorField(3,i), vectorField(4,i), ...
nullVecMap(samplePoint(1),samplePoint(2),samplePoint(3),samplePoint(4),i),'r')
end
title('Normalized Null Violation For Observer Orientations')