Esercitazione 1b

L'esercitazione consiste nello sviluppo dei punti seguenti.

1. Scrivere un programma FORTRAN che implementa il codice FEM descritto in precedenza per l'equazione conservazione_massa.
2. Applicare il codice al problema esercitazione utilizzando una triangolazione generata dal codice RETRAG avendo posto R=1 come valore del raggio esterno, e valori opportuni per i parametri r, , e NL.
3. Verificare il codice di calcolo confrontando in acluni punti significativi la soluzione numerica con quella analitica.
4. Variare i parametri e NL nella costruzione della griglia e calcolare su alcuni punti significativi l'errore commesso nella soluzione numerica, calcolato come differenza in valore assoluto tra la soluzione numerica e quella analitica.
5. Costruire un diagramma in cui, in un punto prescelto, si grafica l'errore in funzione di NL per un valore di costante opportunamente scelto.