East-West J of Mathematics: Vol 22, No (2020) pp 30-51 https://doi.org/10.36853/ewjm.2020.22.01/03 EXISTENCE RESULTS AND ITERATIVE METHOD FOR SOLVING SYSTEMS OF BEAMS EQUATIONS Dang Quang A∗ and Ngo T Kim Quy† ∗ Center for Informatics and Computing Vietnam Academy of Science and Technology (VAST) 18 Hoang Quoc Viet, Cau Giay, Hanoi, Viet Nam e-mail: dangquanga@cic.vast.vn † Posts and Telecommunications Institute of Technology Hanoi, Viet Nam e-mail: quyntk@ptit.edu.vn Abstract In this paper, we propose a method for investigating the solvability and iterative solution of coupled beams equations with fully nonlinear terms Differently from other authors, we reduce the problem to an operator equation for the right-hand side functions The advantage of the proposed method is that it does not require any Nagumo-type conditions for the nonlinear terms Some examples, where exact solution of the problem are known or not, demonstrate the effectiveness of the obtained theoretical results Introduction In the beginning of the 2017 Minh´ os and Coxe [7] for the first time considered the fully fourth order coupled system u(4)(t) v(4) (t) = f(t, u(t), u (t), u (t), u (t), v(t), v (t), v (t), v (t)), = h(t, u(t), u (t), u (t), u (t), v(t), v (t), v (t), v (t)) (1) Key words: Fourth order coupled system; Fixed point theorems; Existence and uniqueness of solution; Iterative method 2010 AMS Mathematics Classification: 34B15, 65L10, 65L20 30 Dang Quang A and Ngo T Kim Quy 31 with the boundary conditions u(0) = u (0) = u (0) = u (1) = 0, v(0) = v (0) = v (0) = v (1) = (2) They gave sufficient conditions for the solvability of the system by using the lower and upper solutions method and the Schauder fixed point theorem The proof of this result is very cumbersome and complicated It requires Nagumotype conditions for the sum of the functions f and h Furthermore, it contained some errors due to the use of non-correct definition of the norm of the space C × C The necessary corrections are made in the Corrigendum in [8] Motivated by the above fact, in this paper we study the system (1)-(2) by another method, namely by reducing it to an operator equation for the pair of nonlinear terms but not for the pair of the functions to be sought (u, v) Without any Nagumo-type conditions and under some easily verified conditions we establish the existence and uniqueness of a solution of the system (1)-(2) Besides, we also prove the property of sign preserving of the solution and the convergence of an iterative method for finding the solution Some examples, where exact solutions of the problem are known or not, demonstrate the effectiveness of the obtained theoretical results The method used here is a further development of the method proposed in our recent works [1, 2, 3, 4] Note that some particular cases of the system (1) were studied before, namely, in [5, 10] the authors considered the equations containing only even order derivatives associated with the boundary conditions different from (2) Under very complicated conditions, by using a fixed point index theorem on cones, the authors obtained the existence of positive solutions But it should be emphasized that the obtained results are of pure theoretical character because no examples of existing solutions are shown The paper is organized as follows In Section we consider the existence and uniqueness of a solution of the problem (1)-(2) and its sign preservation In Section we study an iterative method for solving the problem, where the convergence of iterations is proved Section is devoted to some examples for demonstrating the applicability and efficiency of our approach Finally, Section is Conclusion Existence of a solution To investigate the problem (1)-(2), for u, v ∈ C [0, 1] we set ϕ(t) = f(t, u(t), u (t), u (t), u (t), v(t), v (t), v (t), v (t)), ψ(t) = h(t, u(t), u (t), u (t), u (t), v(t), v (t), v (t), v (t)), T w = (ϕ, ψ) (3) 32 Existence results and iterative method for Then the problem becomes u(4) (t) = ϕ(t), v(4) (t) = ψ(t), < t < 1, 0 such that the functions f(t, U, V ) and h(t, U, V ) are continuous and max{|f(t, U, V )|, |h(t, U, V )|} ≤ M (20) for any (t, U, V ) ∈ DM Then, the problem (1)-(2) has a solution satisfying the estimates M , 12 M |v (t)| ≤ , 12 M , 24 M , |v(t)| ≤ 24 |u (t)| ≤ |u(t)| ≤ M , M |v (t)| ≤ , |u (t)| ≤ M , M |v (t)| ≤ |u (t)| ≤ for any ≤ t ≤ Proof Since the problem (4) is reduced to the operator equation (18), the theorem will be proved if we show that this operator equation has a solution For this purpose, first we show that the operator T defined by (19) maps the closed ball B[0, M ] into itself Let w be an element in B[O, M ] Then, from (8)-(10) it is easy to obtain u ≤ ϕ , 24 u ≤ ϕ , 12 v ≤ ψ , 24 v ≤ ψ 12 (21) 35 Dang Quang A and Ngo T Kim Quy For estimating u and u we notice that the solutions of the problem (12), (14) can be represented in the form u2 (t) = v2 (t) = G2 (t, s)ϕ(s)ds, G2 (t, s)ψ(s)ds, (22) where G2 (t, s) is the Green function −s + st, st − t, G2 (t, s) = ≤ s ≤ t ≤ 1, ≤ t ≤ s ≤ It is easy to verify that max 0≤t≤1 |G2 (t, s)|ds = (23) Therefore, taking into account (22) we have = u2 ≤ u ϕ , v = v2 ≤ ψ (24) Now, rewrite (22) in the form u2 (t) = v2 (t) = t (−s + st)ϕ(s)ds + t (st − t)ϕ(s)ds, t (−s + st)ψ(s)ds + t (st − t)ψ(s)ds (25) From here we obtain u3 (t) = u2 (t) = v3 (t) = v2 (t) = t 1 sϕ(s)ds + t (s − 1)ϕ(s)ds = G3 (t, s)ϕ(s)ds, t 1 sψ(s)ds + t (s − 1)ψ(s)ds = G3 (t, s)ψ(s)ds, (26) where G3 (t, s) is the function continuous in the square [0, 1]2 except for the line s=t s, ≤ s < t ≤ 1, G3 (t, s) = s − 1, ≤ t < s ≤ Hence, M M = v3 ≤ ϕ , v ψ (27) 2 Taking into account (21), (24), (27) and w = max{ ϕ , ψ } ≤ M we have u = u3 ≤ M , 24 M , v ≤ 24 u ≤ M , 12 M ≤ , 12 M , M ≤ , M , M ≤ u1 ≤ u2 ≤ u3 ≤ v1 v2 v3 (28) 36 Existence results and iterative method for Therefore, (t, U, V ) ∈ DM for t ∈ [0, 1] From the definition of T by (19), (17) and the condition (20), we have T w ∈ B[0, M ], i.e., the operator T maps the ball B[0, M ] into itself Next, we prove that the operator T is a compact one in the space F Providing the subscript ϕ for u and ψ for v in the formulas (8), (9), (22) and (26) we have G(t, s)ϕ(s)ds, G(t, s)ψ(s)ds, (29) uϕ (t) = vψ (t) = G1 (t, s)ϕ(s)ds, G (t, s)ψ(s)ds, (30) uϕ (t) = vψ (t) = G2 (t, s)ϕ(s)ds, G2 (t, s)ψ(s)ds (31) uϕ (t) = vψ (t) = G3 (t, s)ϕ(s)ds, G3 (t, s)ψ(s)ds (32) uϕ (t) = vψ (t) = According to [6, Sec 31] the integral operators in (29)-(32) which put each pair of functions (ϕ, ψ) ∈ F in correspondence to the pairs of functions (uϕ , vψ ), (uϕ , vψ ), (uϕ , vψ ), (uϕ , vψ ), are compact operators Therefore, in view of the continuity of the functions f(t, U, V ), h(t, U, V ) it is easy to deduce that the operator T defined by (19) is compact operator in the space F Thus, T is a compact operator from the closed ball B[0, M ]) into itself By the Schauder Fixed Point Theorem [9] the operator equation (18) has a solution The theorem is proved ✷ We now denote ++ DM = {(t, u, u1, u2 , u3 , v, v1 , v2 , v3 )} , where M M , ≤ u1 ≤ , ≤ u2 ≤ 24 12 M M , ≤ v1 ≤ , ≤ v2 ≤ 0≤v≤ 24 12 ≤ t ≤ 1, ≤ u ≤ and M M , |u3 | ≤ , M M , |v3 | ≤ , −− SM = {w ∈ F | − M ≤ ϕ(t) ≤ 0, −M ≤ ψ(t) ≤ 0} −− ++ +− −+ −+ +− Similarly, we introduce the notations DM , SM , DM , SM , DM , SM as follows −− = {(t, u, u1, u2 , u3, v, v1 , v2 , v3 )} , DM Dang Quang A and Ngo T Kim Quy 37 where M M M M ≤ u ≤ 0, − ≤ u1 ≤ 0, − ≤ u2 ≤ 0, |u3 | ≤ , 24 12 M M M M ≤ v ≤ 0, − ≤ v1 ≤ 0, − ≤ v2 ≤ 0, |v3 | ≤ , − 24 12 ≤ t ≤ 1, − and ++ = {w ∈ F | ≤ ϕ(t) ≤ M, ≤ ψ(t) ≤ M } ; SM +− = {(t, u, u1 , u2, u3 , v, v1 , v2 , v3 )} , DM where M M M M , ≤ u1 ≤ , ≤ u2 ≤ , |u3 | ≤ , 24 12 M M M M ≤ v ≤ 0, − ≤ v1 ≤ 0, − ≤ v2 ≤ 0, |v3 | ≤ , − 24 12 ≤ t ≤ 1, ≤ u ≤ and −+ SM = {w ∈ F | − M ≤ ϕ(t) ≤ 0, ≤ ψ(t) ≤ M } ; −+ = {(t, u, u1 , u2, u3 , v, v1 , v2 , v3 )} , DM where ≤ t ≤ 1, − and M M M M ≤ u ≤ 0, − ≤ u1 ≤ 0, − ≤ u2 ≤ 0, |u3 | ≤ , 24 12 M M M M , ≤ v1 ≤ , ≤ v2 ≤ , |v3 | ≤ , 0≤v≤ 24 12 +− = {w ∈ F | ≤ ϕ(t) ≤ M, −M ≤ ψ(t) ≤ 0} SM Now consider some particular cases of Theorem Theorem (Positivity or negativity of solution) ++ (i) Suppose that in DM the functions f, h are continuous and −M ≤ f(t, U, V ) ≤ 0, −M ≤ h(t, U, V ) ≤ (33) Then, the problem (1)-(2) has a solution (u(t), v(t)) with the properties u(t) ≥ 0, u (t) ≥ 0, u (t) ≥ 0, v(t) ≥ 0, v (t) ≥ 0, v (t) ≥ −− (ii) Suppose that in DM the functions f, h are continuous and ≤ f(t, U, V ) ≤ M, ≤ h(t, U, V ) ≤ M (34) 38 Existence results and iterative method for Then, the problem (1)-(2) has a solution (u(t), v(t)) with the properties u(t) ≤ 0, u (t) ≤ 0, u (t) ≤ 0, v(t) ≤ 0, v (t) ≤ 0, v (t) ≤ +− (iii) Suppose that in DM the functions f, h are continuous and −M ≤ f(t, U, V ) ≤ 0, ≤ h(t, U, V ) ≤ M (35) Then, the problem (1)-(2) has a solution (u(t), v(t)) with the properties u(t) ≥ 0, u (t) ≥ 0, u (t) ≥ 0, v(t) ≤ 0, v (t) ≤ 0, v (t) ≤ −+ (iv) Suppose that in DM the functions f, h are continuous and ≤ f(t, U, V ) ≤ M, −M ≤ h(t, U, V ) ≤ (36) Then, the problem (1)-(2) has a solution (u(t), v(t)) with the properties u(t) ≤ 0, u (t) ≤ 0, u (t) ≤ 0, v(t) ≥ 0, v (t) ≥ 0, v (t) ≥ Proof The existence of a solution (u(t), v(t)) of the problem in the case (i) is proved in a similar way as in Theorem 1, where instead of DM and B[0, M ] ++ −− there stand DM and SM The sign of u(t), v(t) and their derivatives are deduced from the representations (8), (9), (22) if taking into account the sign of ϕ(s), ψ(s) and that G(t, s), G1 (t, s), G2 (t, s) are nonpositive functions The proof of the cases (ii), (iii) and (iv) is similar to that of (i), where instead ++ −− −− ++ +− −+ , SM ) there stand the pairs (DM , SM ), (DM , SM ) and of the pair (DM −+ +− (DM , SM ), respectively Now we denote ui1 = (ui ) , ui2 = (ui ) , ui3 = (ui ) ; v1i = (vi ) , v2i = (vi ) , v3i = (vi ) ; U i = (ui , ui1 , ui2 , ui3 ), ϕi = f(t, U i , V i ), V i = (vi , v1i , v2i , v3i ); ψi = h(t, U i , V i ); (i = 1.2) Theorem (Uniqueness of solution) Suppose that there exist numbers ci , di ≥ (i = 0, , 7) such that |f(t, U , V ) − f(t, U , V )| ≤ c0 |u2 − u1 | + c1 |u21 − u11 | + c2 |u22 − u12 | + c3 |u23 − u13 | + c4 |v − v | + c5 |v12 − v11 | + c6 |v22 − v21 | + c7 |v32 (37) − v31 |, |h(t, U , V ) − h(t, U , V )| ≤ d0 |u2 − u1 | + d1 |u21 − u11 | + d2 |u22 − u12 | + d3 |u23 − u13 | + d4 |v − v | + d5 |v12 − v11 | + d6 |v22 − v21 | + d7 |v32 − v31 |, (38) 39 Dang Quang A and Ngo T Kim Quy for any (t, U, V ), (t, U i , V i ) ∈ [0, 1] × R8 (i = 1, 2), and q := max{q1 , q2 } < (39) with c0 + c4 c1 + c5 c2 + c6 c3 + c7 + + + , 24 12 d0 + d4 d1 + d5 d2 + d6 d3 + d7 + + + q2 := 24 12 Then the solution of the problem (1)-(2) is unique if it exists q1 := Proof Suppose the problem has two solutions (u1 (t), v1 (t)) and (u2 (t), v2 (t)) Due to the estimates (28) we have ϕ2 − ϕ1 , 24 ≤ ϕ2 − ϕ1 , ≤ ψ2 − ψ1 , 24 ≤ ψ2 − ψ1 , ϕ2 − ϕ1 , 12 u23 − u13 ≤ ϕ2 − ϕ1 v12 − v11 ≤ ψ2 − ψ1 , 12 v32 − v31 ≤ ψ2 − ψ1 u2 − u1 ≤ u22 − u12 v −v v22 − v21 u21 − u11 ≤ (40) From (37), (38) and the above estimates we have w2 − w1 = max { f(t, U2 , V2 ) − f(t, U1 , V1 ) , h(t, U2, V2 ) − h(t, U1 , V1 ) } ≤ max {q1 max{ ϕ2 − ϕ1 , ψ2 − ψ1 }, q2 max{ ϕ2 − ϕ1 , ψ2 − ψ1 }} ≤ q w2 − w1 (41) with c0 + c4 c1 + c5 c2 + c6 c3 + c7 + + + , 24 12 d0 + d4 d1 + d5 d2 + d6 d3 + d7 + + + , q2 = 24 12 q = max{q1 , q2 } q1 = Since q < the inequality (41) occurs only in the case w2 = w1 This implies u2 = u1 and v2 = v1 Thus, the theorem is proved Theorem Assume that there exist numbers M, ci , di ≥ (i = 0, , 7) such that max{|f(t, U, V )|, |h(t, U, V )|} ≤ M, (42) |f(t, U , V ) − f(t, U , V )| ≤ c0 |u2 − u1 | + c1 |u21 − u11 | + c2 |u22 − u12 | + c3 |u23 − u13 | + c4 |v − v | + c5 |v12 − v11 | + c6 |v22 − v21 | + c7 |v32 − v31 |, (43) ... establish the existence and uniqueness of a solution of the system (1)-(2) Besides, we also prove the property of sign preserving of the solution and the convergence of an iterative method for finding... consider the existence and uniqueness of a solution of the problem (1)-(2) and its sign preservation In Section we study an iterative method for solving the problem, where the convergence of iterations... u1 ≤ u2 ≤ u3 ≤ v1 v2 v3 (28) 36 Existence results and iterative method for Therefore, (t, U, V ) ∈ DM for t ∈ [0, 1] From the definition of T by (19), (17) and the condition (20), we have T