I want to solve my PDE using Homotopy Perturbation Method. I have error while solving i dont know why it is showing error. If anyone knows the what mistake i have done in the code, kindly correct me. I want to find the equation of W(r,z) . dp/dz is a pressure terms which is independent of r. And the initial guess is taken as w0 := (-eta^2 + r^2)/2.  Help me to solve this.

restart;
PDEtools[declare](w(r, z));
               w(r, z) will now be displayed as w

N := 1;

w0 := (-eta^2 + r^2)/2:
w(r,z) := sum(p^i*w[i](r, z), i = 0 .. N)
                   w[0](r, z) + p w[1](r, z)

HPMEq := (1 - p)*diff(w(r, z), r $ 2) + p*(diff(w(r, z), r $ 2) + diff(w(r, z), r)/r - (1 + lambda)*(dp/dz + A2*M^2*w(r, z) + A1*w(r, z)/Da - m^2*UHS*BesselI(0, m*r)/BesselI(0, m*eta) + A3*sin(Phi))/A1);
                                                                /
                                                                |
        // d  / d            \\     / d  / d            \\\     |
(1 - p) ||--- |--- w[0](r, z)|| + p |--- |--- w[1](r, z)||| + p |
        \\ dr \ dr           //     \ dr \ dr           ///     \

  / d  / d            \\     / d  / d            \\
  |--- |--- w[0](r, z)|| + p |--- |--- w[1](r, z)||
  \ dr \ dr           //     \ dr \ dr           //

     / d            \     / d            \                       
     |--- w[0](r, z)| + p |--- w[1](r, z)|      /             /  
     \ dr           /     \ dr           /   1  |             |dp
   + ------------------------------------- - -- |(1 + lambda) |--
                       r                     A1 \             \dz

         2                            
   + A2 M  (w[0](r, z) + p w[1](r, z))

                                       2                    
     A1 (w[0](r, z) + p w[1](r, z))   m  UHS BesselI(0, m r)
   + ------------------------------ - ----------------------
                   Da                   BesselI(0, m eta)   

                  \
                \\|
                |||
   + A3 sin(Phi)|||
                ///

for i from 0 to N do
    equ[1][i] := coeff(HPMEq, p, i) = 0;
end do;
                     d  / d            \    
                    --- |--- w[0](r, z)| = 0
                     dr \ dr           /    

                           d                                   
                          --- w[0](r, z)      /             /  
 / d  / d            \\    dr              1  |             |dp
 |--- |--- w[1](r, z)|| + -------------- - -- |(1 + lambda) |--
 \ dr \ dr           //         r          A1 \             \dz

                                          2                    
          2              A1 w[0](r, z)   m  UHS BesselI(0, m r)
    + A2 M  w[0](r, z) + ------------- - ----------------------
                              Da           BesselI(0, m eta)   

                 \\    
                 ||    
    + A3 sin(Phi)|| = 0
                 //    

NULL;
cond[1][0] :=  D*w[0](0,z) = 0, w[0](eta,z) = 0;
for j to N do
    cond[1][j] := D*w[j](0,z) = 0,w[j](eta,z) = 0;
end do;
               D w[0](0, z) = 0, w[0](eta, z) = 0

               D w[1](0, z) = 0, w[1](eta, z) = 0

for i from 0 to N do
    dsolve({cond[1][i], equ[1][i]}, {w[i](r,z)});
    w[i](r,z) := rhs(%);
end do;
w(r,z) := evalf(simplify(sum(w[n](r,z), n = 0 .. N)));
convert(w(r,z), 'rational');
w(r,z) := subs(r = 0, w(r,z));
Error, invalid input: rhs expects 1 argument, but received 2
                    w[0](r, z) + w[1](r, z)

                    w[0](r, z) + w[1](r, z)

                    w[0](0, z) + w[1](0, z)