Student Support Forum: 'NSolve / Solve Problem with Integral' topicStudent Support Forum > General > Archives > "NSolve / Solve Problem with Integral"

 < Previous Comment | Next Comment > Help | Reply To Comment | Reply To Topic
 Author Comment/Response B 06/02/13 1:28pm In[1]:= kg=1; m=1; sek=1; amp=1; (* Dimensions *) c = 299792458 m/sek (* Light Speed *); Gyr = 10^9*365.25*24*3600*sek (* Billion Years *); Glyr = Gyr*c (* Billion Lightyears *); Mpc = 3.085677581*10^22 m (* Megaparsec *); H0 = 67110 m/Mpc/sek ; \[CapitalOmega]R = 4.165*^-5; \[CapitalOmega]M = 0.3175; \[CapitalOmega]\[CapitalLambda] = 0.6825 - \[CapitalOmega]R; \[CapitalOmega]\[CurlyEpsilon] = 0.04; \[CapitalOmega]T = \[CapitalOmega]R + \[CapitalOmega]M + \[CapitalOmega]\[CapitalLambda]; \[CapitalOmega]K = 1 - \[CapitalOmega]T; f[A_Real] := NIntegrate[1/(a Sqrt[\[CapitalOmega]R a^-4 + \[CapitalOmega]M a^-3 + \[CapitalOmega]K a^-2 + \[CapitalOmega]\[CapitalLambda]]), {a, 0, A}, Method -> {"GlobalAdaptive", "MaxErrorIncreases" -> 10000, Method -> "GaussKronrodRule"}, MaxRecursion -> 10000]; Plot[{f[A], 0.03733}, {A, .03, .04}] Out[6]= ...PlotSnipped... In[7]:= Plot[{f[A], 0.03733}, {A, .03, .12}] Out[7]= ...PlotSnipped... In[8]:= FindRoot[f[A] - 0.03733, {A, .1}] Out[8]= {A -> 0.0999987} URL: ,

 Subject (listing for 'NSolve / Solve Problem with Integral') Author Date Posted NSolve / Solve Problem with Integral Yukterez 06/01/13 6:57pm Re: NSolve / Solve Problem with Integral B 06/02/13 1:28pm Re: NSolve / Solve Problem with Integral Yukterez 06/03/13 02:44am Re: Re: NSolve / Solve Problem with Integral Bill Simpson 06/03/13 2:51pm Re: Re: Re: NSolve / Solve Problem with Integral Yukterez 06/04/13 8:55pm Thanks Problem solved Yukterez 06/04/13 10:34pm Re: NSolve / Solve Problem with Integral Yukterez 06/04/13 9:45pm
 < Previous Comment | Next Comment > Help | Reply To Comment | Reply To Topic