I have attached the file for reference. The first two blocks of code are great. They give me what I need and I have confirmed that by plotting them.
My trouble is with the third block. I am trying to do an iterative procedure and during each iteration I want to compute the value sigmaH, and if (sigmaH*Pm)>0.83, I want to perform the subsequent two operations. At the end of the process I need to find the value of SE.
Right now the code does nothing. I do not know what is going on. Can someone help me/point me in a right direction on how to approach this?
Attachment: iteration.nb, URL: ,