FORTRAN95程序设计实验.doc

FORTRAN95程序设计实验班级：核工11姓名：薛院院学号：2110302026实验一一.实验内容 一个边长为8.5m的边形草地。计算八边形草地的面积并输出。二.问题分析（流程图）开始 输入边数与边长计算三角形面积计算多边形面积结束三.程序!核工程11!薛院院!2013.03.11!计算8边形的面积PROGRAM COLUMATEREAL(4) A !定义边长REAL(8) B,PINTEGER(1) NUMBER !定义边数REAL(8) AREAL !定义面积PRINT*,请输入多边形的边数READ*,NUMBERPRINT*,请输入多边形的边长READ*,AB=A/2/sin(3.14/NUMBER)P=A/2+BAREAL=SQRT(P*(P-A)*(P-B)*2)PRINT*,该多边形的面积为：,AREALEND!输入边数 8 ； 输入边长8.5!输出结果：43.631297090078四.实验结果实验二一. 实验内容计算还清贷款的月算二. 问题分析开始输入RDP计算还款月数将M取整加1，给与M1输出M1结束三. 程序 !核工程11!薛院院!2013.03.11!计算还清贷款的月数!PROGRAM CALUCATE REAL P,D,R,MINTEGER M1PRINT*,请分别输入贷款月利息贷款数以及每月尝还数READ*,P,D,RA=LOG10(1.0)PRINT*,AM=(LOG10(P)-LOG10(P-D*R)/LOG10(1+R)M1=INT(M+0.5)PRINT*,M1END四. 实验结果实验三一.实验内容 用case结构和seclect语句实现税值计算二.问题分析（流程图）开始输入企业产量用case语句计算缴纳的税值输出缴纳的总税值结束三.程序!核工11!薛院院!2110302026!计算应该缴纳的税program calucate_moneyinteger incomereal:income_ratereal:rate_1=0.03,rate_2=0.05,rate_3=0.07,rate_4=0.10,rate_5=0.14,rate_6=0.20print*,请输入收入（以万为单位）：read*,incomeselect case(income)case(0:1000)income_rate=income*rate_1print*,缴纳的税为： ,income_ratecase(1001:2000)income_rate=income*rate_2print*,缴纳的税为： ,income_rate case( 2001:5000)income_rate=income*rate_3print*,缴纳的税为： ,income_rate case( 5001:10000)income_rate=income*rate_4print*,缴纳的税为： ,income_rate case( 10001:50000)income_rate=income*rate_5print*,缴纳的税为： ,income_ratecase defaultincome_rate=income*rate_6print*,缴纳的税为： ,income_rateend selectend四.实验结果实验四一.实验内容 计算 1-12+13-14+(-1)n+11n，其中n满足1+22+32.+n2a二.问题分析（流程图）开始输入用case语句计算缴纳的税值输出缴纳的总税值结束三.程序!班级：核工11!姓名:薛院院!日期：4.8!计算级数之和program mainreal:sum=1,term,sign=1,s=0integer:a,i,n=0print*,请输入一个较大的数（如10000）read*,ado while(senglish_ave)n1=n1+1if(math(i)math_ave)n2=n2+1if(physics(i)physics_ave)n3=n3+1if(averages(i)average)nn=nn+1end dop1=real(n1)/np2=real(n2)/np3=real(n3)/npp=real(nn)/nwrite(*,100)英语,n1,p1*100write(*,100)数学,n2,p2*100write(*,100)物理,n3,p3*100write(*,100)三门课,nn,pp*100100 format(1x,高于,a,平均成绩的学生人数有,i2,人,占总人数比例,f8.2,%)end四.实验结果 实验六一.实验内容 输入4个整数放在4个数组中，数组向左循环位移1,2,3,4个位置，第一个元素移动到最后一个元素，结果放在数组中，输出移除前后的4个整数。二.问题分析（流程图）开始输入4个整数对数组进行处理输出处理后的数组结束三.程序 !班级：核工程11!名字：薛院院!日期：4.2program mainparameter(n=4)integer:x(n)print*,请输入,n,个整数do i=1,nread*,x(i)enddoprint 100,xdo i=1,n t=x(1) do j=2,n x(j-1)=x(j) enddo x(n)=t print 100,xenddo100 format(4i2)end四.实验结果实验七一.实验内容 用牛顿迭代法求解方程x3-2x2+x-1=0在x=1.5附近的值二.问题分析（流程图）开始输入学生人数输入每个人学生的学号，三门成绩计算每个学生的平均成绩计算每门课平均成绩计算3门课平均成绩统计高于每门课平均成绩和三门课平均成绩的人数及比例输出上述计算值结束三.程序 !班级：核工11!姓名：薛院院!学号：2110302026!日期：4.8program mainparameter(max=1000)real x(max)i=1x(1)=1.5 do i=1,1000j=i+1if(abs(x(j)-x(i)10e-7) then x(j)=x(i)-f(x(i)/f1(x(i)endifenddoprint*,计算的该数值为:x=,x(j)end四.实验结果 实验八一.实验内容 对一组数字进行处理二.问题分析（流程图）开始输入一批数对数字进行排序统计奇数个数统计偶数个数输出这批数输出奇偶数个数删除相同的整数统计奇数个数统计偶数个数输出数组及奇偶数结束三.程序 !班级：核工11!姓名：薛院院!学号：2110302026!日期:2013.4.7!主程序program mainparameter (max=50)integer a(max),n,n1,n2call input(a,n)call sort(a,n)n1=odd(a,n)n2=even(a,n)print*,这批整数call output(a,n)print*,这批整数的奇数个数为：,n1print*,这批整数的偶数个数为：,n2call dele(a,n)n1=odd(a,n)n2=even(a,n)print* ,剩余整数call output(a,n)print*,剩余整数奇数个数,n1print*,剩余整数偶数个数,n2containssubroutine input(a,n)integer a(max),nprint(a),请输入这批整数个数n：read*,nprint(a,i2,a),请输入,n,个整数：read*,(a(i),i=1,n)end subroutine!数据输出子程序subroutine output(a,n)integer a(n),nprint*,aend subroutine!整数排序子程序subroutine sort(a,n)integer a(n),ninteger : i, j, m, tdo i =1,n m = ido j = i+1, nif (a(j) a(m) thenm = jend ifend do if (i /= m) then t = a(i) a(i) = a(m) a(m) = t end if enddoend subroutinesubroutine dele(a,n)integer a(n),n,i,ji=1do while(in)j=i+1do while(j=n)if (a(i)=a(j) thendo k=j+1,na(k-1)=a(k)enddon=n-1elsej=j+1endifenddoi=i+1enddoend subroutine!统计奇函数function odd(a,n)integer a(n),n,odd,oddlodd