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

int main()
{
double x[3], x1=2.7, x2=-1.32, x3=-5.1;
// Внимание: массив ext_value3
ext_value3 a[3],res;
int i;

x[0]=x1; x[1]=x2; x[2]=x3;
// Инициализация массива:
for(i=0;i<3;i++){
 a[i].set_dim(3); // установка размерности
 a[i].set(x[i]);  // установка значения
 a[i].set(1.0,i); 
 // устанавливаем первую производную - получили независимые переменные
 // по этим переменным и производится дифференцирование
}
   res=atan2(a[2],sqrt(a[0]*a[0]+a[1]*a[1]));

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;
}