1 | initial version |

You can convert the floating point constants to rationals and use Sympy:

```
sage: y = var("y")
sage: a = 123.123
sage: b = 234234.123
sage: solve(QQ(a)^y==QQ(b), y, algorithm="sympy")
[y == (log(234234123) - 3*log(10))/(log(123123) - 3*log(10))]
```

It seems that Sympy can also provide a numerical solution if its domain is restricted:

```
sage: solve(a^y==b, y, algorithm="sympy", domain="real")
[y == 2.56879371128374]
```

Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.