Path: csiph.com!x330-a1.tempe.blueboxinc.net!usenet.pasdenom.info!news.albasani.net!news2.arglkargh.de!news.musoftware.de!wum.musoftware.de!fu-berlin.de!uni-berlin.de!individual.net!not-for-mail
From: Axel Vogt <&noreply@axelvogt.de>
Newsgroups: comp.soft-sys.math.maple
Subject: Re: Multiple plots with stepped parameters
Date: Sun, 19 Jun 2011 21:38:03 +0200
Lines: 28
Message-ID: <9671hcFjfmU1@mid.individual.net>
References: <SalmonEgg-4EE972.11212019062011@news60.forteinc.com>
Reply-To: &noreply@axelvogt.de
Mime-Version: 1.0
Content-Type: text/plain; charset=ISO-8859-1; format=flowed
Content-Transfer-Encoding: 7bit
X-Trace: individual.net KrSeaysoVD5KQoUR9xtEtwCj17PF9NmhMrG6RX1PxquTwXeCk=
Cancel-Lock: sha1:F7RePDm52v8BiDn4Nc/Tybo7hWU=
User-Agent: Mozilla/5.0 (Windows; U; Windows NT 5.1; en-US; rv:1.9.2.17) Gecko/20110414 Thunderbird/3.1.10
In-Reply-To: <SalmonEgg-4EE972.11212019062011@news60.forteinc.com>
Xref: x330-a1.tempe.blueboxinc.net comp.soft-sys.math.maple:185

On 19.06.2011 20:21, Salmon Egg wrote:
> My example will based upon blackbody radiation. Suppose the emission at
> a frequency nu at various temperatures goes like
>
> J := 1/(exp(nu/T)-1)
>
> How do I plot J against nu for all T in a list like
>
> {0.5, 1.0. 2.0, 5.0, 10.0}?
>
> I think I see how to do that if I listed separate functions in the plot
> command. ie
>
> 1/(exp(nu/0.5)-1), 1/(exp(nu/1.0)-1), ...
>

You may have a scaling problem for showing it, but for example:

L:=[0.5, 1.0. 2.0, 5.0, 10.0]; # use an ordered list, not a set

P:=NULL;
for t in L do        # run through the list
   eval(J, nu=t);     # evaluate expression
   plot(%, T=0 .. 1); # plot for each
   P:=P, %;           # remember it
end do:

plots[display](P);   # print