Hi how do I go about fitting a bunch of 3D points to quadric polynomial surface similar to photo? can use some least square process with a vandermonde matrix??

Fitted Values 100 80 60 40 6 20 4 2 O 6 -2 4 2 0 -2 -4 -6 -6