
Boundary value problems for quasilinear differential equations appear in many mathematical models nowadays. Solvability and dependence of the number of solutions on parameters are often asked by the people from praxis. A prototype of such a problem is the Dirichlet problem for the socalled pLaplacian. In this talk we consider pLaplacian in dimension one. The solvability of * equation omitted here can be approached using bifurcation theory (bifurcation from the infinity). In this context we arrive at a linearized problem, solvability of which is determined by linear fredholm alternative. In order to use this approach, we need to know eigenvalues and eigenfunctions to the weighted linear problem * equation omitted here The eigenvalue problem is solved by a shooting method, starting with ψ'(0) = 1 and seeking for λ such that ψ(1) = 0. The weight functions  φ_{k}' ^(p  2) and  φ_{k} ()^(p 2) are singular or degenerated at some points. Another problem is that the function φ_{k}, which appears in the weight functions, is not known explicitly. Moreover asymptotic behavior of  φ_{k}' ^(p  2) and  φ_{k} ()^(p 2) close to points x where either φ_{k}' (x) = 0 or φ_{k} (x) = 0 is known. For that reason we treat Eq.(2) using the event locator  when the values of φ_{k}' or φ_{k} are sufficiently far from zero, we compute these weights by numerical ODE integration, when they are close to zero we use an asymptotic formula. The initial value problem corresponding to Eq.(2) and the initial value problem Eq.(3) are solved simultaneously, and we require that * equation omitted here The talk concentrates on underlying ideas of general interest and their implementation. Theoretical reasons for why to compute these eigenfunctions will be presented on a poster by Benedikt, Čepička, Girg, and Takáč.

