Question: How choose the case after finding parameter

when i finding parameter i want just choose a case for example a_1=a_1  and any other case a_2=0,and remove other case how i can do in maple

restart

with(PDEtools)

undeclare(prime)

`There is no more prime differentiation variable; all derivatives will be displayed as indexed functions`

(1)

with(DEtools)

with(DifferentialAlgebra)

"with(Student[ODEs][Solve]): "

with(IntegrationTools)

with(inttrans)

with(PDEtools)

with(Physics)

with(PolynomialTools)

with(RootFinding)

with(SolveTools)

with(LinearAlgebra)

with(sumtools)

``

ode := F(xi)^5*a[4]+F(xi)^4*a[3]+F(xi)^3*a[2]+(-k^2*a[1]+(diff(diff(F(xi), xi), xi))*a[5]-w)*F(xi)^2+(1/2)*F(xi)*(diff(diff(F(xi), xi), xi))*a[1]-(1/4)*(diff(F(xi), xi))^2*a[1] = 0

NULL

L := convert((cosh(xi)+sinh(xi))/(cosh(xi)-sinh(xi)), trig)

"Q(xi):=L:"

S := sum(A[i]*Q(xi)^i, i = 0 .. 1)+sum(B[i]*Q(xi)^(-i), i = 1 .. 1)

``

(2)

S

K := F(xi) = S

F1 := eval(ode, K)

simplify(%)

P := numer(lhs())*denom(rhs()) = numer(rhs())*denom(lhs())

Warning,  computation interrupted

 

NULL

solve(identity(P, xi), {k, w, A[0], A[1], B[1], a[1], a[2], a[3], a[4], a[5]})

Warning, solutions may have been lost

 

{k = k, w = w, A[0] = 0, A[1] = A[1], B[1] = 0, a[1] = a[1], a[2] = a[2], a[3] = a[3], a[4] = a[4], a[5] = a[5]}, {k = k, w = -4*A[0]*a[5], A[0] = A[0], A[1] = A[1], B[1] = B[1], a[1] = 0, a[2] = -4*a[5], a[3] = 0, a[4] = 0, a[5] = a[5]}, {k = k, w = (1/2)*A[0]*(3*k^2*A[0]^2*a[4]+2*k^2*A[0]*a[3]+k^2*a[2]+4*k^2*a[5]+2*A[0]^2*a[4]+2*A[0]*a[3]+2*a[2]), A[0] = A[0], A[1] = 0, B[1] = 0, a[1] = -(1/2)*A[0]*(3*A[0]^2*a[4]+2*A[0]*a[3]+a[2]+4*a[5]), a[2] = a[2], a[3] = a[3], a[4] = a[4], a[5] = a[5]}, {k = k, w = w, A[0] = A[0], A[1] = 0, B[1] = 0, a[1] = a[1], a[2] = (-A[0]^3*a[4]+k^2*a[1]-A[0]^2*a[3]+w)/A[0], a[3] = a[3], a[4] = a[4], a[5] = a[5]}, {k = k, w = 4*A[1]*a[5]+4*B[1]*a[5], A[0] = -A[1]-B[1], A[1] = A[1], B[1] = B[1], a[1] = 0, a[2] = -4*a[5], a[3] = 0, a[4] = 0, a[5] = a[5]}, {k = k, w = -k^2*a[1]-4*A[0]*a[5]+a[1], A[0] = A[0], A[1] = (1/4)*A[0]^2/B[1], B[1] = B[1], a[1] = a[1], a[2] = -4*a[5], a[3] = 0, a[4] = 0, a[5] = a[5]}, {k = k, w = w, A[0] = 2*B[1], A[1] = B[1], B[1] = B[1], a[1] = a[1], a[2] = (1/2)*(k^2*a[1]+w-a[1])/B[1], a[3] = 0, a[4] = 0, a[5] = -(1/8)*(k^2*a[1]+w-a[1])/B[1]}, {k = k, w = w, A[0] = A[0], A[1] = B[1], B[1] = B[1], a[1] = 0, a[2] = w/A[0], a[3] = 0, a[4] = 0, a[5] = -(1/4)*w/A[0]}, {k = k, w = 0, A[0] = 0, A[1] = B[1], B[1] = B[1], a[1] = 0, a[2] = a[2], a[3] = 0, a[4] = 0, a[5] = -(1/4)*a[2]}

(3)

Download choose_case.mw

Please Wait...