%exm0563_1.m
clear
%（0）從三維坐標初步觀察兩函數圖形相交情況
x=-2:0.05:2;y=x;[X,Y]=meshgrid(x,y);	%產生x-y平面上網點坐標
F1=sin(X-Y);F2=cos(X+Y);
F0=zeros(size(X));
figure(1)
surf(X,Y,F1),
xlabel('x'),ylabel('y'),
view([-31,62]),hold on,
surf(X,Y,F2),surf(X,Y,F0),
shading interp,
hold off  


%（1）在某區域觀察兩函數0等位線的交點情況
clear;
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等位線的存在。
figure(2)
contour(X,Y,F1,v)						%畫F1的三條等位線。
hold on,contour(X,Y,F2,v),hold off	%畫F2的三條等位線。  


%（2）從圖形獲取零點的初始近似值
clc
disp('用鼠標在圖上獲取2個零點初始數據')
[x0,y0]=ginput(2);  					%在圖上取兩個點的坐標
disp([x0,y0])  


%（3）利用fsolve求精確解。以求（0.7926,7843）附近的解為例。
fun='[sin(x(1)-x(2)),cos(x(1)+x(2))]';								%<12>
xy=fsolve(fun,[x0(2),y0(2)])											%<13>  


%（4）檢驗
fxy1=sin(xy(1)-xy(2));
fxy2=cos(xy(1)+xy(2));
disp([fxy1,fxy2])