Inverse problems in general are important in modern-day mathematics as they appear in many situations in physics, engineering, biology and medicine. This paper deals with inverse problems pertaining to a special type of second order difference equation. A comprehensive introduction to difference equations can be found, for example, in  and  , amongst others. In particular, inverse problems for Sturm-Liouville difference equations with Dirichlet boundary conditions have been considered recently by Bohner and Koyunbakan in  where they show that the specification of the eigenvalues and weights uniquely determines the potential. In addition, they also prove that if the potential is symmetric, then it is uniquely determined by the eigenvalues only―this result can also be found in  where it is proved using different methods.
This paper is a sequel to  where the following second order difference equation was considered
where represents the weights associated with a potential function,. The boundary conditions, which were imposed respectively at the initial and terminal end- points, had the general form
Given the weights and the eigenvalues for the above boundary value problem with and, a unique reconstruction of the potential was obtained, see  for details. This can be considered as a generalization of the results obtained in  in that more general boundary conditions are considered.
We now investigate the following inverse problem. Given a spectrum for a boundary value problem of the form (1), (2) and (3), with and, together with the potentials, we prove that provided the number of eigenvalues exceeds, where is the number of unit intervals, it is possible to uniquely reconstruct the weights. If the number of eigenvalues is less than we will require a second spectrum corresponding to a boundary value problem with exactly the same equation and boundary conditions of the same form in order to obtain a unique solution for.
The paper has the following structure. The proof of the above inverse problem is done inductively beginning with the cases of (and) which is considered in Section 2. The inductive process continues in Section 3 with the cases corresponding to. The main result of this paper is a generalization of the results obtained in Sections 2 and 3 and is given in Section 4, in particular Theorem 4.1.
As we do not have experimental data for the two examples presented in Section 3, the eigenvalues used are obtained by first solving the “forward” problem. Consequently, the theoretical results for certain of the inverse problems are then verified using these eigenvalues.
An important result concerning the number of eigenvalues associated with a par- ticular boundary value problem was proved in  and will be used throughout this paper.
Theorem 1.1. Consider the boundary value problem given by Equation (1) for
together with boundary conditions (2) and (3) where. This pro- blem has
1) eigenvalues if,
2) eigenvalues if and,
3) eigenvalues if.
(Note that the number of unit intervals considered is).
2. Results for m = 1, 2
In this section we investigate how to reconstruct the difference boundary value problem using a given spectrum or spectra. That is, how does the given spectrum/spectra, to- gether with the potential function, provide us with a method of computing the associated weights,.
Consider (1) with boundary conditions
The cases of or correspond respectively to or in (2) and (3).
In certain instances it is necessary to consider a second boundary value problem in order to obtain unique results. The second problem will be given by (1) with boundary conditions of the form
The case is a “special” case in that if there is only one eigenvalue then two spectra are required in order to reconstruct uniquely. If there are two or more eigenvalues, then only one spectrum is needed in order to obtain a unique solution for and. The results for can be split into four cases (according to the number of given eigenvalues) and these mirror the results obtained in  except that the values of are known and are uniquely recon- structed.
The cases for follow along the lines of those in  if the weights and potentials are interchanged and with the proviso that if the number of eigenvalues is greater than or equal to 3 then only one spectrum is required but if the number of eigenvalues is less than 3 then two spectra are necessary to produce unique results. In this case, we are given and and can uniquely reconstruct, and.
3. Main Results for m = 3
As the inequalities for differ from one another, we now extend these results to include. There is an increased number of variables that need to be reconstructed which increases the technicality of the proofs. That is, given the spectrum/spectra and the potentials, , can we uniquely reconstruct the weights,
Again, we will split the inverse problem for into four cases depending on the coefficients in the boundary conditions (4), (5) and/or (6), (7). We begin with the cases where we require two spectra to reconstruct the weights uniquely.
Theorem 3.1. Consider the boundary value problems (1), (4), (5) and (1), (6), (7) where
1), , and, ,;
2), , and, ,.
The boundary value problems corresponding to 1) and 2) have three eigenvalues, say, , for (1), (4), (5) and, , for (1), (6), (7). Given, , , , , together with, and, the weights, , and may be uniquely reconstructed in either case.
Proof. From Theorem 1.1, for both cases 1) and 2), it is clear that the boundary value problems each have three eigenvalues.
1) Assume, , and, , in (4), (5) and (6), (7). Then at, (1) together with (4) gives
Next, evaluating (1) at and using (8) we obtain
Also, at we have
Then applying (5) gives the equation
which on simplification yields a polynomial in of degree three namely,
This can be rewritten in the form
The eigencondition is given by
For (1), (6) and (7) we obtain the same third order equation as (12) with replaced by and replaced by their tilde (~) versions.
In addition, the eigencondition in this case is given by
Equating coefficients of, , in (12) and (13) results in
and similarly considering the coefficients of, and yields
Solving the six simultaneous equations gives, , and uni- quely.
2) Suppose that, , and, ,. Then (8), (9) and (10) hold and (11) becomes
This can be rewritten as a cubic polynomial i.e. in the form (12) where, and can be found using Mathematica or similar programmes.
Again for (1), (6) and (7) we obtain the same cubic polynomial as above with replaced by and replaced by their tilde (~) versions. The eigenconditions are given by (13) and (14) thus, Equations (15) and (16) hold. Hence, we can again solve the six simultaneous equations uniquely for, , and.
Theorem 3.2. Consider the boundary value problem (1), (4) and (5) where
1), , , and, , ,;
2), , , and, , ,;
3), , and, ,;
The boundary value problem corresponding to any of the four cases above has four eigenvalues, say, , and for (1), (4), (5) and, , and for (1), (6), (7). The weights, , and may be uniquely reconstructed in any of these cases assuming we are given the eigenvalues together with, and.
Proof. This uses the procedure outlined in Theorem 3.1 above. It is similar to the proof of Theorem 3.2 in  but with increased dependence on as a polynomial of degree four is obtained. The eigencondition is also a fourth order polynomial namely,
and similarly for. Equating relevant coefficients of the powers of and and solving the resulting equations simultaneously, yields unique values for, , and.
Theorem 3.3. Consider the boundary value problem (1), (4) and (5) where
1), , ,;
Given the five eigenvalues, associated with the boundary value problems above, as well as the values, and, it is possible to recon- struct, , and uniquely.
Proof. In all three cases, starting with in (1) and imposing the boundary con- ditions where appropriate, yields a fifth order polynomial in, that is
The associated eigencondition is also a fifth order polynomial given by
By equating relevant coefficients of powers of, unique solutions for, , and result from the simultaneous solution of the resulting equations.
Theorem 3.4. Assume that we have the boundary value problem (1), (4) and (5) with, and. Given the six eigenvalues, associated with this boundary value problem, together with the values, and, a uni- que reconstruction of, , and is possible.
Proof. As per usual we start the evaluation of (1) at together with (4). Following procedures outlined above in (3.2) and  , we obtain an equation which reduces to a sixth order polynomial in. Equating coefficients of this polynomial with those of the eigencondition provides a unique solution for, , and.
To illustrate part (3) of Theorem 3.3, suppose that, ,
, and are the given eigenvalues corre- sponding to the boundary value problem
Clearly, it is seen that,. Using these values together with, , and the method outlined in the theorems, one finds, , , and as given in (18). Equating these with the coefficients of (19) and simultaneously solving the five equations yields, ,
Assume that we are given eigenvalues, ,
, and corresponding to the boundary value problem
To illustrate Theorem 3.2(3), suppose also that, , , and are the eigenvalues corresponding to the boundary value problem
Note that the two boundary value problems above are of the same form i.e. they have exactly the same equation and their boundary conditions are of the same type. In the first problem, while
in the second problem and, ,
for both problems. Using these values routine calculations give, , , , , , and. Simultaneously solving the eight equations yields
, , and.
4. General Case m = r
As mentioned previously, the case of is a special case in that if the number of eigenvalues is less than, then two spectra are required. However, if the number of eigenvalues is greater than or equal to then only one spectrum is required. Using the previous two sections and with, we inductively obtain results for the general case where. That is, given the spectrum/spectra and the potentials,
the theorem below provides an algorithmic method enabling one to uniquely reconstruct the weights,.
Theorem 4.1. For, assume that we are given the spectrum of the boundary value problem (1) with boundary conditions (4) and (5) where the number of eigenvalues k is calculated according to Theorem 1.1. In addition, suppose that the potentials are known,
1) if then the weights can be uniquely reconstructed;
2) if and we are given a second spectrum corresponding to the boundary value problem (1) with boundary conditions (6) and (7) then the weights can be uniquely reconstructed.
Proof. Follows as in (  , Theorem 4.1) where the k simultaneous equations in 1) and 2k simultaneous equations in 2) are now solved to find unique values for
Remark: 1) It is not possible for the number of eigenvalues of (1), (4) and (5) to be less than m as this would imply a Dirichlet boundary condition at, meaning that the interval shrinks and we wish to maintain the number of unit intervals as.
2) It should be noted that because there are more weights than potentials, two spectra are required more often than in  in order to obtain a unique reconstruction of the weights.
We thank the Editor and the referee for their comments. Research of S. Currie is supported by NRF grant no. IFR2011040100017. This support is greatly appreciated.
 Bohner, M. and Koyunbakan, H. (2016) Inverse problems for Sturm-Liouville Difference Equations. Filomat, 30, 1297-1304.
 Hald, O.H. (1976) Discrete Inverse Sturm-Liouville Problems. I. Uniqueness for Symmetric Potentials. Numerische Mathematik, 27, 249-256.