Control

program xxx
implicit none
real*8:: t,a,h,z,ca0,f1,f2,f3,b,cain,J,tw,fw
integer::i,n
!real*8,allocatable:: w1(:),w2(:),w3(:)
real::w1(0:1000),w2(0:1000),w3(0:1000)
CHARACTER(10)::RESP

print*, 'cuanto vale ca(t=0.0)';read*,ca0 !CONDICION INICIAL
print*, 'T(t=0)=?';read*,z
print*, 'Tw(t=0)=?';read*,tw
print*,'Fw=';read*,fw
print*, 'dame el valor de a';read*,a
print*, 'dame el valor de b';read*,b
print*, 'cual es tu incremento h';read*,h
cain=0.5 !CONCENTRACION DE ENTRADAn=(b-a)/h

w1(0)=ca0
w2(0)=z !condicion inicial de temperatura del reactor
w3(0)=tw
J=30.0
do i=0,n
1 CONTINUE
2 CONTINUE
t=a+real(i)*hw1(i+1)=w1(i)+h*f1(t,w1(i),w2(i)) !SE ESTA REALIZANDO EL METODO DE EULER
w2(i+1)=w2(i)+h*f2(t,w2(i),w1(i),w3(i))
w3(i+1)=w3(i)+h*f3(t,w3(i),fw,w2(i))print*,t,'MINUTOS',w1(i),w2(i),w3(i)
IF (t/=J) THEN
ELSE IF (t==J) THEN
cain=0.25*cain+cain
print*;PRINT*, 'CAIN=',cain;print*
J=b+1
GO TO 1
ELSE
GO TO 2
END IF
end doOPEN(10,FILE='E.TXT',STATUS='UNKNOWN')
write(10,*) 't= Ca= T= Tw='
DO i=0,n
t=a+real(i)*h
WRITE(10,*) t,w1(i),w2(i),w3(i)
END DOCLOSE(10)
end program xxx


real*8 function f1(t,w1,w2)
real*8 w,t,cain,K,K0,Ea,R
cain=0.5
K0=7.08E10
Ea=30000
R=1.99
K=K0*EXP(-Ea/(R*w2))f1=(40./48.)*(cain-w1)-K*w1
end
real*8 function f2(t,w2,w1,w3)
real*8 t,w1,w2,w3,tin,k,k0,Ea,R
tin=530
K0=7.08E10
Ea=30000
R=1.99
K=K0*EXP(-Ea/(R*w2))f2=(40./48.)*(tin-w2)-30000./(50.*0.75)*k*w1-(150.*250.)/(50.*48.*0.75)*(w2-w3)
end
real*8 function f3(t,w3,fw,w2)
real*8 t,w3,fw
f3=(150.*250.)/(62.3*3.85)*(w2-w3)+(fw/3.85)*(w3-530.)
END