R fsolve () справка для переменной

Question

У меня есть функция, которая была определена как

ff<-function(x,y){
result<-(2*x*y)+(2*x)-(x^2)-(2*(y^2))
return(result)
}

h=numeric()    
gradient<-function(fun,x0,y0,fx, fy){
    ff2<-replace(of,c("x","y"),c((x0+(fx*h)),(y0+(fy*h))))
  grad<-solve(ff2,0)
  return (grad)
}
gradient(ff, -1,1,6,-6)

Я хотел бы создать новую функцию с именем ff2, заменив все X на (-1 + 6h), а Y на (-1 + -6h) из ff и решите, для Y h это переменная, которая была объявлена ​​как h=numeric() Я сталкиваюсь с двумя проблемами:

1. Когда я пытаюсь заменить X, я получаюошибка x[list] <- values : object of type 'closure' is not subsettable
2. Когда я пытаюсь решить для h с помощью функции решения () r, я получаю ошибку 'a' is 0-diml