Beh, senza sapere cosa deve fare lo script è difficile poterti aiutare concretamente, l'unica cosa che ti posso dire per ora è che la riga incriminata da errore perché assegni ad ogni elemento del vettore u un altro vettore (e la cosa ovviamente non la puoi fare, visto che ogni elemento di un vettore preso in se è uno scalare).
Se vuoi invece che u contenga una matrice allora devi mettere:

u(:,z)=(X-z)./h;

in questo modo avrai una matrice di n righe (73) per 10 colonne.
cambiato questo però ci sono altri errori disseminati per lo script che non saprei come correggere senza sapere quale deve essere il risultato finale.

Ciao e auguri