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

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

	res=sin(a1/a2)/(2+cos(a2+a3));

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