Hello,
there is a number of issues which can influence the solution of a dipole antenna:
1) integration must be precise enough to give a stable system of linear equations (which basically means that if you increase the accuracy of integration it no longer influences the results significantly). Each MoM software should have a parameter that would allow you to control how accurately it calculates integrals. If it doesn't, it's not a serious code.
2) Current approximation over the wire should be accurate enough. If your codes don't take of this automatically than you must check it instead. If they use low-order basis, wire must be segmented into lambda/10 segments or smaller. Decrease this value until you get a stable result. If the code is using higher-order basis, than you are able to increase the order of approximation on the wire. I guess your codes are using low-order, although I didn't work with any of them.
3) What is the wire model? There are several wire models on the market and in the academia, but very popular one is thin-wire approximation. This means that the wire is treated as one linear current along its axis. The wire thickness is then separately taken into account... Anyway, make sure that all of your codes "see" the same thickness of that wire. This is particularly the case for the NEC family of the software, if you decrease the wire thickness, the result for the impedance will approach the theoretical 75 ohms...
4) Excitation model is significant for impedance calculation, so you should be careful to check if the excitation models in all of the programs is equivalent.
At the end, my advice to you is to download a free demo of WIPL-D, which is a really serious commercial code ...
www.wipl-d.com -> register and download the demo.
Regards,
drasko