求前輩指點(diǎn)迷津,感覺橢流線表達(dá)式有問題,又不知道哪里不對(duì)
3gM{lS}h# 2MC\~"L< 程序如下
6%'{Cq1DE LNg1q1P3 clearclc
dqkkA/1 h=1000;N=100;
4 '+)9&g R=100;for i=2:101
RS:0xN\JN theata(i)=asin(sqrt((i-1)/N));end
O]Hg4">f theata(1)=0;for i=1:101
3LTO+>, |" % r(i)=R*sqrt((i-1)/N);r(i)=sqrt(R^2/i)
6JL
7ut end% r(1)=0;
i3~!ofTb a=zeros(101,1);c=zeros(101,1);
^g0 Ig2' x=zeros(101,1);y=zeros(101,1);
ysa"f+/ a(1)=550;x(1)=-50;
eJwr y(1)=0;
5