Отображение поверхности с помощью surf и линии уровня ezcontour

syms x y;
a=4;
b=5;
 
z= x.*y.*exp((-x.^2)./a - (y.^2)/b);
subplot(2,2,1), contour(x,y,z, 20);
>> subplot(2,2,2), ezcontour('x.*y.*exp((-x.^2)./a - (y.^2)/b)',[-5, 5], 30);
subplot(2,2,3), ezcontour('x.*y.*exp((-x.^2)./a - (y.^2)/b)',[-5, 5], 40);
subplot(2,2,4), ezcontour('x.*y.*exp((-x.^2)./a - (y.^2)/b)',[-5, 5], 50);
 
 
z= x.*y.*exp((-x.^2)./a - (y.^2)/b);
subplot(2,2,1), surf( x, y, z );
subplot(2,2,2), ezsurf('x.*y.*exp((-x.^2)./a - (y.^2)/b)',[-5, 5] );
subplot(2,2,3), ezsurf( 'x.*y.*exp((-x.^2)./a - (y.^2)/b)',[-7,7] );
subplot(2,2,4), ezsurf( 'x.*y.*exp((-x.^2)./a - (y.^2)/b)',[-9, 9]);
 
 
Error using ==> mupadmex
Error in MuPAD command: DOUBLE cannot convert the input expression into a double array.
 
If the input expression contains a symbolic variable, use the VPA function instead.
 
Error in ==> sym.sym>sym.double at 927
            Xstr = mupadmex('mllib::double', S.s, 0);
 
Error in ==> datachk at 11
    y = full(double(x));
 
Error in ==> contour>parseargs at 161
            x = datachk(args{1});
 
Error in ==> contour at 72
            [pvpairs,args,msg] = parseargs(args);

clc, clear all, close all
 
[x,y] = meshgrid(linspace(0,5,50));
 
a = 4;
b = 5;
 
z = x.*y.*exp((-x.^2)./a - (y.^2)/b);
 
fun = @(x,y)(x.*y.*exp((-x.^2)./a - (y.^2)/b));
 
subplot(2,2,1), surf(x,y,z);
subplot(2,2,2), ezsurf(fun,[-5, 5] );
subplot(2,2,3), ezsurf(fun,[-7, 7] );
subplot(2,2,4), ezsurf(fun,[-9, 9]);