Find the value x

15 visualizaciones (últimos 30 días)
Pichaya Thochampa
Pichaya Thochampa el 17 de Ag. de 2020
Editada: Matt J el 19 de Ag. de 2020
y = x ^ e - e ^ x is [0 5] when x^e=e^x

Respuesta aceptada

Matt J
Matt J el 17 de Ag. de 2020
Editada: Matt J el 17 de Ag. de 2020
>> [x,fval]=fminbnd(@(x) abs(exp(x)-x^exp(1)),0,5)
x =
2.7183
fval =
4.4743e-11
  4 comentarios
Pichaya Thochampa
Pichaya Thochampa el 19 de Ag. de 2020
Editada: Bruno Luong el 19 de Ag. de 2020
Excuse me, ask me why. When issuing a command of a program that is
>> x = fminbnd (@ (x) (x ^ (exp (1)) - (exp (1) ^ x)), 0,5)
to get the answer 1.0000, but when using
>> x = fminbnd (@ (x) ((exp (1) ^ x) -x ^ exp (1)), 0,5)
is only 2.7183.
Matt J
Matt J el 19 de Ag. de 2020
Editada: Matt J el 19 de Ag. de 2020
Because the first version has a local minimum at x=1. (Note: fminbnd may find a local, rather than a global minimum).
fplot(@ (x) (x ^ (exp (1)) - (exp (1) ^ x)), [0,5]);
whereas the second version has a minimum at x=exp(1).
fplot(@ (x) ((exp (1) ^ x) -x ^ exp (1)), [0,5])

Iniciar sesión para comentar.

Más respuestas (3)

hosein Javan
hosein Javan el 17 de Ag. de 2020
Editada: hosein Javan el 17 de Ag. de 2020
you can rewrite the equation in the form of "x^(1/x) = exp(1/exp(1))". this means that if "f(z)=z^(1/z)" then we have "f(x) = f(exp(1))". giving an inverse function, then it is "x=exp(1)".
syms x
f(x) = x^(1/x) - exp(1/exp(1));
x_ex = solve(f(x),x) % exact value
x_app = vpa(x_ex) % approximation
err_f = subs(f(x),x,x_app) % error
x_ex =
lambertw(0, 51*log(2) - log(3253102820258857))/(51*log(2) - log(3253102820258857))
x_app =
2.7182817878772095927811977987387
err_x =
2.9387358770557187699218413430556e-39
x = 0:0.001:5;
y = x.^exp(1) - exp(x);
plot(x,y)
  2 comentarios
Pichaya Thochampa
Pichaya Thochampa el 17 de Ag. de 2020
Thank you so much:)
hosein Javan
hosein Javan el 17 de Ag. de 2020
don't mention it.

Iniciar sesión para comentar.


KSSV
KSSV el 17 de Ag. de 2020
The value of x would be e i.e exp(1). Check the below:
x = linspace(0,5,10^6) ; % as you said x lies in [0 5]
f = x.^exp(1)-exp(x) ;
[val,idx] = min(abs(f)) ; % pick the value for which f is minimum i.e f = 0
x0 = x(idx) ; % the value of x
f = x0^exp(1)-exp(x0) % check
  1 comentario
Pichaya Thochampa
Pichaya Thochampa el 17 de Ag. de 2020
Thank you so much:)

Iniciar sesión para comentar.


Matt J
Matt J el 17 de Ag. de 2020
x=exp(1)

Categorías

Más información sobre Time Series en Help Center y File Exchange.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by