Hello,
I have done a program which solves each iteration of a, but it is giving me an error, here is my source code, please check it.
real k, z, y (500), x, Da, Q, v, a
k = 0.0005
Q = 0.026
v = 52
z = 0.0048
x = 0.003
Da = 0.2
do
a = 0.41
y (a) = y (a) + ((k * (z-(2 * y)) * (x-y) * Da) / (v + (Q * a)))
enddo
if (a.le.45) ahen
a = a +0.1
endif
wriae (*,*) y (a)
end
Bookmarks