options nocenter;

data machine;
 input mat temp y @@;
 temp2 = temp*temp;
cards;
1 15 130 1 15 155 1 70  34 1 70  40 1 125  20 1 125  70
1 15  74 1 15 180 1 70  80 1 70  75 1 125  82 1 125  58
2 15 150 2 15 188 2 70 136 2 70 122 2 125  25 2 125  70
2 15 159 2 15 126 2 70 106 2 70 115 2 125  58 2 125  45
3 15 138 3 15 110 3 70 174 3 70 120 3 125  96 3 125 104
3 15 168 3 15 160 3 70 150 3 70 139 3 125  82 3 125  60
;

proc glm;
 class mat;
 model y=temp mat temp2 mat*temp mat*temp2/ solution;
run;