On 6/29/2011 6:58 AM, Lars R. wrote:
>
>
> So the program 'maxima' prints after a very short time:
>
>
> (%i1) binomial(1000,402)/(2^1000);
> (%o1) 138375989102164508869118744263820047102478631240107164548096353\
> 437104485769496863239251132205776942412106967230343812796549218291211\
> 502370211070909495252274618853753553269750995170002584212749043301928\
> 705163749324086671038996613530725409460218145152701829040905136436576\
> 746556205530054031875 / 133938575898283415118553131132500226320175601\
> 463191700930468798546293881390617015311649797351961982265949334114694\
> 143353148393160711539255449807219683732185049182097185302887317763432\
> 563279639273474427276913080937294774265842484594489569299325963286432\
> 1399559710817770957553728956578048354650708508672
>
> (I put spaces around the / to make it clearer)
>
> To get a floating point approximation of the quotient:
>
> (%i2) float(binomial(1000,402)/(2^1000));
> (%o2) 1.0331302104275842E-10
>
Cool :) Another Maxima user!
--
~Mike
Post a reply to this message
|