Ćwiczenia

1. Obejrzyj pliki źródłowe trapez.f, calkowanie.f oraz funkcja.f. Pliki te stanowią razem program służący do obliczania całki oznaczonej funkcji, której wzór jest obecnie zapisany w module funkcja.f w przedziale od a do b metodą trapezów:
   $$\int\limits_a^b f(x) dx \approx h\frac{f(a)+f(b)}{2} + h\sum\limits_{i=1}^{n-1} f(a+ih)$$

   gdzie $n$ jest liczbą kroków całkowania a $h=(b-a)/n$. Skompiluj i uruchom program dla funkcji zakodowanej w module funkcja.f a=0, b=1 oraz 100 kroków całkowania. Sprawdź, czy program działa poprawnie.

2. Użyj procedury calkowanie.f do stabelaryzowania wartości funkcji błędu erf(x) w przedziale od 0 do 10 co 0.1. Funkcja błędu jest dana wzorem:
   $${\rm erf}(x) = \frac{2}{\sqrt{\pi}}\int\limits_0^x \exp(-\frac{y^2}{2})dy$$

   Tabelka powinna być wydrukowana w postaci sformatowanej.