MATLAB实验-pi的计算 圆周率的近似计算 数学软件与 及数学实验 .ppt

;vpa(pi,101);１.圆周率π的计算历程;实验时期;几何法时期;分析法时期 ; 接着有多种表达式出现。如沃利斯1650年给出：;发现了下面的公式; 1914年，印度数学家Srinivasa Ramanujan　发表了下面的公式：; 1989年，David 和 Gregory Chudnovsky 发表了下面的公式; 1995年，由David Bailey,Peter Borwein 和 Simon Plouffe 共同发表了下面的圆周率计算公式（简称BBP公式）; 1997年，Fabrice Bellard发表了一个比BBP算法更快的公式;;编写下面的程序: n=10; %选择展开式的次数 s=0; digits(22); %定义计算过程中的精度 for k=1:n s=s+4*(-1)^(k+1)/(2*k-1); end vpa(s,20) %定义显示精度为20位;（2）;;;等分区间数n;数值积分简介;输入： x=0:0.1:1; y= exp(-x.^2); trapz(x,y) 输出： ans = 0.746211; 梯形方法;复化梯形方法;举例;输入初值：;抛物线方法;当区间划分为n(n1)等分时;clear; a=0;b=1; f=inline(4/(1+x*x)); t1=(b-a)/2*(f(a)+f(b)); er=1;s1=0; n=1;;s=3.133333,r=3.133333 s=3.141569,r=3.141569 s=3.141593,r=3.141593 s=3.141593,r=3.141593 s=3.141593,r=3.141593 s=3.141593,r=3.141593 s=3.141593,r=3.141593 s=3.141593,r=3.141593 s=3.141593,r=3.141593 s=3.141593,r=3.141593 s=3.141593,r=3.141593 s=3.141593,r=3.141593 s=3.141593,r=3.141593 s=3.141593,r=3.141593 s=3.141593,r=3.141593 s=3.141593,r=3.141593;cs=0 n=500 %随机取点数 for i=1:n a=rand(1,2); if a(1)^2+a(2)^2=1 cs=cs+1 end end 4*cs/n;3.18400000000000 3.10400000000000 3.13866666666667 3.12080000000000 3.14376000000000