Risoluzione numerica dell'equazione scalare delle onde mediante ...

PDE Laboratory @ Crs4 13Figura 3: Andamento dell’errore relativo al variare di hdi Fourier di h(t), è data da:√ [πH(ω) =2 · β · 3/2 ω2 · exp −ω(4 · β) ·Il dominio computazionale è riportato in figura 6.( ) ] j+ ω(13)t 03.1 Sviluppo dell’algoritmoIl programma a cui fare riferimento per il problema è wave. Vengono ora messe in evidenzasolo le differenze rispetto al programma precedente. Dette x s ed y s le coordinate della sorgente,viene calcolato il nodo della griglia più vicino.diff=1.d90do i=1,nnodterm=(xx(i)-xs)*(xx(i)-xs)+(yy(i)-ys)*(yy(i)-ys)if (term