Is this a question about how to do this in Mathematica? If yes, show your code attempts. If not, post this elsewhere.