Durbin's elliptic relaxation method!
 

hello.i am using Durbin's elliptic relaxation method with the full RSM to calculate the channel flow. i wrote 1 D program with finite difference method.but the resluts diverged.

i think there is something to do with the boundary conditons. could you give me some suggestions about this method? thank you very much!