MathGroup Archive 2005

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

Search the Archive

Re: How to use "Solve" to get only real number solutions?

  • To: mathgroup at smc.vnet.net
  • Subject: [mg55135] Re: [mg55132] How to use "Solve" to get only real number solutions?
  • From: Bob Hanlon <hanlonr at cox.net>
  • Date: Mon, 14 Mar 2005 03:43:35 -0500 (EST)
  • Reply-to: hanlonr at cox.net
  • Sender: owner-wri-mathgroup at wolfram.com

intensity[x_] := 1/4+(2*Cos[(Pi*x)/100])/Pi+
      (4*Cos[(Pi*x)/100]^2)/Pi^2-
      (2*Cos[(3*Pi*x)/100])/(3*Pi)-
      (8*Cos[(Pi*x)/100]*Cos[(3*Pi*x)/100])/(3*Pi^2)+
      (4*Cos[(3*Pi*x)/100]^2)/(9*Pi^2);

Off[Solve::ifun];

Select[Chop[Solve[
      intensity[x] == 0.3*(1/4+32/(9*Pi^2)+
              (4*Sqrt[2])/(3*Pi)),x]],
  FreeQ[x/.#,Complex]&]

{{x -> -47.4207},{x -> 47.4207}}


Bob Hanlon

> 
> From: "PengYu.UT at gmail.com" <PengYu.UT at gmail.com>
To: mathgroup at smc.vnet.net
> Date: 2005/03/13 Sun AM 04:57:54 EST
> To: mathgroup at smc.vnet.net
> Subject: [mg55135] [mg55132] How to use "Solve" to get only real number solutions?
> 
> Hi All
> 
> For example I have the following function and I want to find the
> solution of the quation shown below. But Mathematica showns not only
> real number solutions. Do you know how to restrict the solutions to
> only real numbers. Thanks!
> 
> Best wishes,
> Peng
> 
> intensity[x_, δ_, p_, λ_] := 1/4 + (2*Cos[(Pi*x)/100])/Pi +
> (4*Cos[(Pi*x)/100]^2)/Pi^2 - (2*Cos[(3*Pi*x)/100])/(3*Pi) -
>  (8*Cos[(Pi*x)/100]*Cos[(3*Pi*x)/100])/(3*Pi^2) +
> (4*Cos[(3*Pi*x)/100]^2)/(9*Pi^2)
> 
> Solve[intensity[x, 0, 200, 50] == 0.3*(1/4 + 32/(9*Pi^2) +
> (4*Sqrt[2])/(3*Pi)), x]
> 
> 


  • Prev by Date: Re: Mathematica cannot simplify a product of UnitStep functions
  • Next by Date: Re: How to use "Solve" to get only real number solutions?
  • Previous by thread: Re: How to use "Solve" to get only real number solutions?
  • Next by thread: Re: How to use "Solve" to get only real number solutions?