exm04063_1.m

来自「同上」· M 代码 · 共 12 行

x=-2:0.5:2;y=x;[X,Y]=meshgrid(x,y);	%产生x-y平面上网点坐标
F1=sin(X-Y);F2=cos(X+Y);
v=[-0.2, 0, 0.2];	%指定三个等位值，是为了更可靠地判断0等位线的存在。
contour(X,Y,F1,v)						%画F1的三条等位线。
hold on,contour(X,Y,F2,v),hold off	%画F2的三条等位线。  

[x0,y0]=ginput(2);  					%在图上取两个点的坐标
disp([x0,y0]) 

fun='[sin(x(1)-x(2)),cos(x(1)+x(2))]';								%<12>
[xy,f,exit]=fsolve(fun,[x0(2),y0(2)])								%<13>