You can use “scatter3” function to plot your trajectory and also pass the additional arguments, the size of the point and the color of the point (which in your case is the acceleration vector). "scatter3" automatically normalizes the 5th input argument values to the current colormap and uses it.
Example:
>> t = 0:0.1:10;
>> x = 5*cos(t);
>> y = 5*sin(t);
>> z = sin(t)+cos(t)+t; % A random trajectory
>> a = rand(1,101) %acceleration values generated randomly for example
>> scatter3(x,y,z,5,a)
>> colormap(‘jet’) %change colormap to required the required one
