实验二 Mathematica程序设计-程序控制结构.docVIP

实验二 Mathematica的程序设计—程序控制结构 使用条件结构定义分段函数：，并绘制在[-2，3]上的图形。 方法一：In[ ]:= f[x_]:= If[x0,x+1,Exp[x]] Plot[f[x],{x,-2,2}] 方法二：In[ ]:= f[x_]:=Which[x0,x+1,x=0,Exp[x]] Plot[f[x],{x,-2,2}] 运行如下： 编写一个显示不超过n的全部素数的函数。 In[ ]:= f[n_]:=Do [If[PrimeQ[i],Print[i] ],{i,2,n}]/;IntegerQ[n]n1 注：这儿/;IntegerQ/n1是对函数中的变量的约定条件。运行如下： In[ ]:= f[3] Out[ ]= 2 3 In[ ]:= f[10] Out[ ]= 2 3 5 7 结果是竖着下来的，不是横着排列。如果要横着，则要用方法二的列表的方式。 方法二： In[ ]:= f[n_]:=(t={};Do[If[PrimeQ[i],AppendTo[t,i]],{i,2,n}];t)/;IntegerQ[n]n1 运行如下： In[ ]:= f[2] Out[ ]= {2} In[ ]:= f[10] Out[ ]= {2,3,5,7} In[ ]:= f[100] Out[ ]= {2,3,5,7,11,13,17,19,23,29,31,37,41,43,47,53,59,61,67,71,73,79,83,89,97} 3、用三种循环函数For,While,Do来计算10！。 答案：For函数： （1）In[ ]:= For[n=1;i=1,i=10,i++,n=n*i];n Out[ ]= 3628800 注意：一定要注意For的格式，一共分四个部分，各个部分用逗号，隔开，不是C语言中的分号！各部分中，如是各有几个语句组成，则要用分号；隔开，形成复合语句。第四个部分循环体部分可与循环变量改变量合为一体，要注意次序。见方法二。 （2）In[ ]:= For[n=1;i=1,i=10,n=n*i;i++];n Out[ ]= 3628800 （3）In[ ]:= For[n=1;i=10,i0,n=n*i;i--];n Out[ ]= 3628800 这是倒着下来的。 （4）In[ ]:= n=1;For[i=1,i=10,i++,n=n*i];n Out[ ]= 3628800 这是初值n放在循环体前（外）的。 While函数： In[ ]:= n=1;i=1;While[i=10,n=n*i;i++];n Out[ ]= 3628800 注意：While循环的内部只有两部分，第一部分是条件；第二部分是循环体，两部分间用逗号，隔开。如果循环体不是一个语句可以完成的，则要用分号隔开，不可用逗号隔开！另外，初值不要放在循环体内，以免引起不必要的错误。 Do函数 In[ ]:= n=1;Do[n=n*i,{i,1,10}];n Out[ ]= 3628800 Do循环的内部也是只有两部分：第一部分是循环体，第二部分是循环次数，用{ }表示的部分。通常是确定的已知次数的，如100次 （{100}），不可以用变量来改变它。 4、一根绳子长32768米，每天截去一半，问多少天长度小于1米？ 要求此题用三种循环函数For,While,Do来编写。 For函数： In[ ]:= For[n=32768;i=0,n=1,i++,n=n/2];i Out[ ]= 16 In[ ]:= For[n=1,32768=,n++];n Out[ ]= 16 While函数 In[ ]:= n=32768;i=0;While[n=1,i++;n=n/2];i Out[ ]= 16 In[ ]:= n=0;While[32768=,n++];n Out[ ]= 16 Do函数 循环体一定要有If函数配合才能判断,因为这是确定次数的循环. In[ ]:= n=32768;i=0;Do[If[n=1,i++];n=n/2, {100}];i Out[ ]= 16 说明：由于Do函数内部只有两部分：循环体与确定次数的控制部分两个部分组成，在循环体内改变的量，不可能反映到控制部分，所以在循环体内就得有IF语句来自行控制。 编写程序，列出9*9的乘法表来。 此题大部人可能这样做： In[ ]:= Table[i*j,{i,1,9},{j,1,9}] Out[ ]= {{1,2,3,4,5,6,7,8,9},{2,4,6,8,10,12,14,16,18},{3,6,9,12,15,18,21,24,27}, {4,8,12,16,20