计算材料学考试..docxVIP

分别用前向 2 点,后向 2 点 ,对称 3 点,对称 5 点方法计算 f(x=lnx 在 x=2 处的一阶导数和误差 ,给出不同 h 值(0.5, 0.2, 0.1, 0.05, 0.02, 0.01, 0.005, 0.002, 0.001, 0.0005, 0.0002, 0.0001下各方法的误差 1.1 对称三点法 program lx implicit none doubleprecision,parameter :: x=2.0d0 doubleprecision :: exact,h,f_sym3,diff_sym3 h=0.5d0 exact=x**(-1 do while (h0 write(*,* Enter value of h (=0 to stop read(*,* h if (h=0.0 then stop h=0 endif f_sym3=(dlog(X+h-dlog(x-h/(2*h diff_sym3=exact-f_sym3 write(*,* h=,h,F_SYM3=,f_sym3,E_SYM3=,diff_sym3 enddo end program lx 1.2 前向两点法 program lx implicit none doubleprecision,parameter :: x=2.0d0 doubleprecision :: exact,h,f_sym3,diff_sym3 h=0.5d0 exact=x**(-1 do while (h0 write(*,* Enter value of h (=0 to stop read(*,* h if (h=0.0 then stop h=0 endif f_sym3=(dlog(x+h-dlog(x/(h diff_sym3=exact-f_sym3 write(*,* h=,h,F_SYM3=,f_sym3,E_SYM3=,diff_sym3 enddo end program lx 1.3 后向两点法 program lx implicit none doubleprecision,parameter :: x=2.0d0 doubleprecision :: exact,h,f_sym3,diff_sym3 h=0.5d0 exact=x**(-1 do while (h0 write(*,* Enter value of h (=0 to stop read(*,* h if (h=0.0 then stop h=0 endif f_sym3=(dlog(x-dlog(x-h/(h diff_sym3=exact-f_sym3 write(*,* h=,h,F_SYM3=,f_sym3,E_SYM3=,diff_sym3 enddo end program lx 1.4 对称五点法 program lx implicit none doubleprecision,parameter :: x=2.0d0 doubleprecision :: exact,h,f_sym5,diff_sym5 h=0.5d0 exact=x**(-1 do while (h0 write(*,* Enter value of h (=0 to stop read(*,* h if (h=0.0d0 then stop h=0 endif f_sym5=(dlog(x-2*h-8*dlog(x-h+8*dlog(x+h-dlog(x+2 *h/(12*h diff_sym5=exact-f_sym5 write(*,(5(A,F8.6,2X h=,h,F_SYM5=,f_sym5,E_SYM5=,diff_sym5 enddo end program lx 分别用梯形法 ,Simpson法,Simpson3/8法,Bode 法计算 2 0cos3d x x π ? 给出 8 个不同 N 值下各方法的误差 3.1 梯形法 doubleprecision function func(x implicit none doubleprecision :: x func=cos(3*x endfunction func program integral implicit none integer :: i,n=1 doubleprecision :: pi,limit_down,limit_up doubleprecision :: exact,h,integ,summ,diff,x doubleprecision :: func exact=-1.0/3.0 pi=acos(-1.0d0 limit_down=0.0d0 limit_up=pi/2.