We consider Weierstrass functions and divisor functions arising from $q$-series. Using these we can obtain new identities for divisor functions. Farkas [3] provided a relation between the sums of divisors satisfying congruence conditions and the sums of numbers of divisors satisfying congruence conditions. In the proof he took logarithmic derivative to theta functions and used the heat equation. In this note, however, we obtain a similar result by differentiating further. For any $n{geq}1$, we have $$k{cdot}{ au}_{2;k,l}(n)=2n{cdot}E_{frac{k-l}{2}}(n;k)+l{cdot}{ au}_{1;k,l}(n)+2k{cdot}{sum_{j=1}^{n-1}}E_{frac{k-1}{2}(j;k){ au}_{1;k,l}(n-j)$$. Finally, we shall give a table for $E_1(N;3)$, ${sigma}(N)$, ${ au}_{1;3,1}(N)$ and ${ au}_{2;3,1}(N)$ ($1{leq}N{leq}50$) and state simulation results for them.