What do you think the spice solver is?
Just express your "matrix" in terms  the relationships between node voltages(potential) and branch currents (flow) and let the Analog solver take care of "the matrix" for you.

-- or for SIMPLE cases, you can directly do the calculations..
keeping in mind that there are NO two dimensional Arrays, not matrix manipulation routines.. so you'll have to write a task or function based on how you choose to name the variables  to do that manipulation for you.

-- just LINK the matlab simulation with the verilog simulation.. I know both Modelsim and all the cadence simulators (AMS/Ncsim/spectre/spectreRF & prob Ultrasim and AMSultra have some form of Matlab/simulink interface..

But the analog simulators have already pulled those matrix operations out of the designers control (and headache) so they can concentrate on the Circuit!!
jbd