Mathematica 9 is now available
Services & Resources / Wolfram Forums
-----
 /
MathGroup Archive
2003
*January
*February
*March
*April
*May
*June
*July
*August
*September
*October
*November
*December
*Archive Index
*Ask about this page
*Print this page
*Give us feedback
*Sign up for the Wolfram Insider

MathGroup Archive 2003

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

Search the Archive

Re: Handling expressions. Could any one find a more direct way?

  • To: mathgroup at smc.vnet.net
  • Subject: [mg39928] Re: [mg39911] Handling expressions. Could any one find a more direct way?
  • From: Daniel Lichtblau <danl at wolfram.com>
  • Date: Wed, 12 Mar 2003 02:28:30 -0500 (EST)
  • References: <200303110738.CAA08690@smc.vnet.net>
  • Sender: owner-wri-mathgroup at wolfram.com

Guillermo Sanchez wrote:
> 
> (*A need to obtain   {k10, k21, k12, k13, k31}  knowing that A = b
> where*)
> 
> A = 0.1/(0.1 + s) + 0.2/(1 + s) + 0.7/(5 + s);
> 
> b = (k21*k31 + k21*s + k31*s + s^2)/(k10*k21*k31 + k10*k21*s +
> k13*k21*s +
>      k10*k31*s + k12*k31*s + k21*k31*s + k10*s^2 + k12*s^2 + k13*s^2 +
>      k21*s^2 + k31*s^2 + s^3);
> 
> (*Here is my procedure*)
> 
> A1 = Together[ExpandAll[Factor[A]]]
> 
> n1 = Numerator[A1]
> 
> d1 = Denominator[A1]
> 
> n2 = Collect[Numerator[b], s]
> 
> d2 = Collect[Denominator[b], s]
> 
> (*Now, handing the solution for previous output -they are been
> eliminated in this mail- I can obtain  {k10, k21, k12, k13, k31}*)
> 
> Solve[{k21*k31== 0.67, 2.39 == k21 + k31, k10*k21*k31 == 0.5,
>    k10*k21 + k13*k21 + k10*k31 + k12*k31 + k21*k31 == 5.6,
>    k10 + k12 + k13 + k21 + k31 == 6.1}, {k10, k21, k12, k13, k31}]
> 
>   (* Could any one find a more direct way?*)
> 
>   Guillermo Sanchez

SolveAlways can be used to advantage here.

aa = 0.1/(0.1 + s) + 0.2/(1 + s) + 0.7/(5 + s);
bb = (k21*k31 + k21*s + k31*s + s^2)/(k10*k21*k31 + k10*k21*s +
  k13*k21*s + k10*k31*s + k12*k31*s + k21*k31*s + k10*s^2 +
  k12*s^2 + k13*s^2 + k21*s^2 + k31*s^2 + s^3); 

In[27]:= InputForm[SolveAlways[aa==bb, s]]
Out[27]//InputForm= 
{{k10 -> 0.746268656716418, k12 -> 1.2548842674339165, 
  k13 -> 1.7088470758496648, k21 -> 0.32435368834411277, 
  k31 -> 2.065646311655888}, {k10 -> 0.746268656716418, 
  k12 -> 1.7088470758496648, k13 -> 1.2548842674339165, 
  k21 -> 2.065646311655888, k31 -> 0.32435368834411277}}

Note that with inexact arithmetic you rely on the vagaries of
approximate polynomial manipulation. Hence this would probably be more
reliable in general if you use exact arithmetic, e.g. start with

aa = 1/10/(1/10 + s) + 1/5/(1 + s) + 7/10/(5 + s); 


Daniel Lichtblau
Wolfram Research


  • Prev by Date: RE: SetPrecision & Plot bug?
  • Next by Date: Re: EllipticTheta error message
  • Previous by thread: Handling expressions. Could any one find a more direct way?
  • Next by thread: Re: Handling expressions. Could any one find a more direct way?