# 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")