You need to include the constant function $1$ or it won't work. With this I want to show inductively that
$$\chi_{[k2^{-n},(k+1)2^{-n}]}$$
lies in the span for all $0≤k<2^n$. This implies $\chi_{[a,b]}$ lies in the span (as an element of $L^2$, ie differing from that function on a null-set) whenever $a,b$ are of the form $2^n k$ for some $n,k$. That implies (since such points are dense in $[0,1]$) that the characteristic function of any sub-interval lies in the closure of the span, but the closure of that is $L^2[0,1]$ so the span of the wavelets was dense.
Start with $n=1$,
$$\psi_{0,0}+1=2\chi_{[0,1/2]}\qquad 1-\psi_{0,0}=2\chi_{[1/2,1]}$$
Now $\psi_{n,k}=\chi_{[k2^{-n},(k+1/2)2^{-n}]}-\chi_{[(k+1/2)2^{-n},(k+1)2^{-n}]}$. From induction we are supposing $\chi_{[k2^{-n},(k+1)2^{-n}]}$ is in the span, so it follows that
$$\chi_{[k2^{-n},(k+1/2)2^{-n}]}-\chi_{[(k+1/2)2^{-n},(k+1)2^{-n}]}+\chi_{[k2^{-n},(k+1)2^{-n}]}=2\chi_{[k2^{-n},(k+1/2)2^{-n}]}=2\chi_{[(2k)2^{-(n+1)},(2k+1)2^{-(n+1)}]}$$
Lies in the closure of the span for any $0≤2k<2^{n+1}$. The odd case follows by considering
$$\chi_{[k2^{-n},(k+1)2^{-n}]}-(\chi_{[k2^{-n},(k+1/2)2^{-n}]}-\chi_{[(k+1/2)2^{-n},(k+1)2^{-n}]})=2\chi_{[(k+1/2)2^{-n},(k+1)2^{-n}]}=2\chi_{[(2k+1)2^{-(n+1)},(2k+2)2^{-(n+1)}]}$$