Question

急急……求解MATLAB7.6的一道编程题，不知道哪里错了，运行不出来……麻烦哪位大虾帮帮忙，感激不尽~~

function f=afiberte(x,y)
x=0.0002003:0.017995:6;
y=0.00000001:0.017995:6;
z=(20*pi*x)/3;a=0.02;
cz=1-(100/(x*x));
dz=1-((0.56*x*x)/(x*x-16));
ct=2.25;
dt=1;
c2=0.8;d2=1;
u=(z.*a).*sqrt((cz/ct).*(ct*dt-y.*y));
w=(z.*a).*sqrt(y*y-c2*d2);
f=(cz*(besselj(1,u)))./(u*besselj(0,u))+(c2*besselk(1,w))./(w*(besselk(0,w)));
麻烦指点出哪里错了，运行的结果应该是一副二维的曲线图，谢谢！！！

Answer 1

注意该加点的地方，下面程序运行OK

% not necicerry:
%function f=afiberte(x,y)
x=0.0002003:0.017995:6;
y=0.00000001:0.017995:6;
z=(20*pi*x)/3;
a=0.02;
cz=1-(100./(x.*x)); %vector
dz=1-((0.56*x.*x)./(x.*x-16)); %vector
ct=2.25;
dt=1;
c2=0.8;
d2=1;
u=(z.*a).*sqrt((cz/ct).*(ct*dt-y.*y));
w=(z.*a).*sqrt(y.*y-c2*d2);
f=(cz.*(besselj(1,u)))./(u.*besselj(0,u))+(c2*besselk(1,w))./(w.*(besselk(0,w)));