# Exercici 11.7

kill(all)$

## Apartat 1)

f(x,y,z):= x*y*z+sin(z-2)-x^2*y^2-x-y+1;

'f(1,1,2)=f(1,1,2);

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

'Jf(1,1,2)=Jf112:subst([x=1,y=1,z=2],Jfxyz);

load ("eigen")$     H=H:columnvector([x-1,y-1,z-2]);

Jf112.H = 0;

Per tant l'aproximació lineal és x+y-2z+2=0

## Apartat 2)

Pla tangent en el punt P=(1,1,2):

Jf112.H = 0;

Recta normal a la superficie en el punt P=(1,1,2)  (equacions paramètriques) :

x = 1 + %lambda;
y = 1 + %lambda;
z = 2 - 2*%lambda;

## Apartat 3)

Comprovarem que es compleixin les condicions de la funció implícita

1.La funció f és del tipus Cr amb r>=1

2. El punt compleix l'equacio:

'f(1,1,2)=f(1,1,2);

3. El determinant jacobia de la matriu de derivades respecte la variable que es vol demostrar que es defineix com funció implícita no és nul:

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

'Jf(1,1,2)=Jf112:subst([x=1,y=1,z=2],Jfz);

Per tant es compleixen les condicions del teorema de la funció implícita

depends([z],[x,y]);

DxEq1:diff(x*y*z+sin(z-2)-x^2*y^2-x-y+1,x);
DyEq1:diff(x*y*z+sin(z-2)-x^2*y^2-x-y+1,y);

Considerant el punt (x,y)=(1,1)on z(1,1)=2 ens quedarà:

'DxEq11=DxEq11:cos(2-2)*('diff(z,x,1))+('diff(z,x,1))+2-2-1;
'DyEq12=DyEq12:cos(2-2)*('diff(z,y,1))+('diff(z,y,1))+2-2-1;

Sist01:[DxEq11,DyEq12];

solve(Sist01,['diff(z,x),'diff(z,y)]);

'Jz(1,1)=Jzxy:matrix([1/2,1/2]);

load ("eigen")$
H1=H1:columnvector([x-1,y-1]);

Per tant el polinomi de Taylor de z en el punt (1,1) serà:

T1=2+Jzxy.H1;

ratsimp(%);