matlab程序，运行不出来，其中x1与y1都是1000*1000的矩阵

syms x1 y1;
for o=1:1001;
for p=1:1001;
if o+p==1002
x1(o,p)=0;
else
[x1,y1]=solve('-1*x1*sin(2*pi*cf*(-501+o-501+p))+y1*cos(2*pi*cf*(-501+o-501+p))-skk(o,p)=0','x1^2+y1^2-s0(o,p)^2+s1(o,p)^2=0');
S2=vpa(x1,4);
%精度运算S3=vpa(y1,4);
end
end
end

skk，s0，s1都是1000*1000的矩阵并且已知。其中o+p=1002是因为，我怕sin（2*pi*cf*(-501+o-501+p)）在求解中会变为0，可能会产生分母为0 的情况，所以写了这个，不知道是否需要，如果不需要的话又怎么写？

syms x1 y1;
% cf=1; ssk, s1,....
for o=1:1001;
for p=1:1001;
if o+p==1002
sx1(o,p)=0;
else
[sx1,sy1]=solve(-1*x1*sin(2*pi*cf*(-501+o-501+p))+y1*cos(2*pi*cf*(-501+o-501+p))-skk(o,p)==0,x1^2+y1^2-s0(o,p)^2+s1(o,p)^2==0);
S2=vpa(x1,4);
%精度运算S3=vpa(y1,4);
end
end
end

这个主要是set的用法 它是功能 设置对象的属性。 用法 set(h,'propertyname',propertyvalue,…) 用属性值'propertyvalue'设置关于用参量h 标志的对象（一个或多个）的属性名'propertyname'（一个或多个）。h 可以 为一句柄的向量。在这种情形下，命令set 可以设置所有对象的属性值。 set(h,a) 用指定的属性值设置由h 标志的对象的属性。其中a 为一结构数组，该 结构数组的域名为对象的属性名，域名值为相应属性名的属性值。 set(h,pn,pv…) 对由h 指定的所有对象中指定的细胞数组属性名pn 设置为相应 的细胞数组属性值pv。 set(h,pn,) 对于每m 个图形对象设置n 个属性值，其中 m=length(h)，n 为包含属性名的细胞数组pn 中包含的属性名个数。即允许 用户对每一对象的指定的属性设置不同的属性值。 a= set(h) 返回句柄h 中允许用户设置的属性名与可能的属性值。输出参量a 为一 结构数组，其域名为对象的属性名，域名值为相应的属性名对应的属性值。 若没有指定输出参量a，则系统自动将信息显示于屏幕，h 必须为标量。 a= set(0,'factory') 返回那些用户可以设置缺省值的所有对象的属性名，同时显示 可能的属性值，输出参量a 为一结构数组，其域名为对象的属性名，域名值 为相应的属性名对应的属性值，若没有指定输出参量a，则系统自动将信息 显示于屏幕。 a= set(0,'factoryobjecttypepropertyname')返回指定根对象（0）类型中指定的属性 名objecttypepropertyname 的所有可能的属性值。输入参量是由固定的关键 字factory、对象类型（如axes）与属性名（如position 等）组成。 a= set(h,'default') 返回由h 标记的对象上缺省设置的值，其中h 必须是标量。 a= set(h,'defaultobjecttypepropertyname') 返回指定对象h 的类型中指定的属性 名objecttypepropertyname 的所有可能的属性值。输入参量是由固定的关键 字factory、对象类型（如axes）与属性名（如position 等）组成。