Re: [R] help about solving two equations

2010-03-11 Thread Shaoqiong Zhao
anks a lot. Annie - Original Message - From: "Ravi Varadhan" To: "Shaoqiong Zhao" Cc: r-help@r-project.org Sent: Wednesday, March 10, 2010 9:03:25 PM GMT -06:00 US/Canada Central Subject: Re: [R] help about solving two equations Here is how you can solve: fn <- funct

[R] help about solving two equations

2010-03-10 Thread Shaoqiong Zhao
I have two matrix s1 and s2, each of them is 1000*1. and I have two equations: digamma(p)-digamma(p+q)=s1, digamma(q)-digamma(p+q)=s2, and I want to sovle these two equations to get the value of x and y, which are also two 1000*1 matrices. I write a program like this: f <- function(x) { p<- x[1]