Computation for Formulation 14.16 on page 153

This is the only one of my notebooks created with Version 3.0 that Version 5.1 fails to evaluate correctly.  Namely, the  Coefficient  function in Version 5.1 needs a better implementation.

The version number of Mathematica used to evaluate this notebook:

$VersionNumber

From (14.19), we have

aAA[
        zet_]  =  f'[zet]^2*
          a[f[zet]] + (f'''[zet]/(2*f'[zet])) - (3/
              4) (f''[zet]/f'[zet])^2

The preceding represents  \(TraditionalForm\`\(\(a\^ ** \)\_\(m, 0, 2\)\)(ζ)\)  and we have chosen  A  to suggest asterisk.

We write

PAA[zet_]  =  f'[zet]^p*P[f[zet]]

for (14.20).  In particular, we have

D[PAA[zet],  zet]

for (14.21) and

Collect[ 
      D[PAA[zet], {zet, 2}],  {P''[f[zet]],  P'[f[zet]], 
        P[f[zet]]}, Simplify ]

for (14.22).  Similarly, we introduce

QAA[zet_]  =  f'[zet]^q*Q[f[zet]]

We introduce

v[0]  =  1

v[1]  =  v1

v[2]  =  v2

v[3]  =  v3

v[4]  =  v4

v[5]  =  v5

v[6]  =  v6

v[7]  =  v7

v[8]  =  v8

v[9]  =  v9

v[10]  =  v10

v[11]  =  v11

v[12]  =  v12

v[13]  =  v13

v[14]  =  v14

v[15]  =  v15

v[16]  =  v16

v[17]  =  v17

v[18]  =  v18

v[19]  =  v19

v[20]  =  v20

v[21]  =  v21

v[22]  =  v22

v[23]  =  v23

v[24]  =  v24

v[25]  =  v25

v[26]  =  v26

v[27]  =  v27

v[28]  =  v28

v[29]  =  v29

so that (14.48) yields

\(\(VAA[
          zet_] = \((Sum[
                v[i]*D[PAA[zet], {zet, i}]*D[QAA[zet], {zet, 6 - i}], {i, 0, 
                  6}] + aAA[zet]*
                Sum[\ v[7 + i]*D[PAA[zet], {zet, i}]*
                    D[QAA[zet], {zet, 4 - i}], {i, 0, 4}] + 
              D[aAA[zet], zet]*
                Sum[\ v[12 + i]*D[PAA[zet], {zet, i}]*
                    D[QAA[zet], {zet, 3 - i}], {i, 0, 3}]
 + 
              D[aAA[zet], {zet, 2}]*
                Sum[v[16 + i]*D[PAA[zet], {zet, i}]*
                    D[QAA[zet], {zet, 2 - i}], {i, 0, 2}] + \((aAA[zet])\)^2*
                Sum[\ v[19 + i]*D[PAA[zet], {zet, i}]*
                    D[QAA[zet], {zet, 2 - i}], {i, 0, 2}] + 
              D[aAA[zet], {zet, 3}]*
                Sum[\ v[22 + i]*D[PAA[zet], {zet, i}]*
                    D[QAA[zet], {zet, 1 - i}], {i, 0, 1}] + 
              aAA[zet]*D[aAA[zet], zet]*
                Sum[\ v[24 + i]*D[PAA[zet], {zet, i}]*
                    D[QAA[zet], {zet, 1 - i}], {i, 0, 1}] + 
              v[26]*D[aAA[zet], {zet, 4}]*\ PAA[zet]*QAA[zet] + 
              v[27]*aAA[zet]*D[aAA[zet], {zet, 2}]*\ PAA[zet]*QAA[zet] + 
              v[28]*D[aAA[zet], zet]^2*PAA[zet]*QAA[zet] + 
              v[29]*aAA[zet]^3*\ PAA[zet]*QAA[zet])\) // Expand;\)\)

Length[VAA[zet]]

and

V[
          z_] = (Sum[v[i]*D[P[z], {z, i}]*D[Q[z], {z, 6 - i}], {i, 0, 6}] + 
              a[z]*Sum[ 
                  v[7 + i]*D[P[z], {z, i}]*D[Q[z], {z, 4 - i}], {i, 0, 4}] + 
              D[a[z], z]*
                Sum[ v[12 + i]*D[P[z], {z, i}]*D[Q[z], {z, 3 - i}], {i, 0, 
                    3}] + D[a[z], {z, 2}]*
                Sum[ v[16 + i]*D[P[z], {z, i}]*D[Q[z], {z, 2 - i}], {i, 0, 
                    2}] + (a[z])^2*
                Sum[ v[19 + i]*D[P[z], {z, i}]*D[Q[z], {z, 2 - i}], {i, 0, 
                    2}] + D[a[z], {z, 3}]*
                Sum[ v[22 + i]*D[P[z], {z, i}]*D[Q[z], {z, 1 - i}], {i, 0, 
                    1}] + a[z]*D[a[z], z]*
                Sum[ v[24 + i]*D[P[z], {z, i}]*D[Q[z], {z, 1 - i}], {i, 0, 
                    1}] + v[26]*D[a[z], {z, 4}]* P[z]*Q[z] + 
              v[27]*a[z]*D[a[z], {z, 2}]*P[z]*Q[z] + 
              v[28]*D[a[z], z]^2*P[z]*Q[z] + v[29]*a[z]^3* P[z]*Q[z]) // 
          Expand;

We introduce

differ  =   Expand[ 
        VAA[zet]  -  f'[zet]^(p + q + 6)*V[f[zet]] ]

Date[]

s[i_]  :=  Sum[ c[j]*m[j],  {j, 1, i}]

Date[]

m[1] = a[f[zet]] P[f[zet]] Q[f[zet]] a'[f[zet]] ((f'[zet])^(4 + p + q)) f''[zet]

Date[]

c[1]  =  Coefficient[differ,  m[1]]

Date[]

dif[1]  =  differ  -  s[1]  // Expand;

Date[]

m[2] = a[f[zet]]^2 Q[f[zet]] ((f'[zet])^(4 + p + q)) P'[f[zet]] f''[zet]

Date[]

c[2]  =  Coefficient[ differ,  m[2] ]

Date[]

dif[2]  =  differ  -  s[2]  // Expand;

Date[]

m[3] = a[f[zet]]^2 P[f[zet]] ((f'[zet])^(4 + p + q)) Q'[f[zet]] f''[zet]

Date[]

c[3]  =  Coefficient[ differ,  m[3] ]

Date[]

dif[3]  =  differ  -  s[3]  // Expand;

Date[]

m[4] =  a'[f[zet]] ((f'[zet])^(4 + p + q)) P'[f[zet]] Q'[f[zet]] f''[zet]

Date[]

c[4]  =  Coefficient[ differ, m[4] ]

Date[]

dif[4]  =  differ  -  s[4]  // Expand;

Date[]

m[5] = Q[f[zet]] ((f'[zet])^(4 + p + q)) P'[f[zet]] a''[f[zet]] f''[zet]

Date[]

c[5]  =  Coefficient[ differ,  m[5] ]

Date[]

dif[5]  =  differ  -  s[5]  // Expand;

Date[]

m[6] = P[f[zet]] ((f'[zet])^(4 + p + q)) Q'[f[zet]] a''[f[zet]] f''[zet]

Date[]

c[6]  =  Coefficient[ differ,  m[6] ]

Date[]

dif[6]  =  differ  -  s[6]  // Expand;

Date[]

m[7] = a[f[zet]]^2 P[f[zet]] Q[f[zet]] ((f'[zet])^(2 + p + q)) ((f''[zet])^2)

Date[]

c[7]  =  Coefficient[ differ,  m[7] ]

Date[]

dif[7]  =  differ  -  s[7]  // Expand;

Date[]

m[8] =  Q[f[zet]] a'[f[zet]] ((f'[zet])^(2 + p + q)) P'[f[zet]] ((f''[zet])^2)

Date[]

c[8]  =  Coefficient[ differ,  m[8] ]

Date[]

dif[8]  =  differ  -  s[8] // Expand;

Date[]

m[9] =  P[f[zet]] a'[f[zet]] ((f'[zet])^(2 + p + q)) Q'[f[zet]] ((f''[zet])^2)

Date[]

c[9]  =  Coefficient[ differ,  m[9] ]

Date[]

dif[9]  =  differ  -  s[9]  // Expand;

Date[]

m[10] =  a[f[zet]] ((f'[zet])^(2 + p + q)) P'[f[zet]] Q'[f[zet]] ((f''[zet])^2)

Date[]

c[10]  =  Coefficient[ differ,  m[10] ]

This is as far as Version 5.1 is able to carry the computation.
The  Coefficient  function is still too inefficient (although an error 
for it that was introduced in Version 5.0 has been corrected).

Date[]

dif[10]  =  differ  -  s[10]  // Expand;

m[11] = P[f[zet]] Q[f[zet]] ((f'[zet])^(2 + p + q)) a''[f[zet]] ((f''[zet])^2)

Date[]

c[11]  =  Coefficient[ differ,  m[11] ]

Date[]

dif[11]  =  differ  -  s[11]  // Expand;

Date[]

m[12] = P[f[zet]] Q[f[zet]] a'[f[zet]] ((f'[zet])^(p + q)) ((f''[zet])^3)

Date[]

c[12]  =  Coefficient[ differ,  m[12] ]

Date[]

dif[12]  =  differ  - s[12]  // Expand;

Date[]

m[13] = a[f[zet]] Q[f[zet]] ((f'[zet])^(p + q)) P'[f[zet]] ((f''[zet])^3)

Date[]

c[13]  =  Coefficient[ differ,  m[13] ]

Date[]

dif[13]  =  differ  - s[13]  // Expand;

Date[]

m[14] = a[f[zet]] P[f[zet]] ((f'[zet])^(p + q)) Q'[f[zet]] ((f''[zet])^3)

Date[]

c[14]  =  Coefficient[ differ,  m[14] ]

Date[]

dif[14]  =  differ  -  s[14]  // Expand;

Date[]

m[15] = a[f[zet]] P[f[zet]] Q[f[zet]] ((f'[zet])^(-2 + p + q)) ((f''[zet])^4)

Date[]

c[15]  =  Coefficient[ differ,  m[15] ]

Date[]

dif[15]  =  differ  -  s[15]  // Expand;

Date[]

m[16] =  ((f'[zet])^(-2 + p + q)) P'[f[zet]] Q'[f[zet]] ((f''[zet])^4)

Date[]

c[16]  =  Coefficient[ differ,  m[16] ]

Date[]

dif[16]  =  differ  -  s[16]  // Expand;

Date[]

m[17] =  Q[f[zet]] ((f'[zet])^(-4 + p + q)) P'[f[zet]] ((f''[zet])^5)

Date[]

c[17]  =  Coefficient[ differ,  m[17] ]

Date[]

dif[17]  =  differ  -  s[17]  // Expand;

Date[]

m[18] = P[f[zet]] ((f'[zet])^(-4 + p + q)) Q'[f[zet]] ((f''[zet])^5)

Date[]

c[18]  =  Coefficient[ differ,  m[18] ]

Date[]

dif[18]  =  differ  -  s[18]  // Expand;

Date[]

m[19] = P[f[zet]] Q[f[zet]] ((f'[zet])^(-6 + p + q)) ((f''[zet])^6)

Date[]

c[19]  =  Coefficient[ differ,  m[19] ]

Date[]

dif[19]  =  differ  -  s[19]  // Expand;

Date[]

m[20] = Q[f[zet]] a'[f[zet]] ((f'[zet])^(4 + p + q)) f''[zet] P''[f[zet]]

Date[]

c[20]  =  Coefficient[ differ,  m[20] ]

dif[20]  =  differ  -  s[20]  // Expand;

Date[]

m[21] = a[f[zet]] ((f'[zet])^(4 + p + q)) Q'[f[zet]] f''[zet] P''[f[zet]]

Date[]

c[21]  =  Coefficient[ differ,  m[21] ]

Date[]

dif[21]  =  differ  -  s[21]  // Expand;

Date[]

m[22] = a[f[zet]] Q[f[zet]] ((f'[zet])^(2 + p + q)) ((f''[zet])^2) P''[f[zet]]

Date[]

c[22]  =  Coefficient[ differ,  m[22] ]

Date[]

dif[22]  =  differ  -  s[22]  // Expand;

Date[]

m[23] = ((f'[zet])^(p + q)) Q'[f[zet]] ((f''[zet])^3) P''[f[zet]]

Date[]

c[23]  =  Coefficient[ differ,  m[23] ]

Date[]

dif[23]  =  differ  -  s[23]  // Expand;

Date[]

m[24] = Q[f[zet]] ((f'[zet])^(-2 + p + q)) ((f''[zet])^4) P''[f[zet]]

Date[]

c[24]  =  Coefficient[ differ,  m[24] ]

Date[]

dif[24]  =  differ  -  s[24]  // Expand;

Date[]

m[25] = P[f[zet]] a'[f[zet]] ((f'[zet])^(4 + p + q)) f''[zet] Q''[f[zet]]

Date[]

c[25]  =  Coefficient[ differ,  m[25] ]

Date[]

dif[25]  =  differ  -  s[25]  // Expand;

Date[]

m[26] = a[f[zet]] ((f'[zet])^(4 + p + q)) P'[f[zet]] f''[zet] Q''[f[zet]]

Date[]

c[26]  =  Coefficient[ differ,  m[26] ]

Date[]

dif[26]  =  differ  -  s[26]  // Expand;

Date[]

m[27] = a[f[zet]] P[f[zet]] ((f'[zet])^(2 + p + q)) ((f''[zet])^2) Q''[f[zet]]

Date[]

c[27]  =  Coefficient[ differ,  m[27] ]

Date[]

dif[27]  =  differ  -  s[27]  // Expand;

Date[]

m[28] = ((f'[zet])^(p + q)) P'[f[zet]] ((f''[zet])^3) Q''[f[zet]]

Date[]

c[28]  =  Coefficient[ differ,  m[28] ]

Date[]

dif[28]  =  differ  -  s[28]  // Expand;

Date[]

m[29] = P[f[zet]] ((f'[zet])^(-2 + p + q)) ((f''[zet])^4) Q''[f[zet]]

Date[]

c[29]  =  Coefficient[ differ,  m[29] ]

Date[]

dif[29]  =  differ  -  s[29]  // Expand;

Date[]

m[30] = ((f'[zet])^(2 + p + q)) ((f''[zet])^2) P''[f[zet]] Q''[f[zet]]

Date[]

c[30]  =  Coefficient[ differ,  m[30] ]

Date[]

dif[30]  =  differ  -  s[30]  // Expand;

Date[]

m[31] = P[f[zet]] Q[f[zet]] ((f'[zet])^(4 + p + q)) f''[zet] Derivative[3][a][f[zet]]

Date[]

c[31]  =  Coefficient[ differ,  m[31] ]

Date[]

dif[31]  =  differ  -  s[31]  // Expand;

Date[]

m[32] = a[f[zet]]^2 P[f[zet]] Q[f[zet]] ((f'[zet])^(3 + p + q)) Derivative[3][f][zet]

Date[]

c[32]  =  Coefficient[ differ,  m[32] ]

Date[]

dif[32]  =  differ  -  s[32]  // Expand;

Date[]

m[33] = Q[f[zet]] a'[f[zet]] ((f'[zet])^(3 + p + q)) P'[f[zet]] Derivative[3][f][zet]

Date[]

c[33]  =  Coefficient[ differ,  m[33] ]

Date[]

dif[33]  =  differ  -  s[33]  // Expand;

Date[]

m[34] = P[f[zet]] a'[f[zet]] ((f'[zet])^(3 + p + q)) Q'[f[zet]] Derivative[3][f][zet]

Date[]

c[34]  =  Coefficient[ differ,  m[34] ]

Date[]

dif[34]  =  differ  -  s[34]  // Expand;

Date[]

m[35] = a[f[zet]] ((f'[zet])^(3 + p + q)) P'[f[zet]] Q'[f[zet]] Derivative[3][f][zet]

Date[]

c[35]  =  Coefficient[ differ,  m[35] ]

Date[]

dif[35]  =  differ  -  s[35]  // Expand;

Date[]

m[36] = P[f[zet]] Q[f[zet]] ((f'[zet])^(3 + p + q)) a''[f[zet]] Derivative[3][f][zet]

Date[]

c[36]  =  Coefficient[ differ,  m[36] ]

Date[]

dif[36]  =  differ  -  s[36]  // Expand;

Date[]

m[37] = P[f[zet]] Q[f[zet]] a'[f[zet]] ((f'[zet])^(1 + p + q)) f''[zet] Derivative[3][f][zet]

Date[]

c[37]  =  Coefficient[ differ,  m[37] ]

Date[]

dif[37]  =  differ  -  s[37]  // Expand;

Date[]

m[38] = a[f[zet]] Q[f[zet]] ((f'[zet])^(1 + p + q)) P'[f[zet]] f''[zet] Derivative[3][f][zet]

Date[]

c[38]  =  Coefficient[ differ,  m[38] ]

Date[]

dif[38]  =  differ  -  s[38]  // Expand;

Date[]

m[39] = a[f[zet]] P[f[zet]] ((f'[zet])^(1 + p + q)) Q'[f[zet]] f''[zet] Derivative[3][f][zet]

Date[]

c[39]  =  Coefficient[ differ,  m[39] ]

Date[]

dif[39]  =  differ  -  s[39]  // Expand;

Date[]

m[40] = a[f[zet]] P[f[zet]] Q[f[zet]] ((f'[zet])^(-1 + p + q)) ((f''[zet])^2) Derivative[3][f][zet]

Date[]

c[40]  =  Coefficient[ differ,  m[40] ]

Date[]

dif[40]  =  differ  -  s[40]  // Expand;

Date[]

m[41] = ((f'[zet])^(-1 + p + q)) P'[f[zet]] Q'[f[zet]] ((f''[zet])^2) Derivative[3][f][zet]

Date[]

c[41]  =  Coefficient[ differ,  m[41] ]

Date[]

dif[41]  =  differ  -  s[41]  // Expand;

Date[]

m[42] =  Q[f[zet]] ((f'[zet])^(-3 + p + q)) P'[f[zet]] ((f''[zet])^3) Derivative[3][f][zet]

Date[]

c[42]  =  Coefficient[ differ,  m[42] ]

Date[]

dif[42]  =  differ  -  s[42]  // Expand;

Date[]

m[43] = P[f[zet]] ((f'[zet])^(-3 + p + q)) Q'[f[zet]] ((f''[zet])^3) Derivative[3][f][zet]

Date[]

c[43]  =  Coefficient[ differ,  m[43] ]

Date[]

dif[43]  =  differ  -  s[43]  // Expand;

Date[]

m[44] = P[f[zet]] Q[f[zet]] ((f'[zet])^(-5 + p + 
            q)) ((f''[zet])^4) Derivative[3][f][zet]

Date[]

c[44]  =  Coefficient[ differ,  m[44] ]

Date[]

dif[44]  =  differ  -  s[44]  // Expand;

Date[]

m[45] = a[f[zet]] Q[f[zet]] ((f'[zet])^(3 + p + q)) P''[f[zet]] Derivative[3][f][zet]

Date[]

c[45]  =  Coefficient[ differ,  m[45] ]

Date[]

dif[45]  =  differ  -  s[45]  // Expand;

Date[]

m[46] =  ((f'[zet])^(1 + p + q)) Q'[f[zet]] f''[zet] P''[f[zet]] Derivative[3][f][zet]

Date[]

c[46]  =  Coefficient[ differ,  m[46] ]

Date[]

dif[46]  =  differ  -  s[46]  // Expand;

Date[]

m[47] =  Q[f[zet]] ((f'[zet])^(-1 + p + q)) ((f''[zet])^2) P''[f[zet]] Derivative[3][f][zet]

Date[]

c[47]  =  Coefficient[ differ,  m[47] ]

Date[]

dif[47]  =  differ  -  s[47]  // Expand;

Date[]

m[48] = a[f[zet]] P[f[zet]] ((f'[zet])^(3 + p + q)) Q''[f[zet]] Derivative[3][f][zet]

Date[]

c[48]  =  Coefficient[ differ,  m[48] ]

Date[]

dif[48]  =  differ  -  s[48]  // Expand;

Date[]

m[49] = ((f'[zet])^(1 + p + q)) P'[f[zet]] f''[zet] Q''[f[zet]] Derivative[3][f][zet]

Date[]

c[49]  =  Coefficient[ differ,  m[49] ]

Date[]

dif[49]  =  differ  -  s[49]  // Expand;

Date[]

m[50] = P[f[zet]] ((f'[zet])^(-1 + p + q)) ((f''[zet])^2) Q''[f[zet]] Derivative[3][f][zet]

Date[]

c[50]  =  Coefficient[ differ,  m[50] ]

Date[]

dif[50]  =  differ  -  s[50]  // Expand;

Date[]

m[51] = ((f'[zet])^(3 + p + q)) P''[f[zet]] Q''[f[zet]] Derivative[3][f][zet]

Date[]

c[51]  =  Coefficient[ differ,  m[51] ]

Date[]

dif[51]  =  differ  -  s[51]  // Expand;

Date[]

m[52] =  a[f[zet]] P[f[zet]] Q[f[zet]] ((f'[zet])^(p + q)) ((Derivative[3][f][zet])^2)

Date[]

c[52]  =  Coefficient[ differ,  m[52] ]

Date[]

dif[52]  =  differ  -  s[52]  // Expand;

Date[]

m[53] = ((f'[zet])^(p + q)) P'[f[zet]] Q'[f[zet]] ((Derivative[3][f][zet])^2)

Date[]

c[53]  =  Coefficient[ differ,  m[53] ]

Date[]

dif[53]  =  differ  -  s[53]  // Expand;

Date[]

m[54] = Q[f[zet]] ((f'[zet])^(-2 + p + q)) P'[f[zet]] f''[zet] ((Derivative[3][f][zet])^2)

Date[]

c[54]  =  Coefficient[ differ,  m[54] ]

Date[]

dif[54]  =  differ  -  s[54]  // Expand;

Date[]

m[55] = P[f[zet]] ((f'[zet])^(-2 + p + q)) Q'[f[zet]] f''[zet] ((Derivative[3][f][zet])^2)

Date[]

c[55]  =  Coefficient[ differ,  m[55] ]

Date[]

dif[55]  =  differ  -  s[55]  // Expand;

Date[]

m[56] = P[f[zet]] Q[f[zet]] ((f'[zet])^(-4 + p + q)) ((f''[zet])^2) ((Derivative[3][f][zet])^2)

Date[]

c[56]  =  Coefficient[ differ,  m[56] ]

Date[]

dif[56]  =  differ  -  s[56]  // Expand;

Date[]

m[57] = Q[f[zet]] ((f'[zet])^(p + q)) P''[f[zet]] ((Derivative[3][f][zet])^2)

Date[]

c[57]  =  Coefficient[ differ,  m[57] ]

Date[]

dif[57]  =  differ  -  s[57]  // Expand;

Date[]

m[58] =  P[f[zet]] ((f'[zet])^(p + q)) Q''[f[zet]] ((Derivative[3][f][zet])^2)

Date[]

c[58]  =  Coefficient[ differ,  m[58] ]

Date[]

dif[58]  =  differ  -  s[58]  // Expand;

Date[]

m[59] = P[f[zet]] Q[f[zet]] ((f'[zet])^(-3 + p + q)) ((Derivative[3][f][zet])^3)

Date[]

c[59]  =  Coefficient[ differ,  m[59] ]

Date[]

dif[59]  =  differ  -  s[59]  // Expand;

Date[]

m[60] = a[f[zet]] Q[f[zet]] ((f'[zet])^(4 + p + q)) f''[zet] Derivative[3][P][f[zet]]

Date[]

c[60]  =  Coefficient[ differ,  m[60] ]

Date[]

dif[60]  =  differ  -  s[60]  // Expand;

Date[]

m[61] =  ((f'[zet])^(2 + p + q)) Q'[f[zet]] ((f''[zet])^2) Derivative[3][P][f[zet]]

Date[]

c[61]  =  Coefficient[ differ,  m[61] ]

Date[]

dif[61]  =  differ  -  s[61]  // Expand;

Date[]

m[62] =  Q[f[zet]] ((f'[zet])^(p + q)) ((f''[zet])^3) Derivative[3][P][f[zet]]

Date[]

c[62]  =  Coefficient[ differ,  m[62] ]

Date[]

dif[62]  =  differ  -  s[62]  // Expand;

Date[]

m[63] = ((f'[zet])^(4 + p + q)) f''[zet] Q''[f[zet]] Derivative[3][P][f[zet]]

Date[]

c[63]  =  Coefficient[ differ,  m[63] ]

Date[]

dif[63]  =  differ  -  s[63]  // Expand;

Date[]

m[64] = ((f'[zet])^(3 + p + q)) Q'[f[zet]] Derivative[3][f][zet] Derivative[3][P][f[zet]]

Date[]

c[64]  =  Coefficient[ differ,  m[64] ]

Date[]

dif[64]  =  differ  -  s[64]  // Expand;

Date[]

m[65] = Q[f[zet]] ((f'[zet])^(1 + p + q)) f''[zet] Derivative[3][f][zet] Derivative[3][P][f[zet]]

Date[]

c[65]  =  Coefficient[ differ,  m[65] ]

Date[]

dif[65]  =  differ  -  s[65]  // Expand;

Date[]

m[66] = a[f[zet]] P[f[zet]] ((f'[zet])^(4 + p + q)) f''[zet] Derivative[3][Q][f[zet]]

Date[]

c[66]  =  Coefficient[ differ,  m[66] ]

Date[]

dif[66]  =  differ  -  s[66]  // Expand;

Date[]

m[67] = ((f'[zet])^(2 + p + q)) P'[f[zet]] ((f''[zet])^2) Derivative[3][Q][f[zet]]

Date[]

c[67]  =  Coefficient[ differ,  m[67] ]

Date[]

dif[67]  =  differ  -  s[67]  // Expand;

Date[]

m[68] = P[f[zet]] ((f'[zet])^(p + q)) ((f''[zet])^3) Derivative[3][Q][f[zet]]

Date[]

c[68]  =  Coefficient[ differ,  m[68] ]

Date[]

dif[68]  =  differ  -  s[68]  // Expand;

Date[]

m[69] = ((f'[zet])^(4 + p + q)) f''[zet] P''[f[zet]] Derivative[3][Q][f[zet]]

Date[]

c[69]  =  Coefficient[ differ,  m[69] ]

Date[]

dif[69]  =  differ  -  s[69]  // Expand;

Date[]

m[70] = ((f'[zet])^(3 + p + q)) P'[f[zet]] Derivative[3][f][zet] Derivative[3][Q][f[zet]]

Date[]

c[70]  =  Coefficient[ differ,  m[70] ]

Date[]

dif[70]  =  differ  -  s[70]  // Expand;

Date[]

m[71] = P[f[zet]] ((f'[zet])^(1 + p + q)) f''[zet] Derivative[3][f][zet] Derivative[3][Q][f[zet]]

Date[]

c[71]  =  Coefficient[ differ,  m[71] ]

Date[]

dif[71]  =  differ  -  s[71]  // Expand;

Date[]

m[72] =  P[f[zet]] Q[f[zet]] a'[f[zet]] ((f'[zet])^(2 + p + q)) Derivative[4][f][zet]

Date[]

c[72]  =  Coefficient[ differ,  m[72] ]

Date[]

dif[72]  =  differ  -  s[72]  // Expand;

Date[]

m[73] = a[f[zet]] Q[f[zet]] ((f'[zet])^(2 + p + q)) P'[f[zet]] Derivative[4][f][zet]

Date[]

c[73]  =  Coefficient[ differ,  m[73] ]

Date[]

dif[73]  =  differ  -  s[73]  // Expand;

Date[]

m[74] = a[f[zet]] P[f[zet]] ((f'[zet])^(2 + p + q)) Q'[f[zet]] Derivative[4][f][zet]

Date[]

c[74]  =  Coefficient[ differ,  m[74] ]

Date[]

dif[74]  =  differ  -  s[74]  // Expand;

Date[]

m[75] = a[f[zet]] P[f[zet]] Q[f[zet]] ((f'[zet])^(p + q)) f''[zet] Derivative[4][f][zet]

Date[]

c[75]  =  Coefficient[ differ,  m[75] ]

Date[]

dif[75]  =  differ  -  s[75]  // Expand;

Date[]

m[76] = ((f'[zet])^(p + q)) P'[f[zet]] Q'[f[zet]] f''[zet] Derivative[4][f][zet]

Date[]

c[76]  =  Coefficient[ differ,  m[76] ]

Date[]

dif[76]  =  differ  -  s[76]  // Expand;

Date[]

m[77] = Q[f[zet]] ((f'[zet])^(-2 + p + q)) P'[f[zet]] ((f''[zet])^2) Derivative[4][f][zet]

Date[]

c[77]  =  Coefficient[ differ,  m[77] ]

Date[]

dif[77]  =  differ  -  s[77]  // Expand;

Date[]

m[78] = P[f[zet]] ((f'[zet])^(-2 + p + q)) Q'[f[zet]] ((f''[zet])^2) Derivative[4][f][zet]

Date[]

c[78]  =  Coefficient[ differ,  m[78] ]

Date[]

dif[78]  =  differ  -  s[78]  // Expand;

Date[]

m[79] = P[f[zet]] Q[f[zet]] ((f'[zet])^(-4 + p + q)) ((f''[zet])^3) Derivative[4][f][zet]

Date[]

c[79]  =  Coefficient[ differ,  m[79] ]

Date[]

dif[79]  =  differ  -  s[79]  // Expand;

Date[]

m[80] = ((f'[zet])^(2 + p + q)) Q'[f[zet]] P''[f[zet]] Derivative[4][f][zet]

Date[]

c[80]  =  Coefficient[ differ,  m[80] ]

Date[]

dif[80]  =  differ  -  s[80]  // Expand;

Date[]

m[81] = Q[f[zet]] ((f'[zet])^(p + q)) f''[zet] P''[f[zet]] Derivative[4][f][zet]

Date[]

c[81]  =  Coefficient[ differ,  m[81] ]

Date[]

dif[81]  =  differ  -  s[81]  // Expand;

Date[]

m[82] = ((f'[zet])^(2 + p + q)) P'[f[zet]] Q''[f[zet]] Derivative[4][f][zet]

Date[]

c[82]  =  Coefficient[ differ,  m[82] ]

Date[]

dif[82]  =  differ  -  s[82]  // Expand;

Date[]

m[83] = P[f[zet]] ((f'[zet])^(p + q)) f''[zet] Q''[f[zet]] Derivative[4][f][zet]

Date[]

c[83]  =  Coefficient[ differ,  m[83] ]

Date[]

dif[83]  =  differ  -  s[83]  // Expand;

Date[]

m[84] = Q[f[zet]] ((f'[zet])^(-1 + p + q)) P'[f[zet]] Derivative[3][f][zet] Derivative[4][f][zet]

Date[]

c[84]  =  Coefficient[ differ,  m[84] ]

Date[]

dif[84]  =  differ  -  s[84]  // Expand;

Date[]

m[85] = P[f[zet]] ((f'[zet])^(-1 + p + q)) Q'[f[zet]] Derivative[3][f][zet] Derivative[4][f][zet]

Date[]

c[85]  =  Coefficient[ differ,  m[85] ]

Date[]

dif[85]  =  differ  -  s[85]  // Expand;

Date[]

m[86] = P[f[zet]] Q[f[zet]] ((f'[zet])^(-3 + p + q)) f''[zet] Derivative[3][f][zet] Derivative[4][f][zet]

Date[]

c[86]  =  Coefficient[ differ,  m[86] ]

Date[]

dif[86]  =  differ  -  s[86]  // Expand;

Date[]

m[87] = Q[f[zet]] ((f'[zet])^(2 + p + q)) Derivative[3][P][f[zet]] Derivative[4][f][zet]

Date[]

c[87]  =  Coefficient[ differ,  m[87] ]

Date[]

dif[87]  =  differ  -  s[87]  // Expand;

Date[]

m[88] =  P[f[zet]] ((f'[zet])^(2 + p + q)) Derivative[3][Q][f[zet]] Derivative[4][f][zet]

Date[]

c[88]  =  Coefficient[ differ,  m[88] ]

Date[]

dif[88]  =  differ  -  s[88]  // Expand;

Date[]

m[89] = P[f[zet]] Q[f[zet]] ((f'[zet])^(-2 + p + q)) ((Derivative[4][f][zet])^2)

Date[]

c[89]  =  Coefficient[ differ,  m[89] ]

Date[]

dif[89]  =  differ  -  s[89]  // Expand;

Date[]

m[90] = ((f'[zet])^(4 + p + q)) Q'[f[zet]] f''[zet] Derivative[4][P][f[zet]]

Date[]

c[90]  =  Coefficient[ differ,  m[90] ]

Date[]

dif[90]  =  differ  -  s[90]  // Expand;

Date[]

m[91] =  Q[f[zet]] ((f'[zet])^(2 + p + q)) ((f''[zet])^2) Derivative[4][P][f[zet]]

Date[]

c[91]  =  Coefficient[ differ,  m[91] ]

Date[]

dif[91]  =  differ  -  s[91]  // Expand;

Date[]

m[92] =  Q[f[zet]] ((f'[zet])^(3 + p + q)) Derivative[3][f][zet] Derivative[4][P][f[zet]]

Date[]

c[92]  =  Coefficient[ differ,  m[92] ]

Date[]

Date[]

dif[92]  =  differ  -  s[92]  // Expand;

Date[]

m[93] = ((f'[zet])^(4 + p + q)) P'[f[zet]] f''[zet] Derivative[4][Q][f[zet]]

Date[]

c[93]  =  Coefficient[ differ,  m[93] ]

Date[]

dif[93]  =  differ  -  s[93]  // Expand;

m[94] = P[f[zet]] ((f'[zet])^(2 + p + q)) ((f''[zet])^2) Derivative[4][Q][f[zet]]

Date[]

c[94]  =  Coefficient[ differ,  m[94] ]

Date[]

dif[94]  =  differ  -  s[94]  // Expand;

Date[]

m[95] = P[f[zet]] ((f'[zet])^(3 + p + q)) Derivative[3][f][zet] Derivative[4][Q][f[zet]]

Date[]

c[95]  =  Coefficient[ differ,  m[95] ]

Date[]

dif[95]  =  differ  -  s[95]  // Expand;

Date[]

m[96] = a[f[zet]] P[f[zet]] Q[f[zet]] ((f'[zet])^(1 + p + q)) Derivative[5][f][zet]

Date[]

c[96]  =  Coefficient[ differ,  m[96] ]

Date[]

dif[96]  =  differ  -  s[96]  // Expand;

Date[]

m[97] = ((f'[zet])^(1 + p + q)) P'[f[zet]] Q'[f[zet]] Derivative[5][f][zet]

Date[]

c[97]  =  Coefficient[ differ,  m[97] ]

Date[]

dif[97]  =  differ  -  s[97]  // Expand;

Date[]

m[98] =  Q[f[zet]] ((f'[zet])^(-1 + p + q)) P'[f[zet]] f''[zet] Derivative[5][f][zet]

Date[]

c[98]  =  Coefficient[ differ,  m[98] ]

Date[]

dif[98]  =  differ  -  s[98]  // Expand;

Date[]

m[99] = P[f[zet]] ((f'[zet])^(-1 + p + q)) Q'[f[zet]] f''[zet] Derivative[5][f][zet]

Date[]

c[99]  =  Coefficient[ differ,  m[99] ]

Date[]

dif[99]  =  differ  -  s[99]  // Expand;

Date[]

m[100] = P[f[zet]] Q[f[zet]] ((f'[zet])^(-3 + p + q)) ((f''[zet])^2) Derivative[5][f][zet]

Date[]

c[100]  =  Coefficient[ differ,  m[100] ]

Date[]

dif[100]  =  differ  -  s[100]  // Expand;

Date[]

m[101] = Q[f[zet]] ((f'[zet])^(1 + p + q)) P''[f[zet]] Derivative[5][f][zet]

Date[]

c[101]  =  Coefficient[ differ,  m[101] ]

Date[]

dif[101]  =  differ  -  s[101]  // Expand;

Date[]

m[102] = P[f[zet]] ((f'[zet])^(1 + p + q)) Q''[f[zet]] Derivative[5][f][zet]

Date[]

c[102]  =  Coefficient[ differ,  m[102] ]

Date[]

dif[102]  =  differ  -  s[102]  // Expand;

Date[]

m[103] = P[f[zet]] Q[f[zet]] ((f'[zet])^(-2 + p + q)) Derivative[3][f][zet] Derivative[5][f][zet]

Date[]

c[103]  =  Coefficient[ differ,  m[103] ]

Date[]

dif[103]  =  differ  -  s[103]  // Expand;

Date[]

m[104] = Q[f[zet]] ((f'[zet])^(4 + p + q)) f''[zet] Derivative[5][P][f[zet]]

Date[]

c[104]  =  Coefficient[ differ,  m[104] ]

Date[]

dif[104]  =  differ  -  s[104]  // Expand;

Date[]

m[105] =  P[f[zet]] ((f'[zet])^(4 + p + q)) f''[zet] Derivative[5][Q][f[zet]]

Date[]

c[105]  =  Coefficient[ differ,  m[105] ]

Date[]

dif[105]  =  differ  -  s[105]  // Expand;

Date[]

m[106] = Q[f[zet]] ((f'[zet])^(p + q)) P'[f[zet]] Derivative[6][f][zet]

Date[]

c[106]  =  Coefficient[ differ,  m[106] ]

Date[]

dif[106]  =  differ  -  s[106]  // Expand;

Date[]

m[107] = P[f[zet]] ((f'[zet])^(p + q)) Q'[f[zet]] Derivative[6][f][zet]

Date[]

c[107]  =  Coefficient[ differ,  m[107] ]

Date[]

dif[107]  =  differ  -  s[107]  // Expand;

Date[]

m[108] = P[f[zet]] Q[f[zet]] ((f'[zet])^(-2 + p + q)) f''[zet] Derivative[6][f][zet]

Date[]

c[108]  =  Coefficient[ differ,  m[108] ]

Date[]

dif[108]  =  differ  -  s[108]  // Expand;

Date[]

m[109] = P[f[zet]] Q[f[zet]] ((f'[zet])^(-1 + p + q)) Derivative[7][f][zet]

Date[]

c[109]  =  Coefficient[ differ,  m[109] ]

Date[]

dif[109]  =  differ  -  s[109]  // Expand

Date[]

Thus, we have

Expand[  differ  -  Sum[ c[i]*m[i],  {i, 1, 109}] ]

Date[]

eqs  =  Table[ c[i] ==  0,  {i, 1, 109} ];

unknowns  =  Table[ v[i],  {i, 1, 29}]

In order to have  differ = 0, it is necessary and sufficient that   y1, y2, y3, y4, y5, y6, y7, y8, y9, y10, y11   be solutions of the system

sol  = Solve[ eqs,  unknowns]

Date[]

As a check, we set

v1  =  v1   /.  sol[[1]]

v2  =  v2  /.  sol[[1]]  // Factor

v3  =  v3  /.  sol[[1]]

v4  =  v4  /.  sol[[1]]  // Factor

v5  =  v5  /.  sol[[1]]  // Factor

v6  =  v6  /.  sol[[1]]  // Factor

v7  =  v7  /.  sol[[1]]  // Factor

v8  =  v8  /.  sol[[1]]  // Factor

v9  =  v9  /.  sol[[1]]  // Factor

v10  =  v10  /.  sol[[1]]  // Factor

v11  =  v11  /.  sol[[1]]  // Factor

v12  =  v12  /.  sol[[1]]   // Factor

v13  =  v13  /.  sol[[1]]   // Factor

v14  =  v14  /.  sol[[1]]   // Factor

v15  =  v15  /.  sol[[1]]   // Factor

v16  =  v16  /.  sol[[1]]   // Factor

v17  =  v17  /.  sol[[1]]   // Factor

v18  =  v18  /.  sol[[1]]   // Factor

v19  =  v19  /.  sol[[1]]   // Factor

v20  =  v20  /.  sol[[1]]   // Factor

v21 =  v21  /.  sol[[1]]   // Factor

v22 =  v22  /.  sol[[1]]   // Factor

v23 =  v23  /.  sol[[1]]   // Factor

v24 =  v24  /.  sol[[1]]   // Factor

v25 =  v25  /.  sol[[1]]   // Factor

v26 =  v26  /.  sol[[1]]   // Factor

v27 =  v27  /.  sol[[1]]   // Factor

v28 =  v28  /.  sol[[1]]   // Factor

v29 =  v29  /.  sol[[1]]   // Factor

Date[]

differ // FullSimplify

Date[]

V[z]

This yields (14.48) of page 153.

MaxMemoryUsed[]

System Information report written at:02/22/2005 04:09:43 PM
[System Summary]

Item Value
OS Name Microsoft Windows 2000 Professional
Version 5.0.2195 Service Pack 4 Build 2195
OS Manufacturer Microsoft Corporation
System Name CHALKLEY
System Manufacturer Dell Computer Corporation
System Model XPS-Z
System Type X86-based PC
Processor x86 Family 6 Model 8 Stepping 3 GenuineIntel~930 Mhz
BIOS Version 08/29/00
Windows Directory C:\WINNT
System Directory C:\WINNT\system32
Boot Device\Device\Harddisk0\Partition1
Locale United States
User Name CHALKLEY\chalklr
Time Zone Eastern Standard Time
Total Physical Memory 523,568 KB
Available Physical Memory 260,668 KB
Total Virtual Memory 1,801,932 KB
Available Virtual Memory 1,256,860 KB
Page File Space 1,278,364 KB
Page File C:\pagefile.sys