Can you convert this code to the C$++$ language?

import matplotlib.pyplot as plt

from mpl$\_$toolkits.mplot$3$d$.$art$3$d import Poly$3$DCollection

import numpy as np

# Create figure and $3$D axis

fig $=$ plt$.$figure$()$

ax $=$ fig.add$\_$subplot$(111,$ projection$=\text{'}3$d$\text{'})$

# Function to draw a cube

def draw$\_$cube$($ax$,$ center, size, color$)$:

r $=$ size $/2\mathrm{.}0$

points $=$ np$.$array$([[$center$[0]-$ r$,$ center$[1]-$ r$,$ center$[2]-$ r$],$

$[$center$[0]+$ r$,$ center$[1]-$ r$,$ center$[2]-$ r$],$

$[$center$[0]+$ r$,$ center$[1]+$ r$,$ center$[2]-$ r$],$

$[$center$[0]-$ r$,$ center$[1]+$ r$,$ center$[2]-$ r$],$

$[$center$[0]-$ r$,$ center$[1]-$ r$,$ center$[2]+$ r$],$

$[$center$[0]+$ r$,$ center$[1]-$ r$,$ center$[2]+$ r$],$

$[$center$[0]+$ r$,$ center$[1]+$ r$,$ center$[2]+$ r$],$

$[$center$[0]-$ r$,$ center$[1]+$ r$,$ center$[2]+$ r$]])$

edges $=[[$points$[$j$]$ for j in $[0,1,2,3]],$

$[$points$[$j$]$ for j in $[4,5,6,7]],$

$[$points$[$j$]$ for j in $[0,1,5,4]],$

$[$points$[$j$]$ for j in $[2,3,7,6]],$

$[$points$[$j$]$ for j in $[1,2,6,5]],$

$[$points$[$j$]$ for j in $[4,7,3,0]]]$

ax$.$add$\_$collection$3$d$($Poly$3$DCollection$($edges$,$ facecolors$=$color, linewidths$=1,$ edgecolors$=\text{'}$r$\text{'},$ alpha$=\mathrm{.}25))$

# Function to draw a cylinder

def draw$\_$cylinder$($ax$,$ center, radius, height, color$)$:

z $=$ np$.$linspace$($center$[2]-$ height $/2,$ center$[2]+$ height $/2,50)$

theta $=$ np$.$linspace$(0,2*$ np$.$pi$,50)$

theta$\_$grid, z$\_$grid $=$ np$.$meshgrid$($theta$,$ z$)$

x$\_$grid $=$ radius $*$ np$.$cos$($theta$\_$grid$)+$ center$[0]$

y$\_$grid $=$ radius $*$np$.$sin$($theta$\_$grid$)+$ center$[1]$

ax$.$plot$\_$surface$($x$\_$grid, y$\_$grid, z$\_$grid, color$=$color$)$

# Function to draw a sphere

def draw$\_$sphere$($ax$,$ center, radius, color$)$:

u$,$ v $=$ np$.$mgrid$[0$:$2*$np$.$pi:$20$j$,0$:np$.$pi:$10$j$]$

x $=$ radius $*$ np$.$cos$($u$)*$ np$.$sin$($v$)+$ center$[0]$

y $=$ radius $*$ np$.$sin$($u$)*$ np$.$sin$($v$)+$ center$[1]$

z $=$ radius $*$ np$.$cos$($v$)+$ center$[2]$

ax$.$plot$\_$surface$($x$,$ y$,$ z$,$ color$=$color$)$

# Drawing shapes based on the provided $2$D image

draw$\_$cube$($ax$,$ center$=[1,1,0\mathrm{.}5],$ size$=1,$ color$=$'white'$)$

draw$\_$cube$($ax$,$ center$=[2\mathrm{.}5,1,1],$ size$=2,$ color$=$'white'$)$

draw$\_$cube$($ax$,$ center$=[-1,1,0\mathrm{.}5],$ size$=0\mathrm{.}5,$ color$=$'white'$)$

draw$\_$cylinder$($ax$,$ center$=[1,1,1\mathrm{.}5],$ radius$=0\mathrm{.}5,$ height$=1,$ color$=$'sandybrown'$)$

draw$\_$sphere$($ax$,$ center$=[1\mathrm{.}5,0,0\mathrm{.}5],$ radius$=0\mathrm{.}5,$ color$=$'pink'$)$

# Set the limits of the axes

ax$.$set$\_$xlim$([-2,4])$

ax$.$set$\_$ylim$([-2,4])$

ax$.$set$\_$zlim$([0,2])$

# Hide the axes

ax$.$axis$(\text{'}$off$\text{'})$

# Show plot

plt$.$show$()$