Use Python $($specifically sympy$)$ to solve each problem:

$1.$ Given f$($x$)=$x and g$($x$)=($x $3)^2,$

$($a$)$ Plot both functions on the same axes, with x$-$interval $[0,5].$

$($b$)$ Find the volume of the solid generated by rotating the region bounded by the two curves around the line x $=1.$

$($c$)$ Find the volume of the solid generated by rotating the region bounded by the two curves around the line y $=4.($The attached image is my code for $($a$)$