How do perform polynomial features and linear regression on 3 dimensional data (x,y,z)?

Question

I need help performing polynomial features on 3 dimensional data and performing linear regression to create a line of best fit on the 3 dimensional polynomial.

I have a r