function [] = lagrangian  ()
	X = -10:0.1:10;
	Y = -10:0.1:10;
for i=X                                
	for j=Y
		M2(round(10*i+101), round(10*j+101)) = i^2*3 + j^2*7;
	end;
end;                                             

surf(X,Y,M2)  