tubeThickness:=0.05
tubeStretchFactor:=10

graph:=Plot3D[Gamma[2+1]/Gamma[2-(2-v)+1]*Gamma[Ceiling[k]+2-(2-v)]/Gamma[Ceiling[k]],{v,0,2},{k,0,7}, Filling->Bottom,Boxed->False,PlotRange->{0,60},PlotTheme->"FilledSurface",Exclusions->None,MaxRecursion->8]


constantPlot:=Plot3D[2+tubeStretchFactor*Sqrt[(tubeThickness/2)^2-x^2],{x,0,tubeThickness/2},{y,0,7},PlotTheme->"FilledSurface"]
linearPlot:=Plot3D[2*Ceiling[y]+tubeStretchFactor*Sqrt[(tubeThickness/2)^2-(x-1)^2],{x,1-tubeThickness/2,1+tubeThickness/2},{y,0,7},PlotTheme->"FilledSurface"]
quadraticPlot:=Plot3D[rising[Ceiling[y],2]+tubeStretchFactor*Sqrt[(tubeThickness/2)^2-(x-2)^2],{x,2-tubeThickness/2,2},{y,0,7},PlotTheme->"FilledSurface"]

output=Show[graph,constantPlot,linearPlot,quadraticPlot]
Export[NotebookDirectory[]<>"Fractional Discrete Quadratic.stl",output]