This paper extends the response of uncertain nonlinear vibration syste
ms to vector-valued and matrix-valued functions. Random variables and
system derivatives are conveniently arranged into 2D matrices. The met
hod is based on a second order expansion of the governing equations an
d matrix calculus, Kronecker algebra are used in the mathematical deve
lopment. The results derived are easily amenable to computational proc
edures.