如何使用R解方程

Question

我想分别将“ dum1”和“ dum2”分别等于5和1.5并求解这两个方程，并获得“ alp”和“ lam”的值。以下是R代码，但错误弹出。 “ fun（x，...）中的错误：缺少参数“ lam”，没有默认值” 有人可以帮我解决吗？

F<-function(alp,lam){
dum1=((1/lam)*(digamma(alp+1)-digamma(1)))
dum2=((1/lam^2)*(psigamma(1)-psigamma(alp+1)))

return(c(dum1-5,dum2-1.5))
}

x<-fsolve(F,c(0,2))