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

int main()
{
double x1=2.1, x2=-1.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; // Здесь будет вычисленный результат
    res=exp(-sqrt(a1*a1+a2*a2+a3*a3));

double s,s2,s3,s4,s5,e,w;
s2=x1*x1+x2*x2+x3*x3;
s=sqrt(s2);
s3=s*s2;
s4=s2*s2;
s5=s3*s2;
e=exp(-s);
w=e;
cout<<"res   "<<res.get()<<"\t"<<w-res.get()<<endl;
w=-x1*e/s;
cout<<"dx1   "<<res.get(0)<<"\t"<<w-res.get(0)<<endl;
w=-x2*e/s;
cout<<"dx2   "<<res.get(1)<<"\t"<<w-res.get(1)<<endl;
w=-x3*e/s;
cout<<"dx3   "<<res.get(2)<<"\t"<<w-res.get(2)<<endl;
w=-e/s+x1*x1*e/s2+x1*x1*e/s3;
cout<<"dx11  "<<res.get(0,0)<<"\t"<<w-res.get(0,0)<<endl;
w=x1*x2*e/s2+x1*x2*e/s3;
cout<<"dx12  "<<res.get(0,1)<<"\t"<<w-res.get(0,1)<<endl;
w=x1*x3*e/s2+x1*x3*e/s3;
cout<<"dx13  "<<res.get(0,2)<<"\t"<<w-res.get(0,2)<<endl;
w=-e/s+x2*x2*e/s2+x2*x2*e/s3;
cout<<"dx22  "<<res.get(1,1)<<"\t"<<w-res.get(1,1)<<endl;
w=x2*x3*e/s2+x2*x3*e/s3;
cout<<"dx23  "<<res.get(1,2)<<"\t"<<w-res.get(1,2)<<endl;
w=-e/s+x3*x3*e/s2+x3*x3*e/s3;
cout<<"dx33  "<<res.get(2,2)<<"\t"<<w-res.get(2,2)<<endl;
w=3*x1*e/s2-x1*x1*x1*e/s3+3*x1*e/s3-3*x1*x1*x1*e/s4-3*x1*x1*x1*e/s5;
cout<<"dx111 "<<res.get(0,0,0)<<"\t"<<w-res.get(0,0,0)<<endl;
w=x2*e/s2-x1*x1*x2*e/s3+x2*e/s3-3*x1*x1*x2*e/s4-3*x1*x1*x2*e/s5;
cout<<"dx112 "<<res.get(0,0,1)<<"\t"<<w-res.get(0,0,1)<<endl;
w=x3*e/s2-x1*x1*x3*e/s3+x3*e/s3-3*x1*x1*x3*e/s4-3*x1*x1*x3*e/s5;
cout<<"dx113 "<<res.get(0,0,2)<<"\t"<<w-res.get(0,0,2)<<endl;
w=x1*e/s2-x2*x2*x1*e/s3+x1*e/s3-3*x2*x2*x1*e/s4-3*x2*x2*x1*e/s5;
cout<<"dx122 "<<res.get(0,1,1)<<"\t"<<w-res.get(0,1,1)<<endl;
w=-x1*x2*x3*e/s3-3*x1*x2*x3*e/s4-3*x1*x2*x3*e/s5;
cout<<"dx123 "<<res.get(0,1,2)<<"\t"<<w-res.get(0,1,2)<<endl;
w=x1*e/s2-x3*x3*x1*e/s3+x1*e/s3-3*x3*x3*x1*e/s4-3*x3*x3*x1*e/s5;
cout<<"dx133 "<<res.get(0,2,2)<<"\t"<<w-res.get(0,2,2)<<endl;
w=3*x2*e/s2-x2*x2*x2*e/s3+3*x2*e/s3-3*x2*x2*x2*e/s4-3*x2*x2*x2*e/s5;
cout<<"dx222 "<<res.get(1,1,1)<<"\t"<<w-res.get(1,1,1)<<endl;
w=x3*e/s2-x2*x2*x3*e/s3+x3*e/s3-3*x2*x2*x3*e/s4-3*x2*x2*x3*e/s5;
cout<<"dx223 "<<res.get(1,1,2)<<"\t"<<w-res.get(1,1,2)<<endl;
w=x2*e/s2-x3*x3*x2*e/s3+x2*e/s3-3*x3*x3*x2*e/s4-3*x3*x3*x2*e/s5;
cout<<"dx233 "<<res.get(1,2,2)<<"\t"<<w-res.get(1,2,2)<<endl;
w=3*x3*e/s2-x3*x3*x3*e/s3+3*x3*e/s3-3*x3*x3*x3*e/s4-3*x3*x3*x3*e/s5;
cout<<"dx333 "<<res.get(2,2,2)<<"\t"<<w-res.get(2,2,2)<<endl;
return 0;
}