Problem recreating example of normal distribution from TI-Nspire [closed]

asked 2017-04-27 19:13:59 +0100

Ross1856
23 ●1 ●5

I have to recreate this formula for calculating normal distribution in Sage:

http ://ntaj.dreamhosters.com/ti.jpg

The picture is from TI-Nspire CAS.

What I'm doing is:

sage: f=1/(0.1sqrt(2pi)) * e^((-1/2)*((x-4)/0.1))^2

sage: float(integral(f, 3.9, 4.1))

0.8696735925295498

I can't get the expected result 0.682689 that Nspire gives. Can anyone see why?

edit retag flag offensive reopen merge delete

Closed for the following reason the question is answered, right answer was accepted by Ross1856
close date 2017-04-28 11:24:30.118994

add a comment

0

answered 2017-04-27 22:04:14 +0100

eric_g
8494 ●10 ●62 ●168 https://luth.obspm.fr/...

You did not get the expected result because the last ^2 is ill placed; it should be inside the parentheses:

sage: f = 1/(0.1*sqrt(2*pi)) * e^((-1/2)*((x-4)/0.1)^2)

Then

sage: numerical_integral(f, 3.9, 4.1)
(0.6826894921370853, 7.579375928402468e-15)

(the second number is an evaluation of the numerical error).

edit flag offensive delete link

add a comment

0

answered 2017-04-27 21:58:16 +0100

dan_fulea
5537 ●5 ●43 ●95

Pari/gp delivers:

? s = 0.1;

? intnum( x=3.9, 4.1, exp( -(x-4)^2/2/s^2 ) / sqrt(2*Pi*s^2) )                                                                
%2 = 0.6826894921370858971704650911

We can do the same in sage, i.e. call gp:

sage: gp( "intnum( x=3.9, 4.1, exp( -(x-4)^2/2/%s^2 ) / sqrt(2*Pi)/%s)" % ( 0.1, 0.1 ) )
0.6826894921370858971704650911

Or call the numerical integral from sage:

sage: integral_numerical( lambda x :    exp( -(x-4)^2/2/0.1^2 ) / sqrt(2*pi)/0.1, 3.9, 4.1 )
(0.6826894921370853, 7.579375928402468e-15)

which computes the approximative value 0.6826894921370853 with an error less than $7,6\cdot 10^{-15}$, as given in the second argument.

edit flag offensive delete link

Comments

I'll have to read up on that. Thank you for the reference.

Ross1856 ( 2017-04-28 11:28:03 +0100 )edit

add a comment

Problem recreating example of normal distribution from TI-Nspire [closed]

Closed for the following reason the question is answered, right answer was accepted by Ross1856
close date 2017-04-28 11:24:30.118994

2 Answers

Comments

Question Tools

Stats

Problem recreating example of normal distribution from TI-Nspire [closed] edit

Closed for the following reason the question is answered, right answer was accepted by Ross1856 close date 2017-04-28 11:24:30.118994

2 Answers

Comments

Question Tools

Stats

Problem recreating example of normal distribution from TI-Nspire [closed]

Closed for the following reason the question is answered, right answer was accepted by Ross1856
close date 2017-04-28 11:24:30.118994