powell 优化算法的过程

Powell优化算法是利用仪器测井理建立误差函数（非相关函数），借助Powell方向加速法求出非相关函数达到最小时的解，对于气，水两相流动，从预设的气，水流量初始值出发，沿不同的广向进行搜索，可求出气，水两相流动中可能最大产量。与目前常用的生产测井解释方法相比，文中提出的方法具有精度高，实用性强等优点，在测井曲线有缺陷时，仍有可能得到较好的结果。
powell.c代码如下：
CODE:

＃i nclude "hjfgf.c"
double oneoptim(double x0[],double s[],double h0,double epsg,int n,double x[])
{double *a,*b,ff;
a=(double *)malloc(n*sizeof(double));
b=(double *)malloc(n*sizeof(double));
jtf(x0,h0,s,n,a,b);
ff=gold(a,b,epsg,n,x);
free(a);
free(b);
return (ff);
}
double powell(double p[],double h0,double eps,double epsg,int n,double x[])
{int i,j,m;
double *xx[4],*ss,*s;
double f,f0,f1,f2,f3,fx,dlt,df,sdx,q,d;
ss=(double *)malloc(n*(n+1)*sizeof(double));
s=(double *)malloc(n*sizeof(double));
for(i=0;i<n;i++)
{for(j=0;j<=n;j++)
*(ss+i*(n+1)+j)=0;
*(ss+i*(n+1)+i)=1;
}
for(i=0;i<4;i++)
xx[i]=(double *)malloc(n*sizeof(doubl