Problem with Lid driven cavity collocated grid SIMPLE
 

Hello everyone,
I am trying to write code about 2D Lid driven cavity in collocated grid with SIMPLE algorithm in MATLAB.

I have tried to write code in staggered grid with SIMPLE and it is sucessful.

But in the collocated, I run the code and the values seem not to be convergent.
I discretize the governing equations with FVM and interpolate using the first-order upwind method.
I follow the steps in the book about the collocated grid SIMPLE but it doesn't work.
I also check my code with the flow chart by @Freewill1 Link:https://www.cfd-online.com/Forums/ma...ated-grid.html.
It seems my code is logical correct but it can't run well.

I hope someone can take a look at my code and help me point out my errors.

Thank you!