# Ejercicio 11.2

kill(all)$

## Apartado a)

f1(x,y):= exp(2*x)-exp(y);   
f2(x,y):= exp(y);
f(x,y):=  [f1(x,y) , f2(x,y)];

Las funciones componentes son de tipo C1 en su campo de existencia  y, en consecuencia, lo será la función f(x,y)=(f1,f2)

'Jf(x,y) = Jfxy: jacobian( f(x,y) , [x,y]) ;

'detJf(x,y) = detJfxy: determinant(Jfxy);

El determinante es > 0 para todo valor de (x,y) de R2.
Por lo tanto, f es localmente invertible y la inversa es diferenciable.

## Apartado b)

Eq1: u = f1(x,y);   
Eq2: v = f2(x,y);
solve([Eq1,Eq2], [x,y]);

Debemos transformar las ecuaciones para poder resolver el sistema

Eq1a: log(u+v) = 2*x;   
Eq2a: log(v) = y;
solve([Eq1a,Eq2a], [x,y]);

Por lo tanto la f^{-1}(u,v)=((1/2)*log(u+v), log(v))

## Apartado c)

finv1(u,v):= (1/2)*log(u+v) ;
finv2(u,v):= log(v) ;
finv(u,v):= [finv1(u,v) , finv2(u,v)];

'Jf_inv(u,v)=Jfinv :jacobian ( finv(u,v) , [u,v]) ;

'Jf_inv(x,y)=Jfinv2: subst ( [u=exp(2*x)-exp(y), v=exp(y)] , Jfinv) ;

Jfxy.Jfinv2 ;

Por lo tanto, queda comprobado que las matrices jacobianas en (x,y) son inversas.