import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
import numpy as np
import matplotlib.pyplot as plt

fig = plt.figure()
ax = fig.gca(projection='3d')
X = np.arange(1, 8, np.pi)
Y = [-8,-7,-6,-5,-4,-3,-2,-1,1,2,3,4,5,6,7,8]
X, Y = np.meshgrid(X, Y)
Z = ((np.cos(Y)+np.sin(Y))/(Y))
surf = ax.plot_surface(X, Y, Z,linewidth=0, antialiased=False)

fig.suptitle('Wykres funkcji f(x,y)=(np.cos(X)+3)/((X**2)+4)', fontsize=10, fontweight='bold')
ax.set_xlabel('X')
ax.set_ylabel('Y')
ax.set_zlabel('Z')
legend = ax.legend(loc='upper left', shadow=True, fontsize='large')
plt.show()