In this work we study the degenerate diffusion equation $partial_{t}=x^{alpha}aleft(xright)partial_{x}^{2}+bleft(xright)partial_{x}$ for $left(x,tright)inleft(0,inftyright)^{2}$, equipped with a Cauchy initial data and the Dirichlet boundary condition at $0$. We assume that the order of degeneracy at 0 of the diffusion operator is $alphainleft(0,2right)$, and both $aleft(xright)$ and $bleft(xright)$ are only locally bounded. We adopt a combination of probabilistic approach and analytic method: by analyzing the behaviors of the underlying diffusion process, we give an explicit construction to the fundamental solution $pleft(x,y,tright)$ and prove several properties for $pleft(x,y,tright)$; by conducting a localization procedure, we obtain an approximation for $pleft(x,y,tright)$ for $x,y$ in a neighborhood of 0 and $t$ sufficiently small, where the error estimates only rely on the local bounds of $aleft(xright)$ and $bleft(xright)$ (and their derivatives). There is a rich literature on such a degenerate diffusion in the case of $alpha=1$. Our work extends part of the existing results to cases with more general order of degeneracy, both in the analysis context (e.g., heat kernel estimates on fundamental solutions) and in the probability view (e.g., wellposedness of stochastic differential equations).