[x, y] = solve(x==y*tan(b*pi/180), (y-1000)/(x-45)==tan(-2*pi/180), y/(x-20)=tan(2*(-2)*pi/180+b*pi/180+pi/2), x, y);
直接 solve(' '), 引号内填你要解的方程的形式。