Quantcast
Viewing all articles
Browse latest Browse all 1085

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

Image may be NSFW.
Clik here to view.
restart

Image may be NSFW.
Clik here to view.
with(PDEtools)

Image may be NSFW.
Clik here to view.
undeclare(prime)

Image may be NSFW.
Clik here to view.
`There is no more prime differentiation variable; all derivatives will be displayed as indexed functions`

(1)

Image may be NSFW.
Clik here to view.
with(DEtools)

Image may be NSFW.
Clik here to view.
with(DifferentialAlgebra)

Image may be NSFW.
Clik here to view.
"with(Student[ODEs][Solve]): "

Image may be NSFW.
Clik here to view.
with(IntegrationTools)

Image may be NSFW.
Clik here to view.
with(inttrans)

Image may be NSFW.
Clik here to view.
with(PDEtools)

Image may be NSFW.
Clik here to view.
with(Physics)

Image may be NSFW.
Clik here to view.
with(PolynomialTools)

Image may be NSFW.
Clik here to view.
with(RootFinding)

Image may be NSFW.
Clik here to view.
with(SolveTools)

Image may be NSFW.
Clik here to view.
with(LinearAlgebra)

Image may be NSFW.
Clik here to view.
with(sumtools)

Image may be NSFW.
Clik here to view.
``

Image may be NSFW.
Clik here to view.
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

Image may be NSFW.
Clik here to view.
NULL

Image may be NSFW.
Clik here to view.
L := convert((cosh(xi)+sinh(xi))/(cosh(xi)-sinh(xi)), trig)

Image may be NSFW.
Clik here to view.
"Q(xi):=L:"

Image may be NSFW.
Clik here to view.
S := sum(A[i]*Q(xi)^i, i = 0 .. 1)+sum(B[i]*Q(xi)^(-i), i = 1 .. 1)

Image may be NSFW.
Clik here to view.
``

(2)

Image may be NSFW.
Clik here to view.
S

Image may be NSFW.
Clik here to view.
K := F(xi) = S

Image may be NSFW.
Clik here to view.
F1 := eval(ode, K)

Image may be NSFW.
Clik here to view.
simplify(%)

Image may be NSFW.
Clik here to view.
P := numer(lhs())*denom(rhs()) = numer(rhs())*denom(lhs())

Warning,  computation interrupted

 

Image may be NSFW.
Clik here to view.
NULL

Image may be NSFW.
Clik here to view.
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

 

Image may be NSFW.
Clik here to view.
{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


Viewing all articles
Browse latest Browse all 1085

Trending Articles