Люди хелп! Нужно разрешить уравнение log(x) = 1-sx

AnnaSM

Вообщем проблема следующая. Нужно разрешить типовое уравнение как в субже. Желательно аналитмчески, потому что
мне потом с решением еще работать надо будет.
Вот что мне сказал матлаб.
> solve(log(x)-1+s*x)

ans =

lambertw(s*exp(1/s

>> help lambertw
LAMBERTW Lambert's W function.
W = LAMBERTW(X) is the solution to w*exp(w) = x.
W = LAMBERTW(K,X) is the K-th branch of this multi-valued function.
Reference: Robert M. Corless, G. H. Gonnet, D. E. G. Hare,
D. J. Jeffrey, and D. E. Knuth, "On the Lambert W Function",
Advances in Computational Mathematics, volume 5, 1996, pp. 329-359.
Also available from:
http://pineapple.apmaths.uwo.ca/~rmc/papers/LambertW/index.html
Overloaded methods
help sym/lambertw.m
Я что то не могк понять, что мне дальше делать.
Ссылка что то не открывается. Блин...

anestez

Maple дает более обширное описание (что странно, ибо в Matlab полностью встроено Maple-ядро для символьных вычислений). Привести его здесь?
Да, аналитического решения, очевидно, не существует. Поэтому и придумали специальную функцию.

AnnaSM

да напиши, что там?

AnnaSM

А ты какой Maple используешь?

anestez

LambertW - The Lambert W function
Calling Sequence
LambertW(x)
LambertW(k, x)
Parameters
x - algebraic expression
k - algebraic expression, understood to be an integer
Description
The LambertW function satisfies
LambertW(x) * exp(LambertW(x = x .
As the equation y exp(y) = x has an infinite number of solutions y for each (non-zero) value of x, LambertW has an infinite number of branches. Exactly one of these branches is analytic at 0. In Maple this branch is referred to as the principal branch of LambertW, and is denoted by LambertW(x). The other branches all have a branch point at 0, and these branches are denoted in Maple by LambertW(k, x where k is any non-zero integer. (The principal branch can also be referred to as LambertW(0, x.
The principal branch and the pair of branches LambertW(-1, x) and LambertW(1, x) share an order 2 branch point at -exp(-1). The branch cut dividing these branches is the subset of the real line from -infinity to -exp(-1 and the values of the branches of LambertW on this branch cut are assigned using the rule of counter-clockwise continuity around the branch point. This means that LambertW(x) is real-valued for x in the range -exp(-1) .. infinity, while the image of -infinity .. -exp(-1) under LambertW(x) is the curve -y cot(y) + y*I, for y in 0 .. Pi.
Similarly, the branch corresponding to -1, LambertW(-1, x is real-valued on the interval -exp(-1) .. 0, while the image of -infinity .. -exp(-1) under this branch is the curve -y cot(y) + y*I, for y in -Pi .. 0.
For all the branches other than the principal branch, the branch cut dividing them is the negative real axis. The branches are numbered up and down from the real axis (this is very similar to the way the branches of the logarithm are indexed by the multiple of 2*Pi*I which must be subtracted from the imaginary part to recover the principal branch). Again, the values of the branches of LambertW along the branch cut are determined by the rule of counter-clockwise continuity around the branch point at 0. Thus, the image of the negative real axis under the branch LambertW(k, x) is the curve -y cot(y) + y*I, for y in 2*k*Pi .. (2*k + 1)*Pi if k > 0 and y in (2*k + 1)*Pi .. (2*k + 2)*Pi if k < -1. These curves, therefore, bound the ranges of the branches of LambertW, and in each case, the upper boundary of the region is included in the range of the corresponding branch.
The asymptotic behavior of LambertW at complex infinity and at 0 (for the non-principal branches) is given by
LambertW(k,x) ~ log(k,x) - log(log(k,x +
sum(c(m,n)*log(log(k,x^(m+1)/log(k,x)^(m+n+1 0 <= m,n <= infinity)
where log(x) denotes the principal branch of the logarithm, log(k, x) = log(x) + 2*k*Pi*I and the c(m, n) are constants independent of k. The expansion for LambertW(-1, x) is not valid for x tending to 0 along the negative real axis (the effect of the branch point at -exp(-1) must be considered but holds otherwise.
The LambertW function is closely related to the tree generating function T(x) popularized in the analysis of algorithms discipline. When T[n] counts the number of distinct oriented trees with n labelled vertices and T(x) = sum( T[n]*x^n/n!, n=1..infinity then T(x) = -LambertW(-x).
Reference: R.M. Corless, G.H. Gonnet, D.E.G. Hare, D.J. Jeffrey, and D.E. Knuth. "On The Lambert W Function". Advances in Computational Mathematics 5 (1996): 329-359.
Examples
> LambertW(0);
0
> LambertW(-exp(-1;
-1
> LambertW(1.5+2.5*I);
0.9698925542 + 0.5301457549 I
> LambertW(-1,-.1);
-3.577152064
> LambertW(1,-5/2*Pi);
5/2 I Pi
> LambertW(100, Pi+exp(1)*I);
LambertW(100, Pi + exp(1) I)
> evalf(%);
-5.017543738 + 627.4530224 I
alias can be used to shorten the name, if desired
> alias(W=LambertW);
W
> solve(exp(x) > 3*x, x);
RealRange(-infinity, Open(-W(-1/3
RealRange(Open(-W(-1, -1/3 infinity)
> diff(W(x x);
W(x)
------------
(1 + W(x x
> int(cos(W(k,x x);
(1/2 W(k, x) + 1/2) x cos(W(k, x
---------------------------------- + 1/2 x sin(W(k, x
W(k, x)
See Also
inifcns

anestez

Maple 9.01

AnnaSM

дай плиз скачать дистрибутив maple 9

anestez

А в Лориене поискать для начала не судьба?
http://lorien.local/smb/10.0.2.13/linux/maple/MapleSoft.Maple.v9.01-ISO.ShareReactor.iso

AnnaSM

E:\>ping 10.0.2.13 -t
Pinging 10.0.2.13 with 32 bytes of data:
Request timed out.
Request timed out.
Request timed out.
Request timed out.
Request timed out.
Request timed out.
Request timed out.
Request timed out.
Request timed out.
Судьба, только я думал, у тебя есть этот образ и ты в онлайне, что весьма важно
Оставить комментарий
Имя или ник:
Комментарий: