xy = read.table("data030501.txt")

x1 = xy[,2]
x2 = xy[,3]
y = xy[,1]

nreg = nls(y~g0*x1^g1*x2^g2, start=list(g0=1, g1 = 1, g2=1))
 
summary(nreg)