// russian : cp1251 encoding
#include <iostream>
#include <cmath>
#include "ext_value.h"
using namespace std;

int main()
{
double x1=1.1, x2=2.3, x3=-2, w;
// Внимание! Самая важная строчка!
// Инициализация независимых переменных
// Всего 3 независимых переменных (значение первого параметра): 
// значит индексы 0,1,2
ext_value3 a1(3,x1,0),a2(3,x2,1),a3(3,x3,2);
ext_value3 res(3); // Здесь будет вычисленный результат
// Реализуется тот же самый пример, но с использованием других операторов.
ext_value3 arg1(3),arg2(3);

arg1=a1;
arg1+=a2;
arg1=sin(arg1);
arg1*=a1;
arg1+=1;
arg1=arg1*2;
arg2=a2;
arg2*=a1;
arg2=cos(arg2);
arg2*=a3;
arg2-=1;
arg2=2*arg2;
res=arg1+arg2;
res*=0.5;

w=cos(x1*x2)*x3 + x1*sin(x2+x1);
cout<<"res   "<<res.get()<<"\t"<<w-res.get()<<endl;
w=-x2*sin(x1*x2)*x3 + sin(x2+x1) + x1*cos(x2+x1);
cout<<"dx1   "<<res.get(0)<<"\t"<<w-res.get(0)<<endl;
w=x1*cos(x2+x1) - x1*sin(x1*x2)*x3;
cout<<"dx2   "<<res.get(1)<<"\t"<<w-res.get(1)<<endl;
w=cos(x1*x2);
cout<<"dx3   "<<res.get(2)<<"\t"<<w-res.get(2)<<endl;
w=-x2*x2*cos(x1*x2)*x3 - x1*sin(x2+x1) + 2*cos(x2+x1);
cout<<"dx11  "<<res.get(0,0)<<"\t"<<w-res.get(0,0)<<endl;
w=-sin(x1*x2)*x3 - x1*x2*cos(x1*x2)*x3 - x1*sin(x2+x1) + cos(x2+x1);
cout<<"dx12  "<<res.get(0,1)<<"\t"<<w-res.get(0,1)<<endl;
w=-x2*sin(x1*x2);
cout<<"dx13  "<<res.get(0,2)<<"\t"<<w-res.get(0,2)<<endl;
w=-x1*x1*cos(x1*x2)*x3 - x1*sin(x2+x1);
cout<<"dx22  "<<res.get(1,1)<<"\t"<<w-res.get(1,1)<<endl;
w=-x1*sin(x1*x2);
cout<<"dx23  "<<res.get(1,2)<<"\t"<<w-res.get(1,2)<<endl;
w=0;
cout<<"dx33  "<<res.get(2,2)<<"\t"<<w-res.get(2,2)<<endl;
w=x2*x2*x2*sin(x1*x2)*x3 - 3*sin(x2+x1) - x1*cos(x2+x1);
cout<<"dx111 "<<res.get(0,0,0)<<"\t"<<w-res.get(0,0,0)<<endl;
w=x1*x2*x2*sin(x1*x2)*x3 - 2*x2*cos(x1*x2)*x3 - 2*sin(x2+x1) - x1*cos(x2+x1);
cout<<"dx112 "<<res.get(0,0,1)<<"\t"<<w-res.get(0,0,1)<<endl;
w=-x2*x2*cos(x1*x2);
cout<<"dx113 "<<res.get(0,0,2)<<"\t"<<w-res.get(0,0,2)<<endl;
w=x1*x1*x2*sin(x1*x2)*x3 - 2*x1*cos(x1*x2)*x3 - sin(x2+x1) - x1*cos(x2+x1);
cout<<"dx122 "<<res.get(0,1,1)<<"\t"<<w-res.get(0,1,1)<<endl;
w=-sin(x1*x2) - x1*x2*cos(x1*x2);
cout<<"dx123 "<<res.get(0,1,2)<<"\t"<<w-res.get(0,1,2)<<endl;
w=0;
cout<<"dx133 "<<res.get(0,2,2)<<"\t"<<w-res.get(0,2,2)<<endl;
w=x1*x1*x1*sin(x1*x2)*x3 - x1*cos(x2 + x1);
cout<<"dx222 "<<res.get(1,1,1)<<"\t"<<w-res.get(1,1,1)<<endl;
w=-x1*x1*cos(x1*x2);
cout<<"dx223 "<<res.get(1,1,2)<<"\t"<<w-res.get(1,1,2)<<endl;
w=0;
cout<<"dx233 "<<res.get(1,2,2)<<"\t"<<w-res.get(1,2,2)<<endl;
w=0;
cout<<"dx333 "<<res.get(2,2,2)<<"\t"<<w-res.get(2,2,2)<<endl;
return 0;
}