Next: Implementation with gradient
Up: Implementation without gradient
Previous: Using an expansion
Contents
The main program for the previous system of 2 equations
(omitting the procedure F that computes the interval evaluation
of the equations) will be
written as
int main()
{
int num,i;
INTERVAL_MATRIX Solution(100,2);
INTERVAL_VECTOR x(2);
INTEGER_VECTOR Type(2);
Clear(Type);
for(i=1;i<=2;i++)x(i)=INTERVAL(-10,10);
Single_Bisection=1;
num=Solve_Simplex(2,2,0,Type,F,NonLinear,CoeffLinear,0.001,0,x,
0,40000,0,0.,0.0001,0.0001,Solution,100);
for(i=1;i<=num;i++)
{
x(1)=Solution(i,1);x(2)=Solution(i,2);
cout<<"Solution "<<i<<":"<<x<<endl;
}
}
Jean-Pierre Merlet
2012-12-20