[x,y] = meshgrid(-1:.1:1, -1:.1:1);
z=(x.^2+y.^2).*cos(x.^2+y.^2);
surf(x,y,z)
Free and Online
https://octave-online.net/
