# source("http://macosa.dima.unige.it/r.R")
# Questa e' una "funzione di Bessel del I tipo". Non ci preoccupiamo del suo signficato
B = function(x) besselJ(x,1)
# Consideriamo la serie di Fourier di cui i seguenti sono polinomi di grado n
p = function(x) { f = pi/4; for(i in 1:n) f = f+(-1)^i*B(i*pi)/i*cos(i*pi*x); f}
# Si vede dai grafici seguenti che essi tendono a F che ha per grafico un semicerchio
F = function(x) sqrt(1-(x-1)^2)
BF=5; HF=3
PLANE(0,2, 0,1.2); graph2(F, 0,2, "brown")
n = 1; graph1(p,0,2,"blue")
n = 2; graph1(p,0,2,"red")
n = 3; graph1(p,0,2,"seagreen")
n = 4; graph1(p,0,2,"brown")
n =200; graph1(p,0,2,"orange")