Multiple pins with phased shift

Hello everybody,

Is it possible to generate 40 KHz square waves from multiple digital pins with Arduino (Uno or MEGA)? while be able to adjust a phase shift to each pin? (each pin generates wave at the predefined delay or phase )

Regards