# Determination of unknown coefficients in a nonlocal parabolic problem from neumann type boundary measured data.

1 Introduction

The mathematical model of sludge particles settling in the water treatment plant (settler) is given by the transport-diffusion operator Au := [u.sub.t] - [(D(x)[u.sub.x]).sub.x] + v[(c(x)u(x)).sub.x]. In the case of the residence time of sludge particles in the settler, the model leads to a nonlinear age-dependent transport-diffusion with nonlocal additional condition. Identification of the diffusion coefficient D(x) is considered. This problem is formulated as an identification/optimal control problem, where the boundary flux measurements is assumed to be a control. For the case of constant ("average") velocity, as a characterization of the optimal control problem two necessary conditions are obtained. These conditions permit to reduce the nonlinear coupled two-dimensional problem to the two-point boundary value problem for the second-order nonlinear ordinary differential equation, and then, to a nonlinear equation, with respect to sludge concentration. For the approximate solution of the problem, an inverse polynomial method is used. This method is based on the Lagrange-Hermite interpolation of the solution of the corresponding direct problem and use of the flux information as well as singular points. The point x [member of] (-L, L) will be defined as a singular point for the inverse problem if the first and/or second derivative of the direct problem's solution vanishes at that point x.

The determination of the leading unknown coefficient in ordinary and partial differential equations is one of actual problems in inverse problem theory and practice (see [1-3,7,8,13] and references therein). For ordinary differential equations, the most widely studied inverse problems were formulated as classical Sturm-Liouville problems [4,12], where the unknown coefficient needs to be determined from spectral data. However such a data is not available in many physical problems. The class of inverse problems considered here is based only on boundary measurements [5-11].

2 Problem Formulation

We consider the following identification / optimal control problem: Find a triple <u,D,c>, which satisfies the nonlocal problem

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2.1)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2.2)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2.3)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2.4)

We assume boundary flux measurement(s) as additional conditions:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2.5)

Although the homogeneous boundary conditions (2.3) are not realistic from the physical point of view, we assume here that the density of sludge particles on the boundary is small enough at the boundary x = [+ or -]L. We also assume that the initial data and coefficients are continuous functions, i.e., [phi](x), D(x), c(x) [member of] C[-L, L], and

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2.6)

For given coefficients c, D [member of] C[-L, L], denote by u = u(x, t; D, c) the unique classical solution of the parabolic initial value problem (2.1)-(2.3). The coefficients c(x), D(x) are assumed to be unknown and the nonlocal condition (2.4) and measured data (2.5) can be treated as observations for determination of the coefficients c(x), D(x). For this reason, the problem (2.1) to (2.5) is defined as a nonlocal optimal control or identification problem, and the triple <u, D, c> is called a solution of the identification problem (2.1)-(2.5). In this case, the problem (2.1)-(2.3) is called a direct problem corresponding to inverse problem (2.1)-(2.5).

3 Necessary Conditions for Optimality

Assume that <u, D, c> is a solution of the problem (2.1)-(2.5). Then integrating equation (2.1) on [0, [infinity]) with respect to the variable t > 0, we obtain

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3.1)

By using the nonlocal condition (2.4) and assuming differentiability of the function u(x, t) under the integrals, we may write

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Taking into account the initial condition (2.2), we can rewrite the integro-differential equation (3.1) in the form

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3.2)

Hence we obtain the following

Proposition 3.1. If <u, D, c> is a solution of the identification problem (2.1)-(2.5), then the couple of the functions < c(x), D(x) > satisfy the two-point boundary value problem (3.2) for the second-order nonlinear ordinary differential equation.

Proposition 3.2. The identification problem (2.1)-(2.5) is equivalent to the following inverse problem: Find a pair (D, c), which satisfies the two-point boundary value problem (3.2) and the additional condition(s)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3.3)

The reduced measured data

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

is obtained from nonlocal condition (2.4) and measured data (2.5). The function c = c(x; D) [member of] [C.sup.2] (-L,L) n C 1[--L,L] is the solution to the nonlinear boundary value problem (3.2) for a given coefficient D = D(x). The problem (3.2) is called forward problem corresponding to the inverse problem (3.2)-(3.3).

According to Proposition 3.2, we can now consider the reduced inverse problem (3.2)-(3.3). Rewriting of equation in (3.2) with respect to the unknown function D = D(x) and using the conditions in (3.3) as initial condition, we have the Cauchy problem

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3.4)

When the function c = c(x) is known, the solution of the Cauchy problem (3 4) has the integral representation

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3.5)

It is seen from (3.4) that ill-posedness of the inverse problem is a result of vanishing of the first or/and the second-order derivatives of c = c(x) at one or more points of (-L, L). This suggests to classify the inverse problem (3.2)-(3.3) as follows (see [13], [5]):

(i1) c'(x) [not equal to] 0, c"(x) [not equal to] 0, [for all]x [member of] [-L, L];

(i2) c'([x.sub.0]) = 0, [x.sub.0] [member of] (-L, L); c'(x) [not equal to] 0, [for all]x [not equal to] [x.sub.0]; c"(x) [not equal to] 0, [for all]x [member of] [-L, L];

(i3) c'([x.sub.1]) = c"([x.sub.1]) = 0, [x.sub.1] [member of] (-L, L); c'(x) [not equal to] 0, c"(x) [not equal to] 0, [for all]x [not equal to] [x.sub.1] .

In the case (i1) the inverse problem is well-conditioned. The cases (i2), (i3) correspond to the situations, when the problem is moderately ill-conditioned and severely ill-conditioned, respectively. The points [x.sub.0], [x.sub.1] [member of] (-L, L) will be defined as singular points for the inverse problem. The geometry of singular points is given in Figure 3.1.

[FIGURE 3.1 OMITTED]

Corollary 3.3. The inverse problem (3.2)-(3.3) is ill-conditioned in the above sense (i2) or (i3).

We apply the inverse polynomial method (IPM), given in [5], for moderately and severely ill-conditioned inverse problems. The distinguished feature of this method is that it permits to determine the unknown coefficient in the neighborhood of the singular point also.

4 Inverse Polynomial Method (IPM)

By the inverse polynomial method as introduced in [5], we adopt the view that good approximation of the forward problem solution yields good solution of the corresponding inverse problem. The method consists of two steps: polynomial approximation of the direct problem solution c = c(x) and solution of the Cauchy problem for the unknown coefficient D = D(x).

For each considered ill-conditioned situation, the polynomial [c.sub.p] = [c.sub.p](x) will be defined from the boundary conditions and information about singular points. As a result, we need to construct the polynomial function

[c.sub.p](x) = [[alpha].sub.0] + [[alpha].sub.1]x + [[alpha].sub.2][x.sup.2] + ... + [[alpha].sub.m][x.sup.m], x [member of] [-L, L], (4.1)

where the parameter m > 1 will be defined from the number of additional conditions and the unknown coefficients [[alpha].sub.i], i = [bar.0, m] will be determined from the corresponding linear algebraic equations system:

Ax = B.

The coefficient matrix A = {([a.sub.ij])} and the column vector B = {([b.sub.j])} will be determined separately for the cases (i2) and (i3).

In the second part of the method the following Cauchy problem, obtained from (3.4), will be solved for the first-order differential equation with respect to the unknown coefficient [D.sub.appr.] = [D.sub.appr.](x):

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4.2)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4.3)

It is seen from equation (4.2), if [c'.sub.p](x) [not equal to] 0, [c".sub.p(]x) [not equal to] 0, x [member of] [-L, L], then for a given function [c.sub.p](x), the Cauchy problem (4.2)-(4.3) has the unique solution

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4.4)

In the presented modification of the inverse method, we use the additional information to increase the order of the polynomial approximation for the function c = c(x).

Lemma 4.1 (See [6]). Suppose <D(x),c(x)> satisfies Eq. (3.2) and conditions (3.3). Then the following compatibility condition holds:

[[integral].sup.L.sub.-L] [??](x)dx = [[??].sub.0] - [[??].sub.1] (conservation law).

Proof. This is obtained by integrating equation in (3.2) on the interval [-L,L] and applying homogeneous boundary conditions in (3.2).

Lemma 4.2 (See [6]). Suppose <[D.sub.1](x), [c.sub.1](x)> and <[D.sub.2](x), [c.sub.2](x)> satisfyEq. (3.2) and conditions (3.3). If [x.sub.0] [member of] (-L, L) is a singular point for the functions [c.sub.1](x) and [c.sub.2](x), i.e., if [c'.sub.1]([x.sub.0]) = [c'.sub.2]([x.sub.0]) = 0, then the following transmission conditions hold:

[c.sub.1]([x.sub.a]) = [c.sub.2]([x.sub.a]), [c.sub.1]([x.sub.b]) = [c.sub.2]([x.sub.b]), [x.sub.a] [member of] (-L, [x.sub.0]), [x.sub.b] [member of] ([x.sub.0], L).

Proof. This is obtained by integrating equation in (3.2) on the interval [x, L] and applying the hypothesis of the lemma.

In this case the points [x.sub.a] and [x.sub.b] are called transmission points of the inverse problem (2.1)-(2.5).

5 Ill-Conditioned Situations and Numerical Examples

5.1 The Moderately Ill-Conditioned Situation

If we assume that the situation (i2) holds,

c'([x.sub.0]) = 0, [x.sub.0] [member of] (-L, L); c'(x) [not equal to] 0, [for all]x [not equal to] [x.sub.0]; c"(x) [not equal to] 0, [for all]x [member of] [-L, L],

then the polynomial function [c.sub.p] = [c.sub.p](x), given by (4.1), will be defined from the four conditions

[c.sub.p](-L) = 0, cp(L) = 0, [c'.sub.p]([x.sub.0]) = 0, [[integral].sup.[I.sub.0].sub.-L] [??](x)dx - [[??].sub.0] - [vc.sup.2.sub.p]([x.sub.0]) = 0

As a result, we obtain the Lagrange-Hermite polynomial function [c.sub.p](x) = [a.sub.0] + [a.sub.1]x + [a.sub.2][x.sup.2] + [a.sub.3][x.sup.3], x [member of] [-L,L] (5.1)

of degree 3, and the unknown coefficients [a.sub.i], i = 0, 3 are determined from the system

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (5.2)

where the number [gamma]= [square root of ([[integral].sup.[x.sub.0].sub.-L][phi](x)dx - [[psi].sub.0]/v)] > 0 is obtained from the fourth condition above. The approximate solution [D.sub.appr.] = [D.sub.appr.](x) of the inverse problem (3.2)-(3.3) can be computed from the integral representations (4.4), simultaneously (except the point x = [x.sub.0]). At the singular point, the numerical solution is computed from

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (5.3)

The error estimation of the inverse polynomial method in the considered situation is obtained as follows.

Lemma 5.1. Let c(x) [member of] [C.sup.4][-L, L] and [c.sub.p] = [c.sub.p](x) be the Lagrange-Hermite polynomial given by (5.1). If [x.sub.a] and [x.sub.b] are transmission points given in Lemma 4.2, then the following estimate holds: [for all]x [member of] [-L, L],

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Proof. The points -L, L, [x.sub.a] and [x.sub.b] are zeros of the function F(x) = c(x) - [c.sub.p](x) + [kappa][delta](x), [delta](x) = (x + L)(x - [x.sub.a])(x - [x.sub.b])(x - L). On the other hand, the constant [kappa] can be determined for all x [member of] (-L, L) such that F(x) = 0. Thus we get [kappa] = ([c.sub.p](x) - c(x))/[delta](x). Then by the Rolle theorem, there exists [xi] [member of] (-L,L) such that [F.sup.(4)]([xi]) = 0. The proof is complete.

Corollary 5.2. If the conditions of Lemma 5.1 hold, then the following error formula holds:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Lemma 5.3. Let c(x) [member of] [C.sup.5][0,1] and [c.sub.p] = [c.sub.p](x) be the polynomial given by (5.1). If [x.sub.0] is a moderately singular point of the inverse problem (3.2)-(3.3), then the following estimate holds:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Proof. The points -L, L, [x.sub.a] and [x.sub.b] are zeros of the function W(x) := c(x) - [c.sub.p](x), so by the Rolle theorem there exists [[gamma].sub.i] [member of] (-L,L) such that W'([[gamma].sub.i]) = 0, i = [bar.1, 3]. On the other hand, W'([x.sub.0]) = 0. Introduce the function G(x) = W'(x) + [kappa][delta] (x), [delta](x) = (x - [[gamma].sub.1])(x - [[gamma].sub.2])(x - [[gamma].sub.3])(x - [x.sub.0]), and [kappa] can be determined for all x [member of] (-L, L) such that G(x) = 0. Thus we get [kappa] = - W'(x)/[delta](x). By applying the Rolle theorem again, we find that there exists [xi] [member of] (-L, L) such that [G.sup.(4)]([xi]) = 0. Thus the proof is complete.

Corollary 5.4. If the conditions of Lemma 5.3 hold, then the following error formula holds:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

In a first series of computational experiments, we consider the moderately ill-conditioned situation.

Example 5.5. Let

c(x) = exp(-0.5[x.sub.2]), D(x) = 1 + exp(-[x.sup.2])

be exact solutions of the inverse problem (3.2) for input data

v = 0.5, [??](x) = -D'(x)c'(x) - D(x)c"(x) + 2vc(x)c'(x), x [member of] (-L, L).

Here L = 1.5 and the value [[??].sub.0] is assumed to be synthetic noise free data given on the left boundary of the interval (-L, L), and [x.sub.0] = 0 is a moderately singular point. Figure 5.1 illustrates the exact solutions c(x), D(x) and their polynomial approximations [c.sub.appr.](x), [D.sub.appr.](x) on the the interval [-L,L]. The relative supnorm-errors (in absence of noise in data)

[[epsilon].sub.c] = [[parallel]c - [c.sub.appr.]/c[parallel].sub.[infinity]], [[epsilon].sub.D] = [[parallel]D - [D.sub.appr.]/D[parallel].sub.[infinity]]

are 0.87 and 0.61, respectively. Here we used n = 45 mesh points in the interval [L, L] . Figure 5.2 shows the exact and approximate solutions (on top c(x), on bottom D(x)) for the above test functions with noisy data. The curves -o- in both figures correspond to synthetic noisy Neumann data D(-L)c(L) = [[psi].sub.[gamma]], [[psi].sub.[gamma]] = [psi][+ or -] [gamma][??] with 7 = 0.05. Here we used the same input parameters as before. The relative supnorm-error, defined above, for the noise factor [gamma] = 0.05 are 0.87 and 0.62, respectively. Here we used n = 45 mesh points in the interval [L, L].

5.2 The Severely Ill-Conditioned Situation

Let us assume situation (i3) holds and [x.sub.1] is a singular point of the inverse problem (3.2)-(3.3)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Since [x.sub.1] is a severely singular point of the problem, we get

[phi]([x.sub.1]) = 0.

We can split Cauchy problem (3.4) into two subproblems

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (5.4)

and

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (5.5)

[FIGURE 5.1 OMITTED]

where [delta] > 0 is a tolerance parameter. Due to the conditions above, we can determine the polynomial function [c.sub.p] = [c.sub.p](x), given in (4.1), from the five conditions

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

As a result we obtain the Lagrange-Hermite polynomial function

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (5.6)

of degree 4, and the unknown coefficients a are determined from the system

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (5.7)

[FIGURE 5.2 OMITTED]

where the number [gamma] = [square root of ([[integral].sup.[x.sub.0].sub.-L] [phi]dx - [[psi].sub.0]/v)] is obtained from the fifth condition above. Contrary to previous situation, the approximate solution cp = cp(x) is not defined at the singular point x = [x.sub.1], but this approximate solution can be computed in a sufficient [delta]-neighborhood of the point [x.sub.1] from the Cauchy problems (5.4) and (5.5).

The error estimation of the inverse polynomial method in the considered situation is obtained as follows.

Lemma 5.6. Let c(x) [member of] [C.sup.4][-L, L] and [c.sub.p] = [u.sub.p](x) be the Lagrange-Hermite polynomial of degree 4 given by (5.6). Then for all x [member of] (-L, L),

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

where [x.sub.a] and [x.sub.b] are transmission points which are defined in Lemma 4.2.

Proof. The proof is similar to the proof of Lemma 5.3 and hence is omitted.

Corollary 5.7. If the conditions of Lemma 5.6 hold, then the following error formula holds:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Lemma 5.8. Let c(x) [member of] [C.sup.5][-L, L] and [c.sub.p] = [c.sub.p](x) be the Lagrange-Hermite polynomial of degree 4 given by (5.6). If [x.sub.1] is a severely singular point of the inverse problem (3.2)-(3.3), then the following estimate holds:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Proof. The proof is similar to the proof of Lemma 5.3 and hence is omitted.

Corollary 5.9. If the conditions of Lemma 5.8 hold, then the following error formula holds:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Now we consider an example for the severely ill-conditioned situation.

Example 5.10. Let [x.sub.1] = 0 be a severely singular point and

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

be exact solutions of the inverse problem (3.2)-(3.3) in case (i3) holds for given input data

v = 0.5, [??](x) = -D'(x)c'(x) - D(x)c"(x) + 2vc(x)c'(x), x [member of] (-L, L).

Here L = 1 .5 and the value [??]0 = D(-L)c'(L) is assumed to be synthetic noise free data given on the left boundary of the interval (-L, L). Figure 5.3 illustrates the exact solution c(x) and its polynomial approximations [c.sub.appr.](x) of the direct problem and exact solution D(x) and the approximate solution [D.sub.appr.](x) of the inverse problem on the interval [-L,[x.sub.1] - [delta]) [union] ([x.sub.1] + [delta], L], [delta] = 0.25. The relative supnorm-errors (in absence of noise in data) are 0.28 and 2.4 and 3.3, respectively. Here we used n = 36 mesh points in the interval [-L, [x.sub.1] [delta]) [union] ([x.sub.1] + [delta], L]. Consider now this example in the case of noisy data. The curve in Figure 5.4 corresponds to the synthetic noisy Neumann data D(-L)c'(-L) = [[??].sub.[gamma]], [[??].sub.[gamma]] = [[??].sub.0] [+ or -] [gamma][[??].sub.0]. The relative supnorm-errors of the functions <c(x); [D.sup.1](x), [D.sup.2](x)> are 1.04, 1.03 and 2.34 for the noise factor 7 = 0.05, respectively. Here the approximate solutions [D.sup.1](x), [D.sup.2](x) are solutions of the Cauchy problems (5.4) and (5.5), respectively. We used n = 36 mesh points in the interval [-L, [x.sub.1] - [delta]) [union] ([x.sub.1] + [delta], L] for [delta] = 0.25.

[FIGURE 5.3 OMITTED]

6 Conclusions

We have considered inverse coefficient problem (2.1)-(2.5) with nonlocal condition for parabolic equation and reduced it to an inverse problem in ordinary differential equations. We have applied inverse polynomial method for approximate solution of the inverse problem which is based on the Lagrange-Hermite interpolation of the solution of the corresponding direct problem and use of the information about singular points.

[FIGURE 5.4 OMITTED]

An efficiency and applicability of the method is demonstrated on various numerical examples with noisy free and noisy data.

On the other hand, the Dirichlet type conditions (2.3) and Neumann type additional conditions (2.5) are interchangeable in the inverse problem (2.1)-(2.5), i.e., the conditions (2.5) can be taken as boundary conditions of the direct problem, as the condition (2.3) are taken as additional condition of the inverse problem. So we can get different new inverse problems as follows:

Problem 1 (NN-DD)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Additional conditions:

u(-L,t) = [g.sub.0](t), u(L,t) = [g.sub.1](t), [for all]t [member of] [R.sub.+].

Problem 2 (DN-DN)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Additional conditions:

-D(-L)[u.sub.x](-L,t) = [f.sub.0](t), u(L, t) = [g.sub.1](t), [for all]t [member of] [R.sub.+].

Acknowledgments

The author wishes to thank Alemdar Hasanov for the problem formulation and fruitful suggestions on the topic of polynomial approximation. The author also thanks A. Okay (Celebi for carefully reading the paper and the valuable suggestions that improved the paper.

References

[1] J. R. Cannon and P. C. DuChateau, Determining of unknown coefficients in a nonlinear heat conduction problem, SIAMJ. Appl. Math. 24 (1973), 298-314.

[2] M. Choulli and P. Stefano, Reconstruction of the coefficients of the stationary transport equation from boundary measurements, Inverse Problems 12 (1996) L19-L23.

[3] P. C. DuChateau, Introduction to inverse problems in partial differential equations for engineers, physicists and mathematicians, in: Parameter Identification and Inverse Problems in Hydrology, Geology and Ecology (J. Gottlieb and P DuChateau eds.), Netherland, Kluver Academic Publishers (1996) 3-38.

[4] I. M. Gelfand and I. M. Levitan, On the determination of a differential equation from its spectrum, Trans. Amer. Math. Soc. 1 (1955) 253-304.

[5] A. Hasanov, An inverse polynomial method for the identification of the leading coefficient in the Sturm-Liouville operator from boundary measurements, Appl. Math. Comput. 140 (2003) 501-515.

[6] A. Hasanov, Error analysis of a multisingular inverse coefficient problem for the Sturm-Liouville operator based on boundary measurement, Appl. Math. Comput. 150(2004) 493-524.

[7] A. Hasanov, P. C. Duchateau and B. Pekta?, Inverse and ill-posed Problems 14(5) (2006) 435-463.

[8] A. Hasanov and B. Pekta?, Simulation of ill-conditioned situations of inverse coefficient problem for the Sturm-Liouville operator based on boundary measurements, Math. Comput. Simulation 61 (2001) 47-52.

[9] A. Hasanov and Z. Seyidmamedov, Determination of leading coefficients in Sturm-Liouville operator from boundary measurements. I. A stripping algorithm, Appl. Math. Comput. 125 (2001) 1-21.

[10] A. Hasanov and Z. Seyidmamedov, Determination of leading coefficients in Sturm-Liouville operator from boundary measurements. II. Unicity and an engineering approach, Appl. Math. Comput. 125 (2001) 23-34.

[11] A. Hasanov and T. S. Shores, Solution of an inverse coefficient problem for an ordinary differential equation Appl. Anal. 67 (1997) 11-20.

[12] H. Hochstedt, The inverse Sturm-Liouville problem, Comm. Pure Appl. Math. 26 (1973)715-729.

[13] G. R. Richter, An inverse problem for the steady state diffusion equation, SIAM J. Appl. Math 41 (1981) 210-221.

Burhan Pektas

Izmir University

Department of Mathematics and Computer Sciences

35350 Uckuyular, Izmir, Turkey

burhan.pektas@izmir.edu.tr

burhanps@gmail.com

Received April 9, 2010; Accepted May 6, 2010

Communicated by A. Okay Celebi

The mathematical model of sludge particles settling in the water treatment plant (settler) is given by the transport-diffusion operator Au := [u.sub.t] - [(D(x)[u.sub.x]).sub.x] + v[(c(x)u(x)).sub.x]. In the case of the residence time of sludge particles in the settler, the model leads to a nonlinear age-dependent transport-diffusion with nonlocal additional condition. Identification of the diffusion coefficient D(x) is considered. This problem is formulated as an identification/optimal control problem, where the boundary flux measurements is assumed to be a control. For the case of constant ("average") velocity, as a characterization of the optimal control problem two necessary conditions are obtained. These conditions permit to reduce the nonlinear coupled two-dimensional problem to the two-point boundary value problem for the second-order nonlinear ordinary differential equation, and then, to a nonlinear equation, with respect to sludge concentration. For the approximate solution of the problem, an inverse polynomial method is used. This method is based on the Lagrange-Hermite interpolation of the solution of the corresponding direct problem and use of the flux information as well as singular points. The point x [member of] (-L, L) will be defined as a singular point for the inverse problem if the first and/or second derivative of the direct problem's solution vanishes at that point x.

The determination of the leading unknown coefficient in ordinary and partial differential equations is one of actual problems in inverse problem theory and practice (see [1-3,7,8,13] and references therein). For ordinary differential equations, the most widely studied inverse problems were formulated as classical Sturm-Liouville problems [4,12], where the unknown coefficient needs to be determined from spectral data. However such a data is not available in many physical problems. The class of inverse problems considered here is based only on boundary measurements [5-11].

2 Problem Formulation

We consider the following identification / optimal control problem: Find a triple <u,D,c>, which satisfies the nonlocal problem

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2.1)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2.2)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2.3)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2.4)

We assume boundary flux measurement(s) as additional conditions:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2.5)

Although the homogeneous boundary conditions (2.3) are not realistic from the physical point of view, we assume here that the density of sludge particles on the boundary is small enough at the boundary x = [+ or -]L. We also assume that the initial data and coefficients are continuous functions, i.e., [phi](x), D(x), c(x) [member of] C[-L, L], and

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2.6)

For given coefficients c, D [member of] C[-L, L], denote by u = u(x, t; D, c) the unique classical solution of the parabolic initial value problem (2.1)-(2.3). The coefficients c(x), D(x) are assumed to be unknown and the nonlocal condition (2.4) and measured data (2.5) can be treated as observations for determination of the coefficients c(x), D(x). For this reason, the problem (2.1) to (2.5) is defined as a nonlocal optimal control or identification problem, and the triple <u, D, c> is called a solution of the identification problem (2.1)-(2.5). In this case, the problem (2.1)-(2.3) is called a direct problem corresponding to inverse problem (2.1)-(2.5).

3 Necessary Conditions for Optimality

Assume that <u, D, c> is a solution of the problem (2.1)-(2.5). Then integrating equation (2.1) on [0, [infinity]) with respect to the variable t > 0, we obtain

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3.1)

By using the nonlocal condition (2.4) and assuming differentiability of the function u(x, t) under the integrals, we may write

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Taking into account the initial condition (2.2), we can rewrite the integro-differential equation (3.1) in the form

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3.2)

Hence we obtain the following

Proposition 3.1. If <u, D, c> is a solution of the identification problem (2.1)-(2.5), then the couple of the functions < c(x), D(x) > satisfy the two-point boundary value problem (3.2) for the second-order nonlinear ordinary differential equation.

Proposition 3.2. The identification problem (2.1)-(2.5) is equivalent to the following inverse problem: Find a pair (D, c), which satisfies the two-point boundary value problem (3.2) and the additional condition(s)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3.3)

The reduced measured data

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

is obtained from nonlocal condition (2.4) and measured data (2.5). The function c = c(x; D) [member of] [C.sup.2] (-L,L) n C 1[--L,L] is the solution to the nonlinear boundary value problem (3.2) for a given coefficient D = D(x). The problem (3.2) is called forward problem corresponding to the inverse problem (3.2)-(3.3).

According to Proposition 3.2, we can now consider the reduced inverse problem (3.2)-(3.3). Rewriting of equation in (3.2) with respect to the unknown function D = D(x) and using the conditions in (3.3) as initial condition, we have the Cauchy problem

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3.4)

When the function c = c(x) is known, the solution of the Cauchy problem (3 4) has the integral representation

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3.5)

It is seen from (3.4) that ill-posedness of the inverse problem is a result of vanishing of the first or/and the second-order derivatives of c = c(x) at one or more points of (-L, L). This suggests to classify the inverse problem (3.2)-(3.3) as follows (see [13], [5]):

(i1) c'(x) [not equal to] 0, c"(x) [not equal to] 0, [for all]x [member of] [-L, L];

(i2) c'([x.sub.0]) = 0, [x.sub.0] [member of] (-L, L); c'(x) [not equal to] 0, [for all]x [not equal to] [x.sub.0]; c"(x) [not equal to] 0, [for all]x [member of] [-L, L];

(i3) c'([x.sub.1]) = c"([x.sub.1]) = 0, [x.sub.1] [member of] (-L, L); c'(x) [not equal to] 0, c"(x) [not equal to] 0, [for all]x [not equal to] [x.sub.1] .

In the case (i1) the inverse problem is well-conditioned. The cases (i2), (i3) correspond to the situations, when the problem is moderately ill-conditioned and severely ill-conditioned, respectively. The points [x.sub.0], [x.sub.1] [member of] (-L, L) will be defined as singular points for the inverse problem. The geometry of singular points is given in Figure 3.1.

[FIGURE 3.1 OMITTED]

Corollary 3.3. The inverse problem (3.2)-(3.3) is ill-conditioned in the above sense (i2) or (i3).

We apply the inverse polynomial method (IPM), given in [5], for moderately and severely ill-conditioned inverse problems. The distinguished feature of this method is that it permits to determine the unknown coefficient in the neighborhood of the singular point also.

4 Inverse Polynomial Method (IPM)

By the inverse polynomial method as introduced in [5], we adopt the view that good approximation of the forward problem solution yields good solution of the corresponding inverse problem. The method consists of two steps: polynomial approximation of the direct problem solution c = c(x) and solution of the Cauchy problem for the unknown coefficient D = D(x).

For each considered ill-conditioned situation, the polynomial [c.sub.p] = [c.sub.p](x) will be defined from the boundary conditions and information about singular points. As a result, we need to construct the polynomial function

[c.sub.p](x) = [[alpha].sub.0] + [[alpha].sub.1]x + [[alpha].sub.2][x.sup.2] + ... + [[alpha].sub.m][x.sup.m], x [member of] [-L, L], (4.1)

where the parameter m > 1 will be defined from the number of additional conditions and the unknown coefficients [[alpha].sub.i], i = [bar.0, m] will be determined from the corresponding linear algebraic equations system:

Ax = B.

The coefficient matrix A = {([a.sub.ij])} and the column vector B = {([b.sub.j])} will be determined separately for the cases (i2) and (i3).

In the second part of the method the following Cauchy problem, obtained from (3.4), will be solved for the first-order differential equation with respect to the unknown coefficient [D.sub.appr.] = [D.sub.appr.](x):

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4.2)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4.3)

It is seen from equation (4.2), if [c'.sub.p](x) [not equal to] 0, [c".sub.p(]x) [not equal to] 0, x [member of] [-L, L], then for a given function [c.sub.p](x), the Cauchy problem (4.2)-(4.3) has the unique solution

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4.4)

In the presented modification of the inverse method, we use the additional information to increase the order of the polynomial approximation for the function c = c(x).

Lemma 4.1 (See [6]). Suppose <D(x),c(x)> satisfies Eq. (3.2) and conditions (3.3). Then the following compatibility condition holds:

[[integral].sup.L.sub.-L] [??](x)dx = [[??].sub.0] - [[??].sub.1] (conservation law).

Proof. This is obtained by integrating equation in (3.2) on the interval [-L,L] and applying homogeneous boundary conditions in (3.2).

Lemma 4.2 (See [6]). Suppose <[D.sub.1](x), [c.sub.1](x)> and <[D.sub.2](x), [c.sub.2](x)> satisfyEq. (3.2) and conditions (3.3). If [x.sub.0] [member of] (-L, L) is a singular point for the functions [c.sub.1](x) and [c.sub.2](x), i.e., if [c'.sub.1]([x.sub.0]) = [c'.sub.2]([x.sub.0]) = 0, then the following transmission conditions hold:

[c.sub.1]([x.sub.a]) = [c.sub.2]([x.sub.a]), [c.sub.1]([x.sub.b]) = [c.sub.2]([x.sub.b]), [x.sub.a] [member of] (-L, [x.sub.0]), [x.sub.b] [member of] ([x.sub.0], L).

Proof. This is obtained by integrating equation in (3.2) on the interval [x, L] and applying the hypothesis of the lemma.

In this case the points [x.sub.a] and [x.sub.b] are called transmission points of the inverse problem (2.1)-(2.5).

5 Ill-Conditioned Situations and Numerical Examples

5.1 The Moderately Ill-Conditioned Situation

If we assume that the situation (i2) holds,

c'([x.sub.0]) = 0, [x.sub.0] [member of] (-L, L); c'(x) [not equal to] 0, [for all]x [not equal to] [x.sub.0]; c"(x) [not equal to] 0, [for all]x [member of] [-L, L],

then the polynomial function [c.sub.p] = [c.sub.p](x), given by (4.1), will be defined from the four conditions

[c.sub.p](-L) = 0, cp(L) = 0, [c'.sub.p]([x.sub.0]) = 0, [[integral].sup.[I.sub.0].sub.-L] [??](x)dx - [[??].sub.0] - [vc.sup.2.sub.p]([x.sub.0]) = 0

As a result, we obtain the Lagrange-Hermite polynomial function [c.sub.p](x) = [a.sub.0] + [a.sub.1]x + [a.sub.2][x.sup.2] + [a.sub.3][x.sup.3], x [member of] [-L,L] (5.1)

of degree 3, and the unknown coefficients [a.sub.i], i = 0, 3 are determined from the system

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (5.2)

where the number [gamma]= [square root of ([[integral].sup.[x.sub.0].sub.-L][phi](x)dx - [[psi].sub.0]/v)] > 0 is obtained from the fourth condition above. The approximate solution [D.sub.appr.] = [D.sub.appr.](x) of the inverse problem (3.2)-(3.3) can be computed from the integral representations (4.4), simultaneously (except the point x = [x.sub.0]). At the singular point, the numerical solution is computed from

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (5.3)

The error estimation of the inverse polynomial method in the considered situation is obtained as follows.

Lemma 5.1. Let c(x) [member of] [C.sup.4][-L, L] and [c.sub.p] = [c.sub.p](x) be the Lagrange-Hermite polynomial given by (5.1). If [x.sub.a] and [x.sub.b] are transmission points given in Lemma 4.2, then the following estimate holds: [for all]x [member of] [-L, L],

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Proof. The points -L, L, [x.sub.a] and [x.sub.b] are zeros of the function F(x) = c(x) - [c.sub.p](x) + [kappa][delta](x), [delta](x) = (x + L)(x - [x.sub.a])(x - [x.sub.b])(x - L). On the other hand, the constant [kappa] can be determined for all x [member of] (-L, L) such that F(x) = 0. Thus we get [kappa] = ([c.sub.p](x) - c(x))/[delta](x). Then by the Rolle theorem, there exists [xi] [member of] (-L,L) such that [F.sup.(4)]([xi]) = 0. The proof is complete.

Corollary 5.2. If the conditions of Lemma 5.1 hold, then the following error formula holds:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Lemma 5.3. Let c(x) [member of] [C.sup.5][0,1] and [c.sub.p] = [c.sub.p](x) be the polynomial given by (5.1). If [x.sub.0] is a moderately singular point of the inverse problem (3.2)-(3.3), then the following estimate holds:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Proof. The points -L, L, [x.sub.a] and [x.sub.b] are zeros of the function W(x) := c(x) - [c.sub.p](x), so by the Rolle theorem there exists [[gamma].sub.i] [member of] (-L,L) such that W'([[gamma].sub.i]) = 0, i = [bar.1, 3]. On the other hand, W'([x.sub.0]) = 0. Introduce the function G(x) = W'(x) + [kappa][delta] (x), [delta](x) = (x - [[gamma].sub.1])(x - [[gamma].sub.2])(x - [[gamma].sub.3])(x - [x.sub.0]), and [kappa] can be determined for all x [member of] (-L, L) such that G(x) = 0. Thus we get [kappa] = - W'(x)/[delta](x). By applying the Rolle theorem again, we find that there exists [xi] [member of] (-L, L) such that [G.sup.(4)]([xi]) = 0. Thus the proof is complete.

Corollary 5.4. If the conditions of Lemma 5.3 hold, then the following error formula holds:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

In a first series of computational experiments, we consider the moderately ill-conditioned situation.

Example 5.5. Let

c(x) = exp(-0.5[x.sub.2]), D(x) = 1 + exp(-[x.sup.2])

be exact solutions of the inverse problem (3.2) for input data

v = 0.5, [??](x) = -D'(x)c'(x) - D(x)c"(x) + 2vc(x)c'(x), x [member of] (-L, L).

Here L = 1.5 and the value [[??].sub.0] is assumed to be synthetic noise free data given on the left boundary of the interval (-L, L), and [x.sub.0] = 0 is a moderately singular point. Figure 5.1 illustrates the exact solutions c(x), D(x) and their polynomial approximations [c.sub.appr.](x), [D.sub.appr.](x) on the the interval [-L,L]. The relative supnorm-errors (in absence of noise in data)

[[epsilon].sub.c] = [[parallel]c - [c.sub.appr.]/c[parallel].sub.[infinity]], [[epsilon].sub.D] = [[parallel]D - [D.sub.appr.]/D[parallel].sub.[infinity]]

are 0.87 and 0.61, respectively. Here we used n = 45 mesh points in the interval [L, L] . Figure 5.2 shows the exact and approximate solutions (on top c(x), on bottom D(x)) for the above test functions with noisy data. The curves -o- in both figures correspond to synthetic noisy Neumann data D(-L)c(L) = [[psi].sub.[gamma]], [[psi].sub.[gamma]] = [psi][+ or -] [gamma][??] with 7 = 0.05. Here we used the same input parameters as before. The relative supnorm-error, defined above, for the noise factor [gamma] = 0.05 are 0.87 and 0.62, respectively. Here we used n = 45 mesh points in the interval [L, L].

5.2 The Severely Ill-Conditioned Situation

Let us assume situation (i3) holds and [x.sub.1] is a singular point of the inverse problem (3.2)-(3.3)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Since [x.sub.1] is a severely singular point of the problem, we get

[phi]([x.sub.1]) = 0.

We can split Cauchy problem (3.4) into two subproblems

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (5.4)

and

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (5.5)

[FIGURE 5.1 OMITTED]

where [delta] > 0 is a tolerance parameter. Due to the conditions above, we can determine the polynomial function [c.sub.p] = [c.sub.p](x), given in (4.1), from the five conditions

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

As a result we obtain the Lagrange-Hermite polynomial function

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (5.6)

of degree 4, and the unknown coefficients a are determined from the system

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (5.7)

[FIGURE 5.2 OMITTED]

where the number [gamma] = [square root of ([[integral].sup.[x.sub.0].sub.-L] [phi]dx - [[psi].sub.0]/v)] is obtained from the fifth condition above. Contrary to previous situation, the approximate solution cp = cp(x) is not defined at the singular point x = [x.sub.1], but this approximate solution can be computed in a sufficient [delta]-neighborhood of the point [x.sub.1] from the Cauchy problems (5.4) and (5.5).

The error estimation of the inverse polynomial method in the considered situation is obtained as follows.

Lemma 5.6. Let c(x) [member of] [C.sup.4][-L, L] and [c.sub.p] = [u.sub.p](x) be the Lagrange-Hermite polynomial of degree 4 given by (5.6). Then for all x [member of] (-L, L),

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

where [x.sub.a] and [x.sub.b] are transmission points which are defined in Lemma 4.2.

Proof. The proof is similar to the proof of Lemma 5.3 and hence is omitted.

Corollary 5.7. If the conditions of Lemma 5.6 hold, then the following error formula holds:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Lemma 5.8. Let c(x) [member of] [C.sup.5][-L, L] and [c.sub.p] = [c.sub.p](x) be the Lagrange-Hermite polynomial of degree 4 given by (5.6). If [x.sub.1] is a severely singular point of the inverse problem (3.2)-(3.3), then the following estimate holds:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Proof. The proof is similar to the proof of Lemma 5.3 and hence is omitted.

Corollary 5.9. If the conditions of Lemma 5.8 hold, then the following error formula holds:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Now we consider an example for the severely ill-conditioned situation.

Example 5.10. Let [x.sub.1] = 0 be a severely singular point and

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

be exact solutions of the inverse problem (3.2)-(3.3) in case (i3) holds for given input data

v = 0.5, [??](x) = -D'(x)c'(x) - D(x)c"(x) + 2vc(x)c'(x), x [member of] (-L, L).

Here L = 1 .5 and the value [??]0 = D(-L)c'(L) is assumed to be synthetic noise free data given on the left boundary of the interval (-L, L). Figure 5.3 illustrates the exact solution c(x) and its polynomial approximations [c.sub.appr.](x) of the direct problem and exact solution D(x) and the approximate solution [D.sub.appr.](x) of the inverse problem on the interval [-L,[x.sub.1] - [delta]) [union] ([x.sub.1] + [delta], L], [delta] = 0.25. The relative supnorm-errors (in absence of noise in data) are 0.28 and 2.4 and 3.3, respectively. Here we used n = 36 mesh points in the interval [-L, [x.sub.1] [delta]) [union] ([x.sub.1] + [delta], L]. Consider now this example in the case of noisy data. The curve in Figure 5.4 corresponds to the synthetic noisy Neumann data D(-L)c'(-L) = [[??].sub.[gamma]], [[??].sub.[gamma]] = [[??].sub.0] [+ or -] [gamma][[??].sub.0]. The relative supnorm-errors of the functions <c(x); [D.sup.1](x), [D.sup.2](x)> are 1.04, 1.03 and 2.34 for the noise factor 7 = 0.05, respectively. Here the approximate solutions [D.sup.1](x), [D.sup.2](x) are solutions of the Cauchy problems (5.4) and (5.5), respectively. We used n = 36 mesh points in the interval [-L, [x.sub.1] - [delta]) [union] ([x.sub.1] + [delta], L] for [delta] = 0.25.

[FIGURE 5.3 OMITTED]

6 Conclusions

We have considered inverse coefficient problem (2.1)-(2.5) with nonlocal condition for parabolic equation and reduced it to an inverse problem in ordinary differential equations. We have applied inverse polynomial method for approximate solution of the inverse problem which is based on the Lagrange-Hermite interpolation of the solution of the corresponding direct problem and use of the information about singular points.

[FIGURE 5.4 OMITTED]

An efficiency and applicability of the method is demonstrated on various numerical examples with noisy free and noisy data.

On the other hand, the Dirichlet type conditions (2.3) and Neumann type additional conditions (2.5) are interchangeable in the inverse problem (2.1)-(2.5), i.e., the conditions (2.5) can be taken as boundary conditions of the direct problem, as the condition (2.3) are taken as additional condition of the inverse problem. So we can get different new inverse problems as follows:

Problem 1 (NN-DD)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Additional conditions:

u(-L,t) = [g.sub.0](t), u(L,t) = [g.sub.1](t), [for all]t [member of] [R.sub.+].

Problem 2 (DN-DN)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Additional conditions:

-D(-L)[u.sub.x](-L,t) = [f.sub.0](t), u(L, t) = [g.sub.1](t), [for all]t [member of] [R.sub.+].

Acknowledgments

The author wishes to thank Alemdar Hasanov for the problem formulation and fruitful suggestions on the topic of polynomial approximation. The author also thanks A. Okay (Celebi for carefully reading the paper and the valuable suggestions that improved the paper.

References

[1] J. R. Cannon and P. C. DuChateau, Determining of unknown coefficients in a nonlinear heat conduction problem, SIAMJ. Appl. Math. 24 (1973), 298-314.

[2] M. Choulli and P. Stefano, Reconstruction of the coefficients of the stationary transport equation from boundary measurements, Inverse Problems 12 (1996) L19-L23.

[3] P. C. DuChateau, Introduction to inverse problems in partial differential equations for engineers, physicists and mathematicians, in: Parameter Identification and Inverse Problems in Hydrology, Geology and Ecology (J. Gottlieb and P DuChateau eds.), Netherland, Kluver Academic Publishers (1996) 3-38.

[4] I. M. Gelfand and I. M. Levitan, On the determination of a differential equation from its spectrum, Trans. Amer. Math. Soc. 1 (1955) 253-304.

[5] A. Hasanov, An inverse polynomial method for the identification of the leading coefficient in the Sturm-Liouville operator from boundary measurements, Appl. Math. Comput. 140 (2003) 501-515.

[6] A. Hasanov, Error analysis of a multisingular inverse coefficient problem for the Sturm-Liouville operator based on boundary measurement, Appl. Math. Comput. 150(2004) 493-524.

[7] A. Hasanov, P. C. Duchateau and B. Pekta?, Inverse and ill-posed Problems 14(5) (2006) 435-463.

[8] A. Hasanov and B. Pekta?, Simulation of ill-conditioned situations of inverse coefficient problem for the Sturm-Liouville operator based on boundary measurements, Math. Comput. Simulation 61 (2001) 47-52.

[9] A. Hasanov and Z. Seyidmamedov, Determination of leading coefficients in Sturm-Liouville operator from boundary measurements. I. A stripping algorithm, Appl. Math. Comput. 125 (2001) 1-21.

[10] A. Hasanov and Z. Seyidmamedov, Determination of leading coefficients in Sturm-Liouville operator from boundary measurements. II. Unicity and an engineering approach, Appl. Math. Comput. 125 (2001) 23-34.

[11] A. Hasanov and T. S. Shores, Solution of an inverse coefficient problem for an ordinary differential equation Appl. Anal. 67 (1997) 11-20.

[12] H. Hochstedt, The inverse Sturm-Liouville problem, Comm. Pure Appl. Math. 26 (1973)715-729.

[13] G. R. Richter, An inverse problem for the steady state diffusion equation, SIAM J. Appl. Math 41 (1981) 210-221.

Burhan Pektas

Izmir University

Department of Mathematics and Computer Sciences

35350 Uckuyular, Izmir, Turkey

burhan.pektas@izmir.edu.tr

burhanps@gmail.com

Received April 9, 2010; Accepted May 6, 2010

Communicated by A. Okay Celebi

Printer friendly Cite/link Email Feedback | |

Author: | Pektas, Burhan |
---|---|

Publication: | Advances in Dynamical Systems and Applications |

Article Type: | Report |

Geographic Code: | 7TURK |

Date: | Dec 1, 2010 |

Words: | 4110 |

Previous Article: | Leighton--Coles--Wintner type oscillation criteria for half-linear impulsive differential equations. |

Next Article: | Recent developments of Lyapunov-type inequalities. |

Topics: |