
We consider the mixed boundary value problem or Zaremba's problem for the
Laplacian in a bounded Lipschitz domain in R^n. We specify Dirichlet data on
part of the boundary and Neumann data on the remainder of the boundary. We
assume that the boundary between the sets where we specify Dirichlet and
Neumann data is a Lipschitz surface. We require that the Neumann data is in L^p
and the Dirichlet data is in the Sobolev space of functions having one
derivative in L^p for some p near 1. Under these conditions, there is a unique
solution to the mixed problem with the nontangential maximal function of the
gradient of the solution in L^p of the boundary. We also obtain results with
data from Hardy spaces when p=1.

We consider the mixed problem for $L$ the Lam\'e system of elasticity in a
bounded Lipschitz domain $ \Omega\subset\reals ^2$. We suppose that the
boundary is written as the union of two disjoint sets, $\partial\Omega =D\cup
N$. We take traction data from the space $L^p(N)$ and Dirichlet data from a
Sobolev space $ W^{1,p}(D)$ and look for a solution $u$ of $Lu =0$ with the
given boundary conditions. We give a scale invariant condition on $D$ and find
an exponent $ p_0 >1$ so that for $1<p<p_0$, we have a unique solution of this
boundary value problem with the nontangential maximal function of the gradient
of the solution in $L^ p(\partial\Omega)$. We also establish the existence of a
unique solution when the data is taken from Hardy spaces and HardySobolev
spaces with $ p$ in $(p_1,1]$ for some $p_1 <1$.

This paper continues the study of the mixed problem for the Laplacian. We
consider a bounded Lipschitz domain $\Omega\subset \reals^n$, $n\geq2$, with
boundary that is decomposed as $\partial\Omega=D\cup N$, $D$ and $N$ disjoint.
We let $\Lambda$ denote the boundary of $D$ (relative to $\partial\Omega$) and
impose conditions on the dimension and shape of $\Lambda$ and the sets $N$ and
$D$. Under these geometric criteria, we show that there exists $p_0>1$
depending on the domain $\Omega$ such that for $p$ in the interval $(1,p_0)$,
the mixed problem with Neumann data in the space $L^p(N)$ and Dirichlet data in
the Sobolev space $W^ {1,p}(D) $ has a unique solution with the nontangential
maximal function of the gradient of the solution in $L^p(\partial\Omega)$. We
also obtain results for $p=1$ when the Dirichlet and Neumann data comes from
Hardy spaces, and a result when the boundary data comes from weighted Sobolev
spaces.