       • To: mathgroup at smc.vnet.net
• From: "Peltio" <peltio at trilight.zone>
• Date: Wed, 20 Apr 2005 05:32:59 -0400 (EDT)
• References: <d3vn1h\$jm9\$1@smc.vnet.net> <d42lb6\$3ki\$1@smc.vnet.net>
• Reply-to: "Peltio" <peltioNOSPAM at despammed.com.invalid>
• Sender: owner-wri-mathgroup at wolfram.com

```"Jim Martin" wrote

> My end goal is to get that
>solution written into a spreadsheet but I have my doubts about any
>spreadsheet (excel or gnumeric) being able to handle the imaginary numbers.

That's unfortunate, because a solution that seems not to suffer from
some sort of discontinuity for positive values of the parameters is that
involving the logartithm. For example this

intF[a_, ha_, hb_, f_, f0_] := With[
{g = Sqrt[hb^2 - 4a^2 f^2 Pi^2]},
(f0 + ha) hb (1/(2 f)
- hb/(2g f Pi)I  Log[(2a f Pi - I g)/(2 a f Pi + I g)])
- 2a ha
]

This is a real quantity, since the log a ratio of complex conjugate numbers
is imaginary and so -I Log[(x+I y)/(x+I y)] happens to be real.

-I Log[(x-I y)/(x+I y)] =
= -I (Log[Abs[x-I y]/Abs[x+I y]] +I Arg[(x-I y)/(x+I y)])
=  -I Log + Arg[(x-I y)/(x+I y)]
= 0 + Arg[(x-I y)/(x+I y)]
= -2 Arg[x+I y]
= -2 ArcTan[y/x]

The problem is that Arg, ArcTan[] and ArcTan[.,.] are affected by different
forms of discontinuity due to the modular world of angles (or the polidromy
of Log, if you prefer). So some of the different equivalent forms that can
be found by expressing log in some other form, are only valid modulo Pi or 2
Pi. Have a look at this graphics, for example

funs = {
Chop[(-I) Log[(x-I y)/(x+I y)]],
Arg[(x-I y)/(x+I y)],
ArcTan[(2 x y)/(y^2-x^2)],
ArcTan[y^2-x^2, 2 x y],
-2 ArcTan[y/x],
- 2 ArcTan[x, y]};
Block[{\$DisplayFunction = Identity},
plots = (Plot[Evaluate[#1 /. x -> 1], {y, -10, 10},
Frame -> True, PlotRange -> {-Pi, 2*Pi},
PlotLabel -> #1] & ) /@ funs
];
Show[GraphicsArray[Partition[plots, 3]]];

The following definition for your integral solution seems to do well for
positive valued parameters, and it could be plotted in a spreadsheet with a
well implemented arctanget function:

intF[a_, ha_, hb_, f_, f0_] := With[
{g = Sqrt[hb^2 - 4a^2 f^2 Pi^2]},
(f0 + ha) hb (1/(2 f) - hb/(g f Pi) ArcTan[ g/( 2 a f Pi)]) - 2a ha
]

Try to plot it into your spreadsheet to see if you experience some kind of
discontinuity that should not be there. Something like this should come up:

Plot[intF[a, 3, 50, 1, 8], {a, 0, 30}, Frame -> True]

cheers,
Peltio