編寫壹個matlab程序，但是算出了的結果不是壹個數值，而是壹大堆算式……具體程序如下

function output=damnk(a,b)

cita=50;

h=50;

c=20;

gama=18;

fi=35;

n=10;

r=(a^2+b^2)^(1/2);

n1=(a+(r^2-(h-b)^2)^(1/2))/n;

i=1;

Tn=0;

Rn=0;

while(i<=n)

Ri=quad(@(x)(x*tand(cita)+0-b+(r^2-(x-a).^2).^(1/2)).*(r^2-(x-a).^2).^(1/2)/r,(i-1)*n1,i*n1);

Ti=quad(@(x) x*tand(cita)+0-b+(r^2-(x-a).^2).^(1/2),(i-1)*n1,i*n1);

Rn=Rn+Ri;

Tn=Tn+Ti;

i=i+1;

end

jiajiao=acosd(dot([0-a,0-b],[(a+(r^2-(h-b)^2)^(1/2))-a,h-b])/(norm([0-a,0-b])*norm([0-a,h-b])));

l=2*pi*r*(jiajiao/360);

z=(Rn*tand(fi)+c*l/gama)/Tn;

output=z;

但是妳的(r^2-(h-b)^2)^(1/2)這項巨容易成為虛數，註意參數a，b的選擇

其次，妳要算出個數，怎麽能用符號計算呢。。。 壹般直接數值積分~有些可以先符號再賦值計算，這裏貌似不是很好~