Mathematica程序举例1

图3.3里的Mathematica程序(版本8.0)

在Mathematica8.0的命令窗口输入如下命令：Cal[n_,i_]=2i*Binomial[n,i];

Print[ScientificForm[Cal[500,8]*1.0,4]];

Print[ScientificForm[Cal[500,12]*1.0,4]]; Print[ScientificForm[Cal[500,16]*1.0,4]]; Print[ScientificForm[Cal[500,24]*1.0,4]]; Print[ScientificForm[Cal[720,8]*1.0,4]];

Print[ScientificForm[Cal[720,12]*1.0,4]]; Print[ScientificForm[Cal[720,16]*1.0,4]]; Print[ScientificForm[Cal[720,24]*1.0,4]];

按Shift+Enter组合键,运行可得：

2.344×1019

1.828×1027

3.75×1034

9.199×1047

4.41×1020

1.514×1029

1.381×1037

6.913×1051

图4.3.2.4.1里的Mathematica程序(版本8.0)

在Mathematica8.0的命令窗口输入如下命令：Cal[m_,n_]=Binomial[m,n];

Fun1[B_,m_,n_,ρ_]=((Cal[m,ρ]*Cal[B-m,n-ρ])/Cal[B,n]); X1=i 826Fun1104,26,26,i Cal i,80.267810.267i 8

; X2=i 1226Fun1104,26,26,i Cal i,120.2671210.267i 12; X3=i 1626Fun1104,26,26,i Cal i,160.2671610.267i 16; X4=i 2426Fun1104,26,26,i Cal i,240.2672410.267i 24; Y1=i 846Fun1248,46,46,i Cal i,80.267810.267i 8; Y2=i 1246Fun1248,46,46,i Cal i,120.2671210.267i 12; Y3=i 1646Fun1248,46,46,i Cal i,160.2671610.267i 16; Y4=i 2446Fun1248,46,46,i Cal i,240.2672410.267i 24; Print[X1 "\n",X2 "\n",X3 "\n",X4 "\n",Y1 "\n",Y2 "\n",Y3 "\n",Y4 "\n"]; 按Shift+Enter 组合键,运行可得：

0.0000951929

3.92256×10-9

5.19076×10-15

7.7419×10-33

0.00104287

5.9339×10-7

3.76588×10-11

3.96973×10-22

图4.3.2.4.2里的Mathematica 程序(版本8.0)

在Mathematica8.0的命令窗口输入如下命令：

Cal[m_,n_]=Binomial[m,n]; Fun1[B_,m_,n_,_]=((Cal[m,]*Cal[B-m,n-])/Cal[B,n]);

Func1[i_,p_]=Fun1[248,46,46,i]*Cal[i,p]*0.267^p*(1-0.267)^(i-p); Fun[p_]:=i p 46Func1i,p

;

Plot[Fun[p],{p,0,46},AxesOrigin

按Shift+Enter 组合键,运行可得下图：