TransWikia.com

Finding root of symbolic function

Mathematica Asked by Greenteamaniac on December 7, 2020

I’m trying to find a root of the following equation
$$theta^frac{t}{s}-theta-sqrt{1-p},$$
where $0<tleq1leq s$ and $0<p<1$.
I’ve tried FindRoot and Reduce, but somehow they don’t process the output.

Using reduce function, what I tried is

Reduce[theta^(t/s) - theta - Sqrt[1 - p] == 0 && theta > 0 && theta < 1 &&0 < t<=  1<=s && 0 < p < 1, theta]

Any suggestions on how to solve for a symbolic function?

One Answer

The t/s is somehow a nasty in this question. The inequality chain make t/s limited by 0 and 1. 1 means t==s. At this point the equation reduces to Sqrt[1-p] and is satisfied only for p==1.

Following this path replace u=t/s and solve the equation.

p[theta_, u_] := 1 - (theta^u - theta)^2

Plot3D[p[theta, u], {theta, 0, 1}, {u, 0, 1}, AxesLabel -> Automatic]

enter image description here enter image description here

The quotient t/s makes the equation unnecessary complicated. u makes it trivial.

It is too meaningful to replace the Sqrt[1-p] by q for example and

Plot3D[theta^u - theta, {u, 0, 1}, {theta, 0, 1}]

Plot3D

For very large s the quotient t/s->0 and the equation gets

1-theta-Sqrt[1-p]

Reduce[1 - theta - Sqrt[1 - p] == 0, {p, theta}]

(* theta == 1 - Sqrt[1 - p] *)

Plot3D[theta == 1 - Sqrt[1 - p], {p, 0, 1}, {theta, 0, 1}, 
 AxesLabel -> Automatic, PlotRange -> All]

Plot3D

Reduce[theta^u - theta - Sqrt[1 - p] == 0 && theta > 0 && theta < 1 &&
   0 < p < 1, u]
(* Element[C[1], Integers] && 0 < p < 1 && 0 < theta < 1 && 
  u == (2*I*Pi*C[1] + Log[Sqrt[1 - p] + theta])/Log[theta] *)

Plot3D[Log[Sqrt[1 - p] + theta]/Log[theta], {p, 0, .95}, {theta, 0, .9}, 
 AxesLabel -> Automatic, PlotRange -> All]

Plot3D

This is not so informativ and comfortable to look at if p->1 and theta->1.

Answered by Steffen Jaeschke on December 7, 2020

Add your own answers!

Ask a Question

Get help from others!

© 2024 TransWikia.com. All rights reserved. Sites we Love: PCI Database, UKBizDB, Menu Kuliner, Sharing RPP