Dave wrote:
. . .

OK. I'm just a bit suspicious of computer programs some times, as
someone will have to choose an algorithm of some sort.

.. . .

I might be able to attack it with a computer
algebra system - maths never was my strongest subject.

I suppose it's natural to be more suspicious of others' work than your
own. I've personally found the opposite to often be more appropriate.

. . .

But I assume you
are talking of something like NEC which breaks antennas into segments.

Yes. You can find a good description of the method of moments in the
second and later editions of Kraus. The fundamental equation can only be
solved numerically, and the method of moments, used by NEC and MININEC,
is an efficient way to do it.

BTW, I'm also looking for an exact formula for input resistance of a
dipole of arbitrary length. I know is 73.13 Ohms when 0.5 wavelengths
long, but I'm not sure exactly how much it varies when the length
changes (I believe it is not a lot).

There is no exact formula for that, either. Calculating an exact
answer requires knowledge of the current distribution, which is a
function of wire diameter. Assuming a sinusoidal distribution gets you
very close for thin dipoles, but it's not exact. You'll find
calculations based on this assumption in just about any antenna text
such as Balanis or Kraus.

Balanis has it, but leaves it as an integral, without simplifying like
he does for the real part. Yet the formuals for hte real and imaginary
parts look very similar. I might be able to attack it with a computer
algebra system - maths never was my strongest subject.

Hallen's integral equation is exact, but it's not a formula, since you
can't plug numbers into one side and get a result on the other. Nor can
it be solved in closed form at all. That's why so much work was done on
approximate solutions and on developing numerical solution methods. Feel
free to write your own program to solve it, but such programs have
existed for decades and have been verified countless times as well as
being highly optimized.

I thought I'd looked in Krauss and not found it, but perhaps it is
there. I think there is a relatively new version of Kraus, but my copy
is quite old.

Getting the resistance is pretty straightforward once you assume the
shape of the current distribution. Assume some arbitrary current at the
feedpoint which, along with the assumed current distribution, gives you
the field strength in any direction. With the impedance of free space,
this directly gives the power density. Integrate the power density over
all space to get the total radiated power. Then you know how much power
is radiated per ampere of current at the feedpoint, from which you can
calculate the feedpoint resistance.

This calculation is done in all editions of Kraus, I'm sure; I have only
the first and second, but I can't imagine it was deleted in later ones.
Be careful when reading Kraus, however. Unlike many authors, he uses a
uniform, rather than triangular, current distribution for his short
elemental dipole examples. This is equivalent to a very short dipole
with huge end hats, not just a plain short dipole. The half wavelength
and other dipoles in his text are conventional.

. . .

Roy Lewallen, W7EL