Groups | Search | Server Info | Keyboard shortcuts | Login | Register [http] [https] [nntp] [nntps]


Groups > comp.soft-sys.math.mathematica > #1906 > unrolled thread

plot of nmaximized result

Started bytarun dutta <tarunduttaz@gmail.com>
First post2011-04-27 09:39 +0000
Last post2011-04-28 10:35 +0000
Articles 2 — 2 participants

Back to article view | Back to comp.soft-sys.math.mathematica


Contents

  plot of nmaximized result tarun dutta <tarunduttaz@gmail.com> - 2011-04-27 09:39 +0000
    Re: plot of nmaximized result "Sjoerd C. de Vries" <sjoerd.c.devries@gmail.com> - 2011-04-28 10:35 +0000

#1906 — plot of nmaximized result

Fromtarun dutta <tarunduttaz@gmail.com>
Date2011-04-27 09:39 +0000
Subjectplot of nmaximized result
Message-ID<ip8o93$pjd$1@smc.vnet.net>
this is the following program.here I want to plot between x VS max.but
I can not get the result.

Do[d = 2635/33859;
 r = 0.6;
 a1 = 0.7;
 a2 = 0.13;
 g = 33859;
 b = 0.6631;
 x = (b*j*(j + 1))/g;
 f = (1/(a1 + a1*z + a2*z^2))*(b/g)*j*(j + 1) + (1 - Exp[-z])^2 +
   d*Exp[-r*z];
 {max, val} = NMaximize[f, z];
 list1 = Print[(*"x= ",*)x];
 list2 = Print[(*" max= ",*)max], {j, 0, 5}]
ListPlot[{list1, list2}]

any help will be appreciated.
regards,
tarun

[toc] | [next] | [standalone]


#1953

From"Sjoerd C. de Vries" <sjoerd.c.devries@gmail.com>
Date2011-04-28 10:35 +0000
Message-ID<ipbft7$ae1$1@smc.vnet.net>
In reply to#1906
As the error message states your function is unbounded, so
maximization is useless. See for yourself by adding Plot[f, {z, -2,
2}] // Print; after the definition of f in your loop

Cheers -- Sjoerd

On Apr 27, 11:39 am, tarun dutta <tarundut...@gmail.com> wrote:
> this is the following program.here I want to plot between x VS max.but
> I can not get the result.
>
> Do[d = 2635/33859;
>  r = 0.6;
>  a1 = 0.7;
>  a2 = 0.13;
>  g = 33859;
>  b = 0.6631;
>  x = (b*j*(j + 1))/g;
>  f = (1/(a1 + a1*z + a2*z^2))*(b/g)*j*(j + 1) + (1 - Exp[-z])^2 +
>    d*Exp[-r*z];
>  {max, val} = NMaximize[f, z];
>  list1 = Print[(*"x= ",*)x];
>  list2 = Print[(*" max= ",*)max], {j, 0, 5}]
> ListPlot[{list1, list2}]
>
> any help will be appreciated.
> regards,
> tarun

[toc] | [prev] | [standalone]


Back to top | Article view | comp.soft-sys.math.mathematica


csiph-web