Skip to main content

Timeline for Solving an Integral Numerically

Current License: CC BY-SA 3.0

10 events
when toggle format what by license comment
Feb 24, 2014 at 7:46 comment added Pankaj Sejwal well one can help you if you paste your tried code by editing your question but I would suggest you take a look at @ubpdqn's answer as well. Otherwise put your code here.
Feb 24, 2014 at 6:51 comment added user12553 @Rorschach, I have tried it in 9 and got the same answer as yours. But when I tried plotting the answer for some values of Q,N, k1,k2,x and {t,0,1), I got an empty graph.
Feb 23, 2014 at 13:46 comment added Pankaj Sejwal @RahulNarain : Actually OP has problem with integration, there are many problems other than that like definition of function a[t] etc. But its just a check and actually not working for version 8 but for 9.
Feb 23, 2014 at 12:57 comment added user484 I'm pretty sure that exp(...) should be Exp[...] instead.
Feb 23, 2014 at 8:24 comment added user12553 Ok, thanks a lot @Rorchach. I will check it out in 9
Feb 23, 2014 at 7:27 comment added Pankaj Sejwal mine is version 9. With version 8 I get -0.0056121 dt exp k1 Q x^2 \[Integral]1/( k2 t (N t)^(3/2) - 0.0056121 exp Q x^2) \[DifferentialD]t
Feb 23, 2014 at 7:22 comment added user12553 Please are you using Mathematica 8 or Mathematica 9? cos I am using Mathematica 8 version
Feb 23, 2014 at 7:15 comment added Pankaj Sejwal could be version mismatch or some subtle issue. I am getting what I have written.
Feb 23, 2014 at 7:07 comment added user12553 Thanks a lot. I applied: Integrate[ Simplify[-((dt exp k1 Q x^2)/(32 N [Pi]^(3/2) t (N t)^(3/ 2) (k2 - (exp Q x^2)/(32 N [Pi]^(3/2) t (N t)^(3/2)))))], t] and got the answer ->dt exp k1 Q x^2 [Integral]1/(-32 k2 [Pi]^(3/2) (N t)^(5/2) + exp Q x^2) [DifferentialD]t this doesn't give the same answer as yours, I dont know why
Feb 23, 2014 at 6:27 history answered Pankaj Sejwal CC BY-SA 3.0