Метод ломаных нахождения минимума многоэкстремальной функции

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <conio.h>
#pragma hdrstop
#pragma argsused
//---------------------------------------------------------------------------
const int N=26;
const double M=180;
double func(double x)
{
double y;
y=pow(x,4)-8*pow(x,3)-8*pow(x,2)+96*x-3;
return y;
}
int minpos(double a[N+1], int n)
{
  int imin;
  double min;
  imin=1;
  min=a[1];
  for (int i=2; i<=n; i++)
  {
     if (a[i]<min)
     {
       min=a[i];
       imin=i;
     }
  }
  return imin;
}
double min(double a[N+1],int n)
{
  double min;
  min=a[1];
  for (int i=2; i<=n; i++)
  {
     if (a[i]<min)
     {
       min=a[i];
     }
  }
  return min;
}
double InterPoint(double InterPoint[3],double x1, double y1, double x2, double y2)
{
 double Temp;
 Temp = (y1 - y2)/(2*M) + (x1 + x2)/2;
 InterPoint[1] = Temp;
 InterPoint[2] = y1 + M*x1 - M*Temp;
 return *InterPoint;
}
 
int main(int argc, char* argv[])
{
 double a=-3, b=7;
 double x[N+1];
 double y[N+1];
 double ymin;
 int i, imin, j;
 printf ("Metod perebora:\n");
 for (i=1;i<=N;i++)
 {
  x[i]=a+(2*i-1)*(b-a)/(2*N);
 }
 for (i=1; i<=N; i++)
 y[i]=func(x[i]);
 ymin=y[1];
 for (i=1; i<=N; i++)
 {
    if (y[i]<=ymin)
     {
         ymin=y[i];
     imin=i;
     }
 
 }
 printf("Naim. znacheniye f: %f\n\n", ymin);
 printf ("Metod lomanyh:\n");
 double Theor, Fact;
 double Inter[3];
 double minPhix[N], minPhiy[N];
 int minNum;
 x[1] = a;
 x[2] = b;
 y[1] = func(x[1]);
 y[2] = func(x[2]);
 InterPoint(Inter,x[1], y[1], x[2], y[2]);
 minPhix[1] = Inter[1];
 minPhiy[1] = Inter[2];
 for( i = 3; i<=N;i++)
 {
   minNum = minpos(minPhiy,i-2);
   for(j = 1;j<=(i-minNum-1);j++)
   {
     x[i-j+1] = x[i-j];
     y[i-j+1] = y[i-j];
   }
   x[minNum+1] = minPhix[minNum];
   y[minNum+1] = func(x[minNum+1]);
   for (j = 1; j<=(i-minNum-1); j++)
   {
     minPhix[i-j] = minPhix[i-j-1];
     minPhiy[i-j] = minPhiy[i-j-1];
   }
   InterPoint(Inter,x[minNum], y[minNum], x[minNum+1], y[minNum+1]);
   minPhix[minNum] = Inter[1];
   minPhiy[minNum] = Inter[2];
   InterPoint(Inter,x[minNum+1], y[minNum+1], x[minNum+2], y[minNum+2]);
   minPhix[minNum+1] = Inter[1];
   minPhiy[minNum+1] = Inter[2];
 }
 ymin = min(y,N);
 Fact = fabs(ymin-(-147));
 Theor = ymin - min(minPhiy,N);
 printf("Naim. znacheniye f: %f\n", ymin);
 printf("Theor. precision is: %f\n", Theor);
 printf("Fact. precision is: %f\n", Fact);
 getch();
 return 0;
}