I forget to attach the file please see now
Code_paper_1.mw

Dear @Carl Love thanks for your answer. I need the inverse laplace of F(q, x, alpha) that is invlaplace(F(q, x, alpha),q,t) then I want to plot it for fix values of t and alpha for x=0..1. 

Dear @vv; thanks for your precious time and valuable suggestions. 

I know each and everything. If you can solve my problem then reply. I am not free to see your bloody answer.

@Preben Alsholm 

Thanks for your kind answer. Now my one more question is which method (buiton command) is best to solve system of nonlinear ODEs numerically. Beacuse some time I face the error like this mehtod is not converges.

@tomleslie 

Dear Sir! I want to see the values of F(eta), G(eta) and H(eta) form 0 to 1 with 0.1 difference. I need the accuracy till 13 decimal places please correct the program accordingly, Thanks for your psoitive response.

@Carl Love 

Really thanks alot for your kind efforts. But still there is problem when I compare the coefficient having "x power 0" (used the command remove(has,f, x)). Please solve the problem the problem is attached.

Equate.mw

@tomleslie 

Dear friend!

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. 

MHD_Flow_H=4000_R=100_alpha=5.mw

MHD_Flow_H_varies_R=100_alpha=5.mw

With my best regards and sincerely.

Muhammad Usman

School of Mathematical Sciences 
Peking University, Beijing, China

@tomleslie 

Thanks dear your reply really apperciate. Actually the orginal file are attached please fix it I will be really thankful to you.

Help.mw

@Carl Love 

Dear Sir I want to plot the above graph from -1..1 where the change will happen? I am waiting your positive response. 

@John Fredsted I required a sepcial lower triangular matrix as shown in attached fine. I know the built-in command which you sent me. If you fix my problem, i will be grateful to you. Thanks

@AmusingYeti 

First I am very thankful to you about your valuable comments. Now I am facing another problem in integration please see the attachment and fixed it as soon as possible. I am waiting your kind comments. Take crae.

Integration.mw

@tomleslie 

Dear how you get this solution, which command you used? Please correct the attached file and send me back. I am waiting your kind response.

@Markiyan Hirnyk 

I do not understand your question what you want to say.

@tomleslie