Add function for generating random orthogonal symplectic matrices

[apkille]

We already do this when generating random symplectic matrices with randsymplectic, we just simply have to add it to the public API. A function name such as randorthosymplectic might do, although it's long and ugly looking.