from scipy.optimize import fsolve
from numpy import *

mu = 0.097
val = 0.01

def f(x):
    return x**(-2.)*exp(-mu*x) -val

print fsolve(f, 1.)