# Hello World program in Python
import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
 
fig = plt.figure()
ax = Axes3D(fig)
 
# X, Y value
X = np.arange(-4, 4, 0.25)
Y = np.arange(-4, 4, 0.25)
X, Y = np.meshgrid(X, Y)    
R = np.sqrt(X ** 2 + Y ** 2)
# height value
Z = np.sin(R)
 
ax.plot_surface(X, Y, Z, rstride=1, cstride=1, cmap=plt.get_cmap('rainbow'))
 
ax.contourf(X, Y, Z, zdir='z', offset=-2, cmap=plt.get_cmap('rainbow'))
 
 

IyBIZWxsbyBXb3JsZCBwcm9ncmFtIGluIFB5dGhvbgppbXBvcnQgbnVtcHkgYXMgbnAKaW1wb3J0IG1hdHBsb3RsaWIucHlwbG90IGFzIHBsdApmcm9tIG1wbF90b29sa2l0cy5tcGxvdDNkIGltcG9ydCBBeGVzM0QKCmZpZyA9IHBsdC5maWd1cmUoKQpheCA9IEF4ZXMzRChmaWcpCgojIFgsIFkgdmFsdWUKWCA9IG5wLmFyYW5nZSgtNCwgNCwgMC4yNSkKWSA9IG5wLmFyYW5nZSgtNCwgNCwgMC4yNSkKWCwgWSA9IG5wLm1lc2hncmlkKFgsIFkpICAgIApSID0gbnAuc3FydChYICoqIDIgKyBZICoqIDIpCiMgaGVpZ2h0IHZhbHVlClogPSBucC5zaW4oUikKCmF4LnBsb3Rfc3VyZmFjZShYLCBZLCBaLCByc3RyaWRlPTEsIGNzdHJpZGU9MSwgY21hcD1wbHQuZ2V0X2NtYXAoJ3JhaW5ib3cnKSkKCmF4LmNvbnRvdXJmKFgsIFksIFosIHpkaXI9J3onLCBvZmZzZXQ9LTIsIGNtYXA9cGx0LmdldF9jbWFwKCdyYWluYm93JykpCgo=