Текст программы
real L1,L2,L,lambda1,lambda2,k
dimension T(11), A(11),B(11),C(11),D(11),E(11),F(11)
to=20
a1=2.E-5
L1=8.E-1
L2=2.E-1
tau=600
lambda1=70
lambda2=20
alfa1=100
alfa2=100
tcp1=500
tcp2=100
L=L1+L2
Bi1=alfa1*L1/lambda1
Bi2=alfa2*L2/lambda2
T1=to/tcp1
h=L/10
F0n=a1*tau/L**2
deltaF0=F0n/10
F0=0
Tcp=tcp1/tcp1
k=lambda2/lambda1
Tcp12=tcp2/tcp1
do I=1,11
T(I)=T1
END DO
DO i=1,8
A(i)=1/h**2
C(i)=1/h**2
B(i)=(2/h**2)+(1/deltaF0)
END DO
DO i=10,11
A(i)=1/h**2
C(i)=1/h**2
B(i)=(2/h**2)+(1/deltaF0)
END DO
10 F0=F0+deltaF0
DO I=1,11
D(I)=T(I)/deltaF0
END DO
E(1)=1/(1+Bi1*h)
F(1)=h*Bi1*Tcp/(1+Bi1*h)
Do i=2,8
E(i)=A(i)/(B(i)-E(i-1)*C(i))
F(i)=(D(i)+F(i-1)*C(i))/(B(i)-E(i-1)*C(i))
end do
E(9)=k/(1+k-E(8))
F(9)=(F(8))/(1+k-E(8))
Do i=10,11
E(i)=A(i)/(B(i)-E(i-1)*C(i))
F(i)=(D(i)+F(i-1)*C(i))/(B(i)-E(i-1)*C(i))
END DO
T(11)=Bi2*Tcp12/(1-E(10)+Bi2)
DO I=10,1,-1
T(I)=T(I+1)*E(I)+F(I)
END DO
write (*,*) F0, (T(I),I=1,11)
if (F0.Le.F0n) then
go to 10
else
END if
END
Результат
Вывод
График №1
График №2