@tomleslie I am going to compare the coefficiets and used the following command

printlevel := 2; Equation := 4;
for i from 0 by 1 while i <= Equation do 
for j from 0 by 1 while j <= Equation do 
C[i, j] := eval(g(i, j)) 
end do end do

but got the following errer.

g called with arguments: 0, 0
 #(g,1): coeff(collect(algsubs(x^n*y^m = z,select(j -> patmatch(j,a::anything*x^n*y^m),H1)),z),z)
Error, (in algsubs) no variables appearing rationally in pattern
 locals defined as: g = g, pat = 1, rep = z, vars = [], opt = remainder, v = v, bad = bad
Please see the attachement

coeffs.mw

@tomleslie dear thanks for your valuable comments!

@acer @Carl Love @Kitonum @Preben Alsholm

@Carl Love thanks for your answer.

@acer @Carl Love @Preben Alsholm @Kitonum

@tomleslie Thanks for you kind response and fix my problem.

@Joe Riel Thanks for you comments. So how to fix the problem in this case?

@Carl Love Now there is no evalm but face more problem than before.

Help.mw

@Carl Love Please the following file

Help.mw

@Carl Love thanks for your answer! I remove evalm but face same problem.

@acer @Carl Love @Kitonum @Preben Alsholm

@acer @Carl Love @Kitonum @Preben Alsholm

@Carl Love I need the result in trignometry or hyperbolic functions

@Carl Love thanks for your reply now what happend if we integrate the 

int(1/sqrt(u^3+u^2*d[2]+u*d[1]+d[0]), u);

@Carl Love thanks for your answer. I was typing mistake, acutally i want to integrate the following expression

int(1/(a[3]*(u(eta)^3+d[2]*u(eta)^2+d[1]*u(eta)+d[0])), u(eta));