This textbook is intended to supplement the classical theory of uni- and multivariate splines and their approximation and interpolation properties with those of fractals, fractal functions, and fractal surfaces. This synthesis will complement currently required courses dealing with these topics and expose the prospective reader to some new and deep relationships. In addition to providing a classical introduction to the main issues involving approximation and interpolation with uni- and multivariate splines, cardinal and exponential splines, and their connection to wavelets and multiscale analysis, which comprises the first half of the book, the second half will describe fractals, fractal functions and fractal surfaces, and their properties. This also includes the new burgeoning theory of superfractals and superfractal functions. The theory of splines is well-established but the relationship to fractal functions is novel. Throughout the book, connections between these two apparently different areas will be exposed and presented. In this way, more options are given to the prospective reader who will encounter complex approximation and interpolation problems in real-world modeling. Numerous examples, figures, and exercises accompany the material.