Region of existence of multiple solutions for a class of Robin type four-point BVPs

This article aims to prove the existence of a solution and compute the region of existence for a class of four-point nonlinear boundary value problems (NLBVPs) defined as \[\begin{gathered} -u''(x)=\psi(x,u,u'), \quad x\in (0,1),\\ u'(0)=\lambda_{1}u(\xi), \quad u'(1)=\lambda_{2} u(\eta),\end{gathered}\] where \(I=[0,1]\), \(0\lt\xi\leq\eta\lt 1\) and \(\lambda_1,\lambda_2\gt 0\). The nonlinear source term \(\psi\in C(I\times\mathbb{R}^2,\mathbb{R})\) is one sided Lipschitz in \(u\) with Lipschitz constant \(L_1\) and Lipschitz in \(u'\) such that \(|\psi(x,u,u')-\psi(x,u,v')|\leq L_2(x)|u'-v'|\). We develop monotone iterative technique (MI-technique) in both well ordered and reverse ordered cases. We prove maximum, anti-maximum principle under certain assumptions and use it to show the monotonic behaviour of the sequences of upper-lower solutions. The sufficient conditions are derived for the existence of solution and verified for two examples. The above NLBVPs is linearised using Newton's quasilinearization method which involves a parameter \(k\) equivalent to \(\max_u\frac{\partial \psi}{\partial u}\). We compute the range of \(k\) for which iterative sequences are convergent.