13.5.2  约束条件

有时，约束条件不是自变量，而是由积分变量导出的某个参数，数学形式是：

H = F( X ) 

定步长积分生成的一系列数据点可能会跨过终止点，要满足这类终止条件往往需要在终点处进行迭代。例如，积分终点示意图如图13-2所示。

假设计算一个弹头的弹道积分，目标的高程是H0，那么，积分终止条件显然是：

H = H0 

H = H( X ) 

while( 1 )  
{  
    X(t) => Y( t+h )     // 单步积分，Y是临时变量  
    if( h < EPS )        // 步长很小，到达终点，退出循环  
    {  
        X( t + h ) = Y( t + h );  
        t += h;  
        break;  
    }  
 
    if( H(X) < H0 )   // 弹头越过目标，步长需要减半  
        h /= 2.0;  
    else{  
        X( t+h ) = Y( t+h );// 尚未到终点，正常继续积分  
        t = t + h;  
    }  
}