C程序高手进来 帮帮忙 谢谢

#include<math.h>
#include<stdio.h>
#include<stdlib.h>
#define E1 0.0001
#define ep 0.0001
#define n 2
#define k 4
float af=10;
int i,j;
double X0[n],XX[n],X[k][n],FF[k];
void produce(double A[n],double B[n]);
double a[n],b[n];
double F(double C[n]);
int cons (double D[n]);
void bou();

double F(double C[n])
{
double F;
F=25/(C[0]*C[1]*C[1]*C[1]);
return F;
}
int cons(double D[n])
{
if(((30/(D[0]*D[1]*D[1])-50)<=0)&&((0.0004*D[0]*D[1]-0.001)<=0)&&((D[0]>=2)&&(D[0]<=4))&&((D[1]>=0.5)&&(D[1]<=1)))
return 1;
else return 0;
}
void bou()
{ a[0]=2;a[1]=0.5;
b[0]=4;b[1]=1;
}
void produce(double A[],double B[])
{
double rr;
int jj;
do
{
rr=rand();
rr=rr/32767;
for(i=0;i<n;i++)
{