vmrnd
              Random arrays from the von Mises distribution.
 r = vmrnd (mu, k) returns an array of random angles
 chosen from a von Mises distribution with location parameter mu and
 concentration parameter k on the interval [-pi, pi].  The size of
 r is the common size of mu and k.  A scalar input functions
 as a constant matrix of the same size as the other inputs.  Both parameters
 must be finite real numbers and k > 0, otherwise NaN is returned.
 When called with a single size argument, vmrnd returns a square
 matrix with the dimension specified.  When called with more than one scalar
 argument, the first two arguments are taken as the number of rows and columns
 and any further arguments specify additional matrix dimensions.  The size may
 also be specified with a row vector of dimensions, sz.
Further information about the von Mises distribution can be found at https://en.wikipedia.org/wiki/Von_Mises_distribution
Source Code: vmrnd