fneval.m

用MATLAB编写关于飞机设计的matlab程序

function [FY]=FNEVAL(FX, P3, F)
%
%	***  Independent variable is FX, returned variable is FY ***
%
%      Real FX, FY, P3, D1, F(0:10,0:10)
%      Integer M
%
     
      M = 1;
      while (FX>(F(1,M+1) + F(6,M) / P3))  %5     If (FX.gt.(F(0,M+1) + F(5,M) / P3)) Then
             M = M + 1;                               %M = M + 1;
                                                  %Goto 5	
      end                                         %Endif(找到FX的取值范围) M是的区间下限。
      D1 = FX - F(1,M);%半径位置的微变量
      FY = F(2,M) + F(3,M)*D1 + F(4,M)*(D1^2) + F(5,M)*(D1^3);
      %Return
%End