c1 = -100;
c2 = -10;
c3 = 42;
a = @(x) x^2 + c1;
z = @(x, y) a(x) + b(x) + c(y) + d(y) + c3;
[X, Y] = meshgrid([-10:0.1:10], [-10:0.1:10]);
surf(X, Y, z(X, Y));
YzEgPSAtMTAwOwpjMiA9IC0xMDsKYzMgPSA0MjsKYSA9IEAoeCkgeF4yICsgYzE7CmIgPSBAY29zOwpjID0gQCh4KSBzcXJ0KHggLSBjMik7CmQgPSBAZXhwOwp6ID0gQCh4LCB5KSBhKHgpICsgYih4KSArIGMoeSkgKyBkKHkpICsgYzM7CltYLCBZXSA9IG1lc2hncmlkKFstMTA6MC4xOjEwXSwgWy0xMDowLjE6MTBdKTsKc3VyZihYLCBZLCB6KFgsIFkpKTs=