This can be done in EMCoS MoM solver. There you can extract in ASCII files all Z, J and E (named Z, I and U).
To do so include following lines in the input file:
Z = yes
U = yes
I = yes
in the display block.
Here what is stated in User Manual:
Commands Z, U and I control producing and filling auxiliary files with additional
data information: ‘InpFile.Z.out’ (impedance matrix), ‘InpFile.U.out’ (excitation,
or right part, column), and ‘InpFile.I.out’ (currents solution), where InpFile is the
name of input file, but without extension .in. This may be done assigning yes to
the corresponding command. If the command is assigned no or none, no
auxiliary file for this command is produced. Currently, commands Z, U and I are
available only for the solvers AMS and LAPACK, AGGR.
The auxiliary files ‘InpFile.Z.out’, ‘InpFile.U.out’ and ‘InpFile.I.out’ are structured
as follows. First line of each file gives, in series, a total number of unknowns,
number of non-boundary edges, number of wire/triangle junctions, number of
wire/wire junctions, and number of ground attachments (for wire segments).
Next lines give, elementwise, the elements of corresponding matrices (first, the
rows, and then columns).
EMCoS Antenna Virtual Lab - there is free Student Version on request from the site.