Failed to run the code

Hello, I am confronted with the problem that ValueError: shapes (3,200) and (600,1) not aligned: 200 (dim 1) != 600 (dim 0) when execute `gp.fit(theta)`, could you show me your datasets and I want to confirm whether my dataset organization is correct.

By the way, I would like to ask whether the number of units in last layer is equal to the number of output feature.
