resp_surf function

Response surface model