hly-neu
2009-05-14, 10:03
key1=[10000 10000];
i=1;
syms s1=0 s2=0 D=0;
e=4-sqrt(3);
n=1;
R=1;
c=pi./6;
h=4;
for r=[0:0.05:10]
for b=[0:pi./20:2*pi]
for a=[0:pi./20:2*pi]
y=(1./(2.*pi.*pi.*10.*10.*(((n+1)./(4.*pi.*pi)).*((((e-0.*sin(b)).*(sin(c))+h.*cos(c)).^n*((h-R.*sin(a).*sin(c))))./(((0.*cos(b)-R.*cos(a)).^2+((0.*sin(b)-e).*cos(c)+(h.*sin(c))-R.*sin(a)).^2+((e-0.*sin(b)).*sin(c)+h.*cos(c)).^2).^2))).^((n+3)./2))).*(((((n+1)./(4.*pi.*pi)).*((((e-r.*sin(b)).*(sin(c))+h.*cos(c)).^n*((h-R.*sin(a).*sin(c))))./(((r.*cos(b)-R.*cos(a)).^2+((r.*sin(b)-e).*cos(c)+(h.*sin(c))-R.*sin(a)).^2+((e-r.*sin(b)).*sin(c)+h.*cos(c)).^2).^2))).^((n+3)./2)));
s=(pi./20).*sum(y);
end
k=(pi./20).*sum(s);
end
D=(k.^2).*r.*0.05;
key1(i)=D;
i=i+1;
end
??? Error using ==> syms at 58
Not a valid variable name.
为什么运行程序是总是提示上面的错误啊？
哪位高手帮忙解决一下呗，谢谢！！！
i=1;
syms s1=0 s2=0 D=0;
e=4-sqrt(3);
n=1;
R=1;
c=pi./6;
h=4;
for r=[0:0.05:10]
for b=[0:pi./20:2*pi]
for a=[0:pi./20:2*pi]
y=(1./(2.*pi.*pi.*10.*10.*(((n+1)./(4.*pi.*pi)).*((((e-0.*sin(b)).*(sin(c))+h.*cos(c)).^n*((h-R.*sin(a).*sin(c))))./(((0.*cos(b)-R.*cos(a)).^2+((0.*sin(b)-e).*cos(c)+(h.*sin(c))-R.*sin(a)).^2+((e-0.*sin(b)).*sin(c)+h.*cos(c)).^2).^2))).^((n+3)./2))).*(((((n+1)./(4.*pi.*pi)).*((((e-r.*sin(b)).*(sin(c))+h.*cos(c)).^n*((h-R.*sin(a).*sin(c))))./(((r.*cos(b)-R.*cos(a)).^2+((r.*sin(b)-e).*cos(c)+(h.*sin(c))-R.*sin(a)).^2+((e-r.*sin(b)).*sin(c)+h.*cos(c)).^2).^2))).^((n+3)./2)));
s=(pi./20).*sum(y);
end
k=(pi./20).*sum(s);
end
D=(k.^2).*r.*0.05;
key1(i)=D;
i=i+1;
end
??? Error using ==> syms at 58
Not a valid variable name.
为什么运行程序是总是提示上面的错误啊？
哪位高手帮忙解决一下呗，谢谢！！！