No Title  |POLVER05_0 |1
d(lambda) / d(t) =1
lambda(0) = -1
Pe=1
Da=1
Pe1=5
q=(1+4*Da/Pe)^0.5
A=(1+q)*exp(q*Pe*(1-lambda)/2)
B=(1-q)*exp(-q*Pe*(1-lambda)/2)
A0=(1+q)*exp(q*Pe*(1-0)/2)
B0=(1-q)*exp(-q*Pe*(1-0)/2)
C=(1+q)^2*exp(q*Pe/2)
D=(1-q)^2*exp(-q*Pe/2)
Sai=(2*exp(Pe*lambda/2)*(A-B))/(C-D)
Sai0=(2*exp(Pe*0/2)*(A0-B0))/(C-D)
Sai1=1-(1-Sai0)*exp((lambda*Pe1))
Sa=if (lambda=>0) then (Sai) else (Sai1)
t(0) =0
t(f) = 2