厦门理工吴荣彬期末作业

发布 2022-09-09 20:22:28 阅读 5422

1、用牛顿迭代法求在区间(-1,1)内sinx=x/2,x值。

编程如下:read(*,x

n=110 x1=x

f=sin(x1)-0.5*x1

f1=cos(x1)-0.5

x=x1-f/f1

write(*,100)n,x1,x

n=n+1if(abs(x-x1).gt.1e-6) goto 10

100 format(1x,'n=',i3,3x,'x1=',f15.7,3x,'x=',f15.7)

end2、建立一个文件(数组从1-1000),调用文件找出素数并建立新文件。

编程如下:parameter (n=1000)

dimension a(n)

data a/1000*0/

common a

real a,l

open(unit=6,file=''status='new',access='sequential',form='formatted')

do 10,i=1,n

a(i)=i

write(6,*)a(i)

10 continue

close(6)

open(unit=7,file=''status='old',access='sequential',form='formatted')

do 20,i=1,n

read(7,*)a(i)

j=sqrt(a(i))

l=230 if((

l=l+1goto 30

endifif((

open(unit=8,file='prime ',status='new',access='sequential',form='formatted')

write(8,*)a(i)

endif20 continue

close(7)

close(8)

end3、用data赋值任意一个数组(三阶矩阵),求这个三阶矩阵的逆矩阵。注意先判断是否有逆矩阵。

编程如下:dimension a(3,3),b(3,3)

real a,b

data a/1,2,3,2,2,1,3,4,3/

write(*,enter 3x3 matrix:'

write(*,100)a

datn=a(1,1)*a(2,2)*a(3,3)+a(2,1)*a(3,2)*a(1,3)+a(3,1)*a(1,2)*a(2,3)

datm=a(3,1)*a(2,2)*a(1,3)+a(2,1)*a(3,3)*a(1,2)+a(1,1)*a(3,2)*a(2,3)

dat=datn-datm

if('it''s irreversible'

if(b(1,1)=(a(2,2)*a(3,3)-a(3,2)*a(2,3))/dat

b(1,2)=-a(1,2)*a(3,3)-a(3,2)*a(1,3))/dat

b(1,3)=(a(1,2)*a(2,3)-a(2,2)*a(1,3))/dat

b(2,1)=-a(2,1)*a(3,3)-a(3,1)*a(2,3))/dat

b(2,2)=(a(1,1)*a(3,3)-a(3,1)*a(1,3))/dat

b(2,3)=-a(1,1)*a(2,3)-a(2,1)*a(1,3))/dat

b(3,1)=(a(2,1)*a(3,2)-a(3,1)*a(2,2))/dat

b(3,2)=-a(1,1)*a(3,2)-a(3,1)*a(1,2))/dat

b(3,3)=(a(2,2)*a(1,1)-a(1,2)*a(2,1))/dat

write(*,dat=',dat

write(*,200)

write(*,output 3x3 inverse matrix:'

write(*,100)b

endif100 format(1x,3f7.3)

200 format(30('-

end4、用子例行设计一个程序,作用是计算惯性矩i,主程序直接套用子程序计算在纯弯曲条件最大的应力(其中弯矩m,圆截面半径r或者弯矩m,矩形截面长宽a,b的值任定)。先判断截面是圆截面还是矩形截面或者其他,如果是圆截面,输入半径r,如果是矩形截面,输入长宽a,b,其他输出字符’i can’t do it 且程序结束。提示:

可以有两个子程序。

编程如下:program main

real r,m,i

write(*,please choose the shape of section:'

write(*,cir=0 rec=1 other=2)'

read(*,n

if(write(*,please input r,m:'

read(*,r,m

call cir(r,i)

xigema=r*m/i

write(*,r=',r,'m','m=',m,''

write(*,100)xigema

elseif(

write(*,please input a,b,m:'

read(*,a,b,m

call rec(a,b,i)

xigema=2*m/(i*b)

write(*,a=',a,'m','b=',b,'m','m=',m,''

write(*,100)xigema

elseprint*,'i can''t do it'

endif100 format(1x,'xigema=',f15.9,2x,'pa')

endsubroutine cir(x,z)

parameter(pi=3.1415926)

z=pi*x**4/4

endsubroutine rec(x,y,z)

z=x*y**3/12end

厦门理工附中作业 一 答案

厦门理工附中2009级高一暑假作业 一 答案。高一数学答案。一 选择题 每小题5分,共60分 二 填空题 每小题5分,共20分 三 解答题 共6小题,共70分,每题应写出适当的解答步骤 17.解 即 0,解得 5分。当时,为a中元素 当时,当时,这样的实数x存在,是或10分。另法 0且。或10分。1...

阅读吴满荣《火烧云》四年级上册

校本研修过程。火烧云 教学设计方案。黔东南州地区锦屏县高柳小学学校学员姓名 吴满荣。设计意图 这个板书直观简洁的反应了课文内容,凸显了重点,符合学生认知特点。教学反思 火烧云 是一篇优美的写景佳作,描写的是火烧云从上来到下去的过程中,颜色和形状的变化特点,表现了作者热爱大自然 热爱生活的思想感情。本...

2023年度锦荣项目监理工作总结

安全方面 由于整个皮带廊工程大部分处于荒山野岭,不安全因素较多,线路长 施工面广,对于监理履行安全责任有很大难度,因此监理部始终将工程安全放在第一位。前期要求总包方在每个标段配置一名专职安全员,尤其是隧道爆破 开挖工程更是危险系数较大,监理人员也是重点监督。首先监督总包 分包的各项安全施工方案和措施...