implement function for daubechies wavelet

hello, please i want to know if i can implement an equation for the daubechies wavelet.I used to use the second derivative of gaussian like this

fe=1000;%sampling frequency
t=[-25/1000:1/fe:25/1000];
l=length(t);
r=floor(l/2);
si=0.0015;
c = 2/(sqrt(3*si*s)*pi^0.25);
k=sqrt(2)*s*si;
k1=s*si;
h=c*(1-(t/k1).^2).*exp(-(t/k).^2);

is it a way to have a similar function for the daubechies wavelet?? thank you