metodos numericos

Problema 1: graficas de la función dada y maximizar la función
Programa en scilab
deff("z=f(x,y)","z=exp(-(2*y+x+2)^2-(3*x-y)^2)+2*exp(-(2*y+x+1)^2-3*(-y+3*x3)^2)+1/((2*y+x)^2+(3*x-y)^2+1)");fplot3d(-3:0.1:3,-3:0.2:3,f)
scf();//
contour(-3:0.1:3,-3:0.1:3,f,5)

Ilustración 1 curvas de nivel

Ilustración 2 funcion en 3D

Maximizar la función con scilab.
function [f, g,ind]=costf(x, ind)
f=-(exp(-(2*x(2)+x(1)+2)^2-(3*x(1)-x(2))^2)+2*exp(-(2*x(2)+x(1)+1)^2-3*(-x(2)+3*x(1)3)^2)+1/((2*x(2)+x(1))^2+(3*x(1)-x(2))^2+1));g=[(6*(3*x(1)-x(2)+2*(x(1)+2*x(2))))/((3*x(1)-x(2))^2+(x(1)+2*x(2))^2+1)^2-2*exp(3*(3*x(1)-x(2)-3)^2-(x(1)+2+1)^2)*(-18*(3*x(1)-x(2)-3)-2*(x(1)+2*x(2)+1))-exp(-(3*x(1)x(2))^2-(x(1)+2*x(2)+2)^2)*(-6*(3*x(1)-x(2))-2*(x(1)+2*x(2)+2));(4*(x(1)+2*x(2))-2*(3*x(1)x(2)))/((3*x(1-x(2)))^2+(x(1)+2*x(2))^2+1)^2-2*exp(-3*(3*x(1)-x(2)-3)^2(x(1)+2*x(2)+1)^2)*(6*(3*x(1-x(2)-3))-4*(x(1)+2*x(2)+1))-exp(-(3*x(1)-x(2))^2(x(1)+2*x(2)+2)^2)*(2*(3*x(1)-x(2))-4*(x(1)+2*x(2)+2))];
endfunction
x0=[2;1];
[f,xopt]=optim(costf,x0)
--> [f,xopt]=optim(cost,x0)
xopt =
-0.00608377
-0.0182513
f =
1.01987

Problema 2
Encontrar las temperaturas de la placa en ciertospuntos de ella

Metodología y herramienta:
La herramienta utilizada fue el método de interpolación un tema visto en clase y que sirve
para encontrar valores en este caso de temperatura enciertos punto de la placa ya
conocidos.
Primero en el primer inciso fue encontrar el valor de z o temperatura en el punto x=6,
y=4.3 usando interpolación lineal conociendo los valores de latemperatura en x=6
Para el segundo inciso fue lo mismo usando interpolación lineal pero fue un poco más
complejo, primer sacar los valores de la temperatura en x=0 y y=5.7 y así con x=2 etc.…
ya teniendolos valores de la temperatura en el punto de y=5.7 entonces a recurrir a
interpolación lineal para x=6.4
Programa en scilab
Valoresy=[0,2,4,6,8,10];
Temperatura=[70,47.15,33.18,29,26,18];...