Modifier le nom et ajouter les arguments nécessaires à la liste d’appel du sous-programme : Nommons par exemple ce sous-programme ondes1DP1 et modifions les arguments d’appel :
conservons
,
,
,
,
,
,
,
et
,
et ajoutons
c, la vitesse d’onde,
V, le vecteur de la vitesse initiale
aux noeuds
et au temps
.
ce qui donne U=ondes1DP1( nx,x,fs,cldi,clne,c,theta ,dt,nt,Uo,V)