Presenteremo un metodo, associato al nome di Sergey A. Smolyak, per l'interpolazione di funzioni reali definite su un ipercubo d-dimensionale. Se la funzione da interpolare ha una certa struttura tensoriale, allora il metodo di Smolyak permette di costruire, partendo da uno schema di interpolazione uni-dimensionale, un algoritmo di interpolazione a dimensione d che mantiene l'ordine di approssimazione dello schema di partenza, a parte un fattore logaritmico il quale può contenere una dipendenza esponenziale dalla dimensione. L'effetto della "maledizione della dimensionalità" viene in questo modo attenuato per certe classi di funzioni. I punti in cui occorre conoscere i valori dell'interpolanda si trovano su "griglie sparse" anzichè su griglie uniformi. Faremo alcuni commenti sull'applicabilità dell'algoritmo alla programmazione dinamica per la soluzione numerica di problemi di controllo ottimo a tempo continuo.