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

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

	res=(a1*a4-a2*a3)/(1+a1*a1+a2*a2+a3*a3+a4*a4);

double zn,zn2,zn3,zn4,w,q,d;
zn=(1+x1*x1+x2*x2+x3*x3+x4*x4);
zn2=zn*zn;
zn3=zn2*zn;
zn4=zn2*zn2;
w=(x1*x4-x2*x3)/(1+x1*x1+x2*x2+x3*x3+x4*x4);
q=res.get();
d=w-q;
cout<<"res   "<<q<<"\t"<<d<<"\t"<<endl;
w=x4/zn-2*x1*(x1*x4-x2*x3)/zn2;
q=res.get(0);
d=w-q;
cout<<"dx1   "<<q<<"\t"<<d<<"\t"<<endl;
w=-x3/zn-2*x2*(x1*x4-x2*x3)/zn2;
q=res.get(1);
d=w-q;
cout<<"dx2   "<<q<<"\t"<<d<<"\t"<<endl;
w=-x2/zn-2*x3*(x1*x4-x2*x3)/zn2;
q=res.get(2);
d=w-q;
cout<<"dx3   "<<q<<"\t"<<d<<"\t"<<endl;
w=x1/zn-2*x4*(x1*x4-x2*x3)/zn2;
q=res.get(3);
d=w-q;
cout<<"dx4   "<<q<<"\t"<<d<<"\t"<<endl;
w=-2*(x1*x4-x2*x3)/zn2-4*x1*x4/zn2+8*x1*x1*(x1*x4-x2*x3)/zn3;
q=res.get(0,0);
d=w-q;
cout<<"dx11  "<<q<<"\t"<<d<<"\t"<<endl;
w=-2*x2*x4/zn2+2*x1*x3/zn2+8*x1*x2*(x1*x4-x2*x3)/zn3;
q=res.get(0,1);
d=w-q;
cout<<"dx12  "<<q<<"\t"<<d<<"\t"<<endl;
w=-2*x3*x4/zn2+2*x1*x2/zn2+8*x1*x3*(x1*x4-x2*x3)/zn3;
q=res.get(0,2);
d=w-q;
cout<<"dx13  "<<q<<"\t"<<d<<"\t"<<endl;
w=1/zn-2*x4*x4/zn2-2*x1*x1/zn2+8*x1*x4*(x1*x4-x2*x3)/zn3;
q=res.get(0,3);
d=w-q;
cout<<"dx14  "<<q<<"\t"<<d<<"\t"<<endl;
w=-2*(x1*x4-x2*x3)/zn2+4*x2*x3/zn2+8*x2*x2*(x1*x4-x2*x3)/zn3;
q=res.get(1,1);
d=w-q;
cout<<"dx22  "<<q<<"\t"<<d<<"\t"<<endl;
w=-1/zn+2*x3*x3/zn2+2*x2*x2/zn2+8*x2*x3*(x1*x4-x2*x3)/zn3;
q=res.get(1,2);
d=w-q;
cout<<"dx23  "<<q<<"\t"<<d<<"\t"<<endl;
w=2*x3*x4/zn2-2*x1*x2/zn2+8*x2*x4*(x1*x4-x2*x3)/zn3;
q=res.get(1,3);
d=w-q;
cout<<"dx24  "<<q<<"\t"<<d<<"\t"<<endl;
w=-2*(x1*x4-x2*x3)/zn2+4*x2*x3/zn2+8*x3*x3*(x1*x4-x2*x3)/zn3;
q=res.get(2,2);
d=w-q;
cout<<"dx33  "<<q<<"\t"<<d<<"\t"<<endl;
w=2*x2*x4/zn2-2*x1*x3/zn2+8*x3*x4*(x1*x4-x2*x3)/zn3;
q=res.get(2,3);
d=w-q;
cout<<"dx34  "<<q<<"\t"<<d<<"\t"<<endl;
w=-2*(x1*x4-x2*x3)/zn2-4*x1*x4/zn2+8*x4*x4*(x1*x4-x2*x3)/zn3;
q=res.get(3,3);
d=w-q;
cout<<"dx44  "<<q<<"\t"<<d<<"\t"<<endl;
w=-6*x4/zn2+24*x1*(x1*x4-x2*x3)/zn3
  +24*x1*x1*x4/zn3-48*x1*x1*x1*(x1*x4-x2*x3)/zn4;
q=res.get(0,0,0);
d=w-q;
cout<<"dx111 "<<q<<"\t"<<d<<"\t"<<endl;
w=2*x3/zn2+8*x2*(x1*x4-x2*x3)/zn3
  +16*x1*x2*x4/zn3-8*x1*x1*x3/zn3
  -48*x1*x1*x2*(x1*x4-x2*x3)/zn4;
q=res.get(0,0,1);
d=w-q;
cout<<"dx112 "<<q<<"\t"<<d<<"\t"<<endl;
w=2*x2/zn2+8*x3*(x1*x4-x2*x3)/zn3
 +16*x1*x3*x4/zn3-8*x1*x1*x2/zn3
 -48*x1*x1*x3*(x1*x4-x2*x3)/zn4;
q=res.get(0,0,2);
d=w-q;
cout<<"dx113 "<<q<<"\t"<<d<<"\t"<<endl;
w=-6*x1/zn2+16*x1*x4*x4/zn3
  +8*x4*(x1*x4-x2*x3)/zn3+8*x1*x1*x1/zn3
  -48*x1*x1*x4*(x1*x4-x2*x3)/zn4;
q=res.get(0,0,3);
d=w-q;
cout<<"dx114 "<<q<<"\t"<<d<<"\t"<<endl;
w=-2*x4/zn2+8*x1*(x1*x4-x2*x3)/zn3
  +8*x2*x2*x4/zn3-16*x1*x2*x3/zn3
  -48*x1*x2*x2*(x1*x4-x2*x3)/zn4;
q=res.get(0,1,1);
d=w-q;
cout<<"dx122 "<<q<<"\t"<<d<<"\t"<<endl;
w=2*x1/zn2+8*x2*x3*x4/zn3-8*x1*x3*x3/zn3
 -8*x1*x2*x2/zn3-48*x1*x2*x3*(x1*x4-x2*x3)/zn4;
q=res.get(0,1,2);
d=w-q;
cout<<"dx123 "<<q<<"\t"<<d<<"\t"<<endl;
w=-2*x2/zn2+8*x2*x4*x4/zn3
  -8*x1*x3*x4/zn3+8*x1*x1*x2/zn3
  -48*x1*x2*x4*(x1*x4-x2*x3)/zn4;
q=res.get(0,1,3);
d=w-q;
cout<<"dx124 "<<q<<"\t"<<d<<"\t"<<endl;
w=-2*x4/zn2+8*x1*(x1*x4-x2*x3)/zn3
  +8*x3*x3*x4/zn3-16*x1*x2*x3/zn3
  -48*x1*x3*x3*(x1*x4-x2*x3)/zn4;
q=res.get(0,2,2);
d=w-q;
cout<<"dx133 "<<q<<"\t"<<d<<"\t"<<endl;
w=-2*x3/zn2+8*x3*x4*x4/zn3
  -8*x1*x2*x4/zn3+8*x1*x1*x3/zn3
  -48*x1*x3*x4*(x1*x4-x2*x3)/zn4;
q=res.get(0,2,3);
d=w-q;
cout<<"dx134 "<<q<<"\t"<<d<<"\t"<<endl;
w=-6*x4/zn2+8*x4*x4*x4/zn3
  +8*x1*(x1*x4-x2*x3)/zn3+16*x1*x1*x4/zn3
  -48*x1*x4*x4*(x1*x4-x2*x3)/zn4;
q=res.get(0,3,3);
d=w-q;
cout<<"dx144 "<<q<<"\t"<<d<<"\t"<<endl;
w=6*x3/zn2+24*x2*(x1*x4-x2*x3)/zn3
  -24*x2*x2*x3/zn3
  -48*x2*x2*x2*(x1*x4-x2*x3)/zn4;
q=res.get(1,1,1);
d=w-q;
cout<<"dx222 "<<q<<"\t"<<d<<"\t"<<endl;
w=6*x2/zn2+8*x3*(x1*x4-x2*x3)/zn3
  -16*x2*x3*x3/zn3-8*x2*x2*x2/zn3
  -48*x2*x2*x3*(x1*x4-x2*x3)/zn4;
q=res.get(1,1,2);
d=w-q;
cout<<"dx223 "<<q<<"\t"<<d<<"\t"<<endl;
w=-2*x1/zn2+8*x4*(x1*x4-x2*x3)/zn3
  -16*x2*x3*x4/zn3+8*x1*x2*x2/zn3
  -48*x2*x2*x4*(x1*x4-x2*x3)/zn4;
q=res.get(1,1,3);
d=w-q;
cout<<"dx224 "<<q<<"\t"<<d<<"\t"<<endl;
w=6*x3/zn2+8*x2*(x1*x4-x2*x3)/zn3
  -8*x3*x3*x3/zn3-16*x2*x2*x3/zn3
  -48*x2*x3*x3*(x1*x4-x2*x3)/zn4;
q=res.get(1,2,2);
d=w-q;
cout<<"dx233 "<<q<<"\t"<<d<<"\t"<<endl;
w=2*x4/zn2-8*x3*x3*x4/zn3
  -8*x2*x2*x4/zn3+8*x1*x2*x3/zn3
  -48*x2*x3*x4*(x1*x4-x2*x3)/zn4;
q=res.get(1,2,3);
d=w-q;
cout<<"dx234 "<<q<<"\t"<<d<<"\t"<<endl;
w=2*x3/zn2-8*x3*x4*x4/zn3
  +8*x2*(x1*x4-x2*x3)/zn3+16*x1*x2*x4/zn3
  -48*x2*x4*x4*(x1*x4-x2*x3)/zn4;
q=res.get(1,3,3);
d=w-q;
cout<<"dx244 "<<q<<"\t"<<d<<"\t"<<endl;
w=6*x2/zn2
  +24*x3*(x1*x4-x2*x3)/zn3-24*x2*x3*x3/zn3
  -48*x3*x3*x3*(x1*x4-x2*x3)/zn4;
q=res.get(2,2,2);
d=w-q;
cout<<"dx333 "<<q<<"\t"<<d<<"\t"<<endl;
w=-2*x1/zn2+8*x4*(x1*x4-x2*x3)/zn3
  -16*x2*x3*x4/zn3+8*x1*x3*x3/zn3
  -48*x3*x3*x4*(x1*x4-x2*x3)/zn4;
q=res.get(2,2,3);
d=w-q;
cout<<"dx334 "<<q<<"\t"<<d<<"\t"<<endl;
w=2*x2/zn2-8*x2*x4*x4/zn3
  +8*x3*(x1*x4-x2*x3)/zn3+16*x1*x3*x4/zn3
  -48*x3*x4*x4*(x1*x4-x2*x3)/zn4;
q=res.get(2,3,3);
d=w-q;
cout<<"dx344 "<<q<<"\t"<<d<<"\t"<<endl;
w=-6*x1/zn2+24*x1*x4*x4/zn3
  +24*x4*(x1*x4-x2*x3)/zn3
  -48*x4*x4*x4*(x1*x4-x2*x3)/zn4;
q=res.get(3,3,3);
d=w-q;
cout<<"dx444 "<<q<<"\t"<<d<<"\t"<<endl;

return 0;
}