Решение поставленной задачи в программе PLATO IDE

Текст программы

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