I went through the attached file. If you are interested in a numerical integration, then using Matlab or Maple or COMSOL it would be an easy task and I can help you.
Dear Mohammed , thank you very much for your answer. yes I still interested in solving this problem even numerically. i prefer to solve it by FORTRAN, if it is possible, otherwise, no problem to see the result by the softwares you mentioned .
I went through the equations again. I see that the T is an independent variable, and hence, it can be assumed as a value. Then, the corresponding X and Y would be evaluated. Am I right?
If it is ture, then I need the profiles of P and M and other parameters in order to solve the equations for X and Y and send you the numerical results. Please send me message as I may not see your response just when you published, but please let the discussions remains public.
In my searching this problem, I have found the attachment book. The Author says that this kind of equations are ill-posed problem. However, yes I have to provide you the profiles of P and M (at a definite temperature T ) as a unction of z where z is the dimension of the film of thickness L for a film of thickness that appears in the limit of integration. Also, I will provide you the other parameters a,b,c,d, where a and d are functions of T. please see the attachments : the book and his notes.