MathGroup Archive 2000

[Date Index] [Thread Index] [Author Index]

Search the Archive

Complex Function Plot

  • To: mathgroup at smc.vnet.net
  • Subject: [mg25587] Complex Function Plot
  • From: Roberto Brambilla <rlbrambilla at cesi.it>
  • Date: Mon, 9 Oct 2000 21:43:27 -0400 (EDT)
  • Sender: owner-wri-mathgroup at wolfram.com

Thank you Bob.
Your solution is OK in almost
the totality of cases. But with a function suggested
by David Park

f[t_]:=BesselK[3, I t]

I've found (with my old Pentium I, 100MHz):

(*trivial solution*)
Plot[{f[t]//Re,f[t]//Im},{t,0.01,20}]//Timing
{32.57 Second}

(*Bob Hanlon solution*)
Plot[Evaluate[{(func=f[t])//Re,func//Im}],{t,0.01,20}]//Timing
{40.04 Second}

(*David Park module, 100 evaluations of f(t) *)
ReImPlot[f, {0.01, 20, 20/100}]//Timing
{10.22 Second}

So it turns out that computing f(t) twice is quicker
than once!

Bye, Roberto
Roberto Brambilla
CESI
Via Rubattino 54
20134 Milano
tel +39.2.2125.5875
fax +39.2.2125.610
rlbrambilla at cesi.it



  • Prev by Date: tracing lines on surface-graphics
  • Next by Date: Harmonic Analysis (Harmonic Matching) (Symbolic) in Mathematica
  • Previous by thread: Re: Complex Function Plot
  • Next by thread: Re: Q: FastBinaryFiles-package: how to detect EndOfFile?