program
implicit none
integer::num,i,j
integer,parameter::ROW=100,COL=100 !矩陣大小
real::U(ROW)
real::Matrix_Bij(ROW,COL)
real::T(ROW) !數值解
real::Y(ROW) !理論解
real::b,a,G,E,pr,n,L,f,&
xi,xj,bij,bij1,bij2,bij3,bij4,&
uyi,PI,a1,a2,a3,a4,Tyy
open (9,file='input.txt')
read(9,*) f,b,num,PI,E,Pr,uyi
open(10,file='output.txt')
G=E/(2*(1+pr))
L=b*f
write(10,*)"f=",f
write(10,*)"num=",num
write(10,*)"E=",E
write(10,*)"Pr=",pr
write(10,*)"G=",G
write(10,*)"b=",b
write(10,*)"L=",L*f
write(10,*)"uyi:",uyi
call slove_Bij
(num,i,j,xi,xj,a,bij,bij1,bij2,bij3,bij4,pr,G,Matrix_Bij,U,uyi,L,b,a1,a2,a3,a4)
call gauss(U,Matrix_Bij,T,num,uyi)
stop
end program
!