Multivariate fractional phase–type distributions