求问一个fsolve解非线性方程组的问题

设计东西要用，按着书上给的例子改完发现没法用，求大佬解答
function t=position(phi23,phi1)
t = [30*cos(phi1)+130*cos(phi23(1))+100*cos(phi23(2))-80;30*sin(phi1)+130*sin(phi23(1))+100*sin(phi23(2))-100]
end
phi1=(0:2/180:2)*pi
phi23=zeros(length(phi1),2);
options=optimset('display','off')
for m=1:length(phi1)
phi23(m,:)=fsolve('position',[1,1],options,phi1(m))
end
目的是求对应phi1从0到2π对应方程组t的解phi23（包含两个变量）