x0=3;
y0=4;
z0=1;                      
x1=0.8090;
y1=0.5878;
z1=0.0036;      
l=x1-x0;
m=y1-y0;
n=z1-z0;
xx=l/sqrt(l^2+m^2+n^2);
yy=m/sqrt(l^2+m^2+n^2);
zz=n/sqrt(l^2+m^2+n^2); 
[color=FF0000]z=220-sqrt((140+sqrt(6400-y^2))^2-x^2); 
mm=diff(z,'x'); 
nn=diff(z,'y');[/color]
qmn=sqrt(1+mm^2+nn^2);
xxn=1/qmn;
yyn=-mm/qmn;
zzn=-nn/qmn; 
n1=1;n2=1.594122;
cos1=xx*xxn+yy*yyn+zz*zzn;
cos2=sqrt(1-n1^2/n2^2*(1-cos1^2));
g=n1*cos1-n2*cos2;
xx1=(n1*xx+g*xxn)/n2
yy1=(n1*yy+g*yyn)/n2
zz1=(n1*zz+g*zzn)/n2

这个程序存在的问题：红色部分，我想求出的是x=0.8090,y=0.5878时x
与y的导数值，只有求出这个值，下面的计算才能继续，怎样把这个计
算融入到这个程序中，请高手帮忙