|
LinearFunctionalSystems[HomogeneousSystem] - return the homogeneous system equivalent to the given system
|
|
Calling Sequence
|
|
HomogeneousSystem(homo, sys, vars)
HomogeneousSystem(homo, A, b, x, case)
HomogeneousSystem(homo, A, x, case)
|
|
|
Parameters
|
|
|
homo
|
-
|
name; indicates whether the given system is homogeneous
|
|
sys
|
-
|
list of equations; linear functional system
|
|
vars
|
-
|
list of function variables such as ![[y1(x), y2(x), () .. ()]](/support/helpjp/helpview.aspx?si=8752/file02036/math76.png) ; variables to solve for
|
|
A
|
-
|
Matrix with rational elements
|
|
b
|
-
|
Vector with rational elements
|
|
x
|
-
|
independent variable
|
|
case
|
-
|
name indicating the case of the system; one of 'differential', 'difference', or 'qdifference'
|
|
|
|
|
|
Description
|
|
|
•
|
The HomogeneousSystem function returns the homogeneous system that is equivalent to the given linear functional system of equations with polynomial coefficients.
|
|
|
The system parameter is entered either in list form (a list of equations sys and a list of function variables vars to solve for), or in matrix form (matrix A, vector b, and the independent variable x, where the vector b is optional).
|
|
|
The matrix form specifies the system  , where L is the operator (either differential, difference, or q-difference),  is the vector of the functions to solve for, A is a rational matrix, and b is a rational vector (right-hand side).
|
|
|
For the matrix from of the calling sequence, the case of the system must be specified as one of 'differential', 'difference', or 'qdifference'. If b is not specified, the system is assumed to be homogeneous.
|
|
|
The parameter homo is set to true if the given system is homogeneous and false otherwise.
|
|
•
|
All of the base functions of the LinearFunctionalSystems package are designed to handle homogeneous systems. In order to solve inhomogeneous systems, the following technique is used.
|
|
|
For a given inhomogeneous system, HomogeneousSystem constructs the equivalent homogeneous system by adding one more equation for the constant function and transforming all right-hand sides of the given system into the coefficients of this new function.
|
|
•
|
The equivalent homogeneous system that is returned is in the form suitable to be used by other functions of the LinearFunctionalSystems package. Namely, it either returns a sequence of a new list of equations and a list of functions to solve for, or a new matrix, vector, variable, and case.
|
|
•
|
This function is part of the LinearFunctionalSystems package, and so it can be used in the form HomogeneousSystem(..) only after executing the command with(LinearFunctionalSystems). However, it can always be accessed through the long form of the command by using the form LinearFunctionalSystems[HomogeneousSystem](..).
|
|
|
|
Examples
|
|
|
>
|
|
|
>
|
![sys := [diff(y1(x), x)-y2(x), diff(y2(x), x)-y3(x)-y4(x), diff(y3(x), x)-y5(x), diff(y4(x), x)-2*y1(x)-2*x*y2(x)-y5(x), diff(y5(x), x)-x^2*y1(x)-2*x*y3(x)-y6(x), diff(y6(x), x)-x^2*y2(x)+2*y3(x)]](/support/helpjp/helpview.aspx?si=8752/file02036/math225.png)
|
|
>
|
![vars := [y1(x), y2(x), y3(x), y4(x), y5(x), y6(x)]](/support/helpjp/helpview.aspx?si=8752/file02036/math229.png)
|
|
>
|
|
|
![[diff(y1(x), x)-y2(x), diff(y2(x), x)-y3(x)-y4(x), diff(y3(x), x)-y5(x), diff(y4(x), x)-2*y1(x)-2*x*y2(x)-y5(x), diff(y5(x), x)-x^2*y1(x)-2*x*y3(x)-y6(x), diff(y6(x), x)-x^2*y2(x)+2*y3(x)], [y1(x), y2(x), y3(x), y4(x), y5(x), y6(x)]](/support/helpjp/helpview.aspx?si=8752/file02036/math236.png)
| (1) |
|
>
|
|
|
| (2) |
|
>
|
![sys := [(diff(y1(x), x))*(-1+x^2)+(6+2*x)*y1(x)+6*y3(x)-6*y4(x) = 0, 2*(diff(y2(x), x))*x*(1+x^2)^2-(4+x)*y2(x)*(1+x^2)^2+4*y3(x)*(1+x^2)^2 = -4*x^2-(4+x)*(1+x^2), 2*(diff(y3(x), x))*(1+x^2)*x-(4+3*x)*(1+x^2)*y2(x)+(4*(1+x^2))*y3(x) = -4-3*x, 2*(diff(y4(x), x))*(1+x^2)*(-1+x^2)*x-((1+x^2)*2)*x*(-7+x^2-2*x)*y1(x)-y2(x)*(4+3*x)*(1+x^2)*(-1+x^2)+2*y3(x)*(6*x+2*x^2-2)*(1+x^2)-(6*y4(x)*(1+x^2)*2)*x = -(4+3*x)*(-1+x^2)]](/support/helpjp/helpview.aspx?si=8752/file02036/math247.png)
|
|
![sys := [(diff(y1(x), x))*(-1+x^2)+(6+2*x)*y1(x)+6*y3(x)-6*y4(x) = 0, 2*(diff(y2(x), x))*x*(1+x^2)^2-(4+x)*y2(x)*(1+x^2)^2+4*y3(x)*(1+x^2)^2 = -4*x^2-(4+x)*(1+x^2), 2*(diff(y3(x), x))*(1+x^2)*x-(4+3*x)*(1+x^2)*y2(x)+4*(1+x^2)*y3(x) = -4-3*x, 2*(diff(y4(x), x))*(1+x^2)*(-1+x^2)*x-2*(1+x^2)*x*(-7+x^2-2*x)*y1(x)-y2(x)*(4+3*x)*(1+x^2)*(-1+x^2)+2*y3(x)*(6*x+2*x^2-2)*(1+x^2)-12*y4(x)*(1+x^2)*x = -(4+3*x)*(-1+x^2)]](/support/helpjp/helpview.aspx?si=8752/file02036/math250.png)
| (3) |
|
>
|
![vars := [y1(x), y2(x), y3(x), y4(x)]](/support/helpjp/helpview.aspx?si=8752/file02036/math254.png)
|
|
>
|
|
|
| (4) |
|
>
|
|
|
| (5) |
|
|
