%數(shù)學(xué)實(shí)驗(yàn)中的最速落徑問題
%A,B是不在一條鉛垂線上的兩點(diǎn)，在鏈接A,B兩點(diǎn)的所有光滑曲線中，找出一條曲線，使得初速度為零的質(zhì)點(diǎn)，在重力的作用下，自A點(diǎn)下滑到B點(diǎn)所需的時(shí)間最短

clear all
%賦初值
H=2.0;n=100;h=H/n;
g=9.8;f=1.0;
%求常數(shù)c
a=0;
b=2/(sqrt(2*g*(n-1)/n*H)+sqrt(2*g*H));
c=(a+b)/2;
i=1;
while abs(f)>1e-10
   s=0;
   for j=1:n
      v=(sqrt(2*g*j*h)+sqrt(2*g*(j-1)*h))/2;
      s=s+v/sqrt(1.0-c^2*v^2);
   end
   f=c-pi/(h*s);
   if f>0
      b=c;
   else a=c;
   end
   c=(a+b)/2;
   i=i+1;
end
%求各點(diǎn)橫坐標(biāo)
x(1)=sqrt(98*h/2)*c*h/sqrt(1.0-c*c*v*v);
for k=2:n
    v=(sqrt(2*g*k*h)+sqrt(2*g*(k-1)*h))/2;
    x(k)=x(k-1)+c*v*h/sqrt(1.0-c*c*v*v);
end
%畫模擬計(jì)算出的最速落徑（星點(diǎn)表示）
plot(x,-(0.02:.02:2),'r')
%畫理論計(jì)算出的最速落徑（虛線表示）
hold on
y=0:0.01:2;
x1=acos(1-y)-sqrt(2*y-y.^2);
plot(x1,-y,':b')
axis equal
grid on