We design and analyze the first hybridizable discontinuous Galerkin methods for stationary, third-order linear equations in one-space dimension. The methods are defined as discrete versions of characterizations of the exact solution in terms of local problems and transmission conditions. They provide approximations to the exact solution
u
u
and its derivatives
q
:=
u
′
q:=u’
and
p
:=
u
p:=u
which are piecewise polynomials of degree
k
u
k_u
,
k
q
k_q
and
k
p
k_p
, respectively. We consider the methods for which the difference between these polynomial degrees is at most two. We prove that all these methods have superconvergence properties which allows us to prove that their numerical traces converge at the nodes of the partition with order at least
2
k
+
1
2\,k+1
, where
k
k
is the minimum of
k
u
,
k
q
k_u,k_q
, and
k
p
k_p
. This allows us to use an element-by-element post-processing to obtain new approximations for
u
,
q
u, q
and
p
p
converging with order at least
2
k
+
1
2k+1
uniformly. Numerical results validating our error estimates are displayed.