Getting Euler Angles from a Quaternion in Matlab

Question

I have the X,Y,Z and W components of a Quaternion through time, in 4 separated vectors.

  QW             1x346             2768  double       
  QX             1x346             2768  double
  QY             1x346             2768  double
  QZ             1x346             2768  double              

I want to convert to euler angles in order to plot the 3 euler components in three different subplots through time, so I need to have 3 vectors like these.

  heading(t)          1x346 
  attitude(t)         1x346 
  bank(t)             1x346 

Is there an inmediate way in matlab to get this (I mean a function with input my 4 vectors and output the above 3 vectors) or do I have to write some code to make the conversion for each timestep? Thanks in advance.