I am very thankful you kind effort. Your file which you sent me I used it to found the solution for H=[0,1000,2000,3000,4000] using your program (attached). Then I calculate the solution using my program (attached) for the same values of H but they are different. My results are correct but I want to calculate the same result using your program. Please let me know where is the problem because solution dosent match. Your program more better than me, and convient to use. I am waiting your positive response.
With my best regards and sincerely.
School of Mathematical Sciences
Peking University, Beijing, China