I used something similar to your approach, which may be the beginning for you: here is a sample code:
q = [0 1 0 1]; = qnorm (); x = q (1)/sqrt (1-q (4) * q (4))
x =
0
y = q (2)/sqrt (1-q (4) * q (4))
y =
1.0000
z = q (3)/sqrt (1-q (4) * q (4))
z =
0
quiver3 (0, 0, 0, x, y, z)
, , 3D. 3D, , .