We deal with the existence of solutions of the Dirichlet problem for sublinear and superlinear partial differential inclusions considered as generalizations of the Euler-Lagrange equation for a certain integral functional without convexity assumption. We develop a duality theory and variational principles for this problem. As a consequence of the duality theory we give a numerical version of the variational principles which enables approximation of the solution for our problem.