MathGroup Archive 2004

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

Search the Archive

LogIntegral^(-1)

  • To: mathgroup at smc.vnet.net
  • Subject: [mg48620] LogIntegral^(-1)
  • From: "Roger L. Bagula" <rlbtftn at netscape.net>
  • Date: Tue, 8 Jun 2004 00:48:11 -0400 (EDT)
  • Sender: owner-wri-mathgroup at wolfram.com

I have an information theory problem:
Suppose that the noise has a constant expectation value with time
so that when you take the differential you get
DN/dt=0
Mathematica for the Shannon capascity as rate of transmission is:
s1[t_]=D[s[t],t]
n[t_]=s[t]/(Exp[s1[t]/w]-1)
n1[t_]=D[n[t],t]
DSolve[FullSimplify[n1[t]]==0,s[t],t]

I get an answer involving:
LogIntegral^(-1)
specifically:

s[t_]=C[2]*(-1+LogIntegral^(-1)[w*(t-C[1])/C[2]])

I'm pretty sure that C[2] represents the average noise.
The problem is I can't get Mathematica to plot this function.
I replaced c2 and c1 with real constants and used a constant w
while trying to plot Re[s[t]]:
I get errors like:
Plot::"plnr":
     "\!\(Re[\(s[t]\)]\) is not a machine-size real number at \!\(t\) = \
\!\(-0.999999916666666699`\)."
Plot::"plnr":
     \:f3b5\!\(Re[s[t]]\) is not a machine-size real number at t = \
-0.918866016854168421`.
Plot::"plnr":
     "\!\(Re[\(s[t]\)]\) is not a machine-size real number at \!\(t\) = \
\!\(-0.830382400281252586`\)."
General::"stop":
     "Further output of \!\(Plot :: \"plnr\"\) will be suppressed during 
this \
calculation."


  • Prev by Date: Best way to store multiple lists of different sizes?
  • Next by Date: Re: Where does Return return to?
  • Previous by thread: Re: Best way to store multiple lists of different sizes?
  • Next by thread: Re: LogIntegral^(-1)