The problem
in
with
on
is studied,
where
,
. A pair
is called a Fucík eigenvalue if there exists a
function
that solves the problem. Local existence
results are obtained and a numerical method is used to approximate
some of the Fucík eigenvalues on various domains in
. Based on a variational formulation in
, critical points of the functional
d
subject to constraints
d
and
d
with
are sought. The information obtained numerically is
compared to the known analytical results.