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

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

double s,s2,s3,s5,r2,r4,r6,w;
s2=x1*x1+x2*x2; s=sqrt(s2); s3=s*s2; s5=s3*s2;
r2=x1*x1+x2*x2+x3*x3; r4=r2*r2; r6=r4*r2;

w=atan2(x3,s);
cout<<"res   "<<res.get()<<"\t"<<w-res.get()<<endl;
w=-x1*x3/(s*r2);
cout<<"dx1   "<<res.get(0)<<"\t"<<w-res.get(0)<<endl;
w=-x2*x3/(s*r2);
cout<<"dx2   "<<res.get(1)<<"\t"<<w-res.get(1)<<endl;
w=s/r2;
cout<<"dx3   "<<res.get(2)<<"\t"<<w-res.get(2)<<endl;
w=-x3/(s*r2)+x1*x1*x3/(s3*r2)+2*x1*x1*x3/(s*r4) ;
cout<<"dx11  "<<res.get(0,0)<<"\t"<<w-res.get(0,0)<<endl;
w=x1*x2*x3/(s3*r2)+2*x1*x2*x3/(s*r4);
cout<<"dx12  "<<res.get(0,1)<<"\t"<<w-res.get(0,1)<<endl;
w=2*x1*x3*x3/(s*r4)-x1/(s*r2);
cout<<"dx13  "<<res.get(0,2)<<"\t"<<w-res.get(0,2)<<endl;
w=-x3/(s*r2)+x2*x2*x3/(s3*r2)+2*x2*x2*x3/(s*r4);
cout<<"dx22  "<<res.get(1,1)<<"\t"<<w-res.get(1,1)<<endl;
w=2*x2*x3*x3/(s*r4)-x2/(s*r2);
cout<<"dx23  "<<res.get(1,2)<<"\t"<<w-res.get(1,2)<<endl;
w=-2*s*x3/r4;
cout<<"dx33  "<<res.get(2,2)<<"\t"<<w-res.get(2,2)<<endl;
w=3*x1*x3/(s3*r2)-3*x1*x1*x1*x3/(s5*r2)+6*x1*x3/(s*r4)-4*x1*x1*x1*x3/(s3*r4)-8*x1*x1*x1*x3/(s*r6);
cout<<"dx111 "<<res.get(0,0,0)<<"\t"<<w-res.get(0,0,0)<<endl;
w=x2*x3/(s3*r2)-3*x1*x1*x2*x3/(s5*r2)+2*x2*x3/(s*r4)-4*x1*x1*x2*x3/(s3*r4)-8*x1*x1*x2*x3/(s*r6);
cout<<"dx112 "<<res.get(0,0,1)<<"\t"<<w-res.get(0,0,1)<<endl;
w=-1/(s*r2)+x1*x1/(s3*r2)+2*x3*x3/(s*r4)-2*x1*x1*x3*x3/(s3*r4)+2*x1*x1/(s*r4)-8*x1*x1*x3*x3/(s*r6);
cout<<"dx113 "<<res.get(0,0,2)<<"\t"<<w-res.get(0,0,2)<<endl;
w=x1*x3/(s3*r2)-3*x1*x2*x2*x3/(s5*r2)+2*x1*x3/(s*r4)-4*x1*x2*x2*x3/(s3*r4)-8*x1*x2*x2*x3/(s*r6);
cout<<"dx122 "<<res.get(0,1,1)<<"\t"<<w-res.get(0,1,1)<<endl;
w=x1*x2/(s3*r2)-2*x1*x2*x3*x3/(s3*r4)+2*x1*x2/(s*r4)-8*x1*x2*x3*x3/(s*r6);
cout<<"dx123 "<<res.get(0,1,2)<<"\t"<<w-res.get(0,1,2)<<endl;
w=6*x1*x3/(s*r4)-8*x1*x3*x3*x3/(s*r6);
cout<<"dx133 "<<res.get(0,2,2)<<"\t"<<w-res.get(0,2,2)<<endl;
w=3*x2*x3/(s3*r2)-3*x2*x2*x2*x3/(s5*r2)+6*x2*x3/(s*r4)-4*x2*x2*x2*x3/(s3*r4)-8*x2*x2*x2*x3/(s*r6);
cout<<"dx222 "<<res.get(1,1,1)<<"\t"<<w-res.get(1,1,1)<<endl;
w=-1/(s*r2)+x2*x2/(s3*r2)+2*x3*x3/(s*r4)-2*x2*x2*x3*x3/(s3*r4)+2*x2*x2/(s*r4)-8*x2*x2*x3*x3/(s*r6);
cout<<"dx223 "<<res.get(1,1,2)<<"\t"<<w-res.get(1,1,2)<<endl;
w=6*x2*x3/(s*r4)-8*x2*x3*x3*x3/(s*r6);;
cout<<"dx233 "<<res.get(1,2,2)<<"\t"<<w-res.get(1,2,2)<<endl;
w=8*s*x3*x3/r6-2*s/r4;
cout<<"dx333 "<<res.get(2,2,2)<<"\t"<<w-res.get(2,2,2)<<endl;
return 0;
}