#include <cstdlib>
#include <cmath>
#include <iostream>
#include <stdio.h>
#include "ext_value.h"
using namespace std;

void fcn( int n, ext_value1 *ymas, ext_value1 *f);
int dopri8(int n,double x,ext_value1 *y,double xend,double eps,double hmax,double h);
void gauss(double *a, double *b, int n);
int resh(int n1,int n2,ext_value1 p[]);
void f(int n1,int n2,ext_value1 p[],ext_value1 q[]);

#include "gauss.cpp"
#include "resh.cpp"
#include "fcn.cpp"
#include "d87-ext1.cpp"
#include "f.cpp"

int main()
{  
int i,n1=3,n2=2,iter;
double alpha=25.0,beta0=0.5,beta1=0.5; 
ext_value1 p[n1];

printf("alpha=%lf:\n",alpha);

for(i=0;i<n1;i++)  p[i].set_dim(n2); // установка размерности
 p[0].set(beta0); p[0].set(1.0,0); 
 p[1].set(beta1); p[1].set(1.0,1); 
 p[2].set(alpha);

// printf("resh start\n");
iter=resh(n1,n2,p);
// printf("resh end\n");

if(iter<0){
printf("задачу решить не удалось\n");
return 1;
}
// здесь вывод результатов
    printf("Получен ответ iter=%d:\n",iter);
	printf("beta0=%lf\n",p[0].get() );
	printf("beta1=%lf\n",p[1].get() );
return 0;
}