Question

matlab解不定方程x^+y^2+xy=1,

请问应该怎么改
程序为x=-2:2;
y=-2:2;
[x,y]=meshgrid(x,y);
z=x*x+y*y-x*y;
k=find(z==1);
x(k)',y(k)'

Answer 1

修改一行即可，倒数第三句的乘都换成点乘：
z=x.*x+y.*y-x.*y;
------运行结果有6组整数解：
ans =
-1 -1 0 0 1 1
ans =
-1 0 -1 1 0 1

Answer 2

程序改为如下（把下面指令放入m文件中）：
x=-2:2;y=-2:2;
x=solve('x^2+y^2+x*y=1')
结果为：

x =

(4 - 3*y^2)^(1/2)/2 - y/2
- y/2 - (4 - 3*y^2)^(1/2)/2

追问

谢谢，但是我结果要怎么把x,y解出来啊！要具体数值的！

追答

建个m文件放入下列代码，能得到你想要的
x1=-2:2;
x2=-2:2;
y1=-2:2;
y2=-2:2;
x1=((4 - 3*y1.^2).^0.5)/2 - y1/2;
x2= - y2/2 - (4 - 3*y2.^2).^(1/2)/2;
plot(x1,y1,'r*',x2,y2,'r*')
x1
y1
x2
y2
结果如下：
x1 =
Columns 1 through 3
1.0000 + 1.4142i 1.0000 1.0000
Columns 4 through 5
0 -1.0000 + 1.4142i
y1 =
-2 -1 0 1 2
x2 =
Columns 1 through 3
1.0000 - 1.4142i 0 -1.0000
Columns 4 through 5
-1.0000 -1.0000 - 1.4142i
y2 =
-2 -1 0 1 2
注：x1和y1对应，x2和y2对应。都是你要求的界