Marine animals and micro-machines often use wiggling motion to generate thrust. The wiggling motion can be modeled by a progressive wave where its wavelength describes the flexibility of wiggling animals. In the present study, an immersed boundary method is used to simulate the flows around the wiggling hydrofoil NACA 65-010 at low Reynolds numbers. One can find from the numerical simulations that the thrust generation is largely determined by the wavelength. The thrust coefficients decrease with the increasing wavelength while the propulsive efficiency reaches a maximum at a certain wavelength due to the viscous effects. The thrust generation is associated with two different flow patterns in the wake: the well-known reversed Karman vortex streets and the vortex dipoles. Both are jet-type flows where the thrust coefficients associated with the reversed Karman vortex streets are larger than the ones associated with the vortex diploes.