Forums before death by AOL, social media and spammers... "We can't have nice things"
|    sci.math.symbolic    |    Symbolic algebra discussion    |    10,432 messages    |
[   << oldest   |   < older   |   list   |   newer >   |   newest >>   ]
|    Message 9,125 of 10,432    |
|    Nasser M. Abbasi to All    |
|    question on Kamke differential equation     |
|    30 Jul 16 18:37:33    |
      From: nma@12000.org              I am solving each Kamke ODE by hand, and I am up to problem #22       (it is slow process) which I have question on.              Maple gives the answer to Kamke #22. The problem is              eq:=diff(y(x),x)-y(x)^2-y(x)*sin(2*x)-cos(2*x) = 0;              Maple answer is in terms of Heun Confluent function       and its derivative. (see (2) below). This ODE is       Riccati first order non-linear.              But the solution done by hand, and it matches the one given       in Kamke, only contains tan(x) and one unresolved integral               u:=cos(x)        I:=int( exp(-u^2)/u^2, x)              This integral has no closed form. (or does it?) So       the general solution to the ode is in terms of I. Here it is               y:= tan(x) + exp(-u^2)/(u^2 * (I + C)) (1)              The above is also what Kamke shows in the book.              C is integration constant. Maple answer on the other hand,       is very complicated using HeunC and its derivative. Here it is              restart;       eq:=diff(y(x),x)-y(x)^2-y(x)*sin(2*x)-cos(2*x) = 0;       sol:=dsolve(eq,y(x));              sol := y(x) = (2*HeunCPrime(1, 1/2, -1/2, -1, 7/8, (1/2)*       cos(2*x)+1/2)*_C1*cos(2*x)/((2*cos(2*x)+2)^(1/2)*       (_C1*HeunC(1, 1/2, -1/2, -1, 7/8, (1/2)*cos(2*x)+1/2)*       (2*cos(2*x)+2)^(1/2)+HeunC(1, -1/2, -1/2, -1, 7/8, (1/2)       *cos(2*x)+1/2)))+(HeunCPrime(1, -1/2, -1/2, -1, 7/8, (1/2)       *cos(2*x)+1/2)*(2*cos(2*x)+2)^(1/2)+2*       HeunCPrime(1, 1/2, -1/2, -1, 7/8, (1/2)*cos(2*x)+1/2)*_C1+       2*HeunC(1, 1/2, -1/2, -1, 7/8, (1/2)*       cos(2*x)+1/2)*_C1)/((2*cos(2*x)+2)^(1/2)*       (_C1*HeunC(1, 1/2, -1/2, -1, 7/8, (1/2)*cos(2*x)+1/2)*       (2*cos(2*x)+2)^(1/2)+HeunC(1, -1/2, -1/2, -1, 7/8, (1/2)*cos(2*x       +1/2))))*sin(2*x)               (2)              So it must be then that (1) == (2) ?              Does this mean there is analytical solution to "I" in terms       of special function HeunC? Since I do not see the integral "I" in       Maple solution.              If the solution to the ODE is unique (up to arbitrary constant),       then one must be able to show that (1)<=>(2)?              Is this correct?              Can someone shed some light on this? And does this mean       integral "I" can be expressed in closed form using these HeunC functions?              So how did Maple come up with answer it did if this       integral has no closed form solution?              The answer by Maple must be correct, since Maple odetest()       returns 0.              Here is a link to the problem with detailed solution and Maple answer              http://12000.org/my_notes/kamek/version_3/KEsu22.htm#subsection_22              ps. Mathematica unable to solve this one for some reason.              Can your CAS solve it? What solution does it give?       The book solution or something else?              --Nasser              --- SoupGate-Win32 v1.05        * Origin: you cannot sedate... all the things you hate (1:229/2)    |
[   << oldest   |   < older   |   list   |   newer >   |   newest >>   ]
(c) 1994, bbs@darkrealms.ca