#! /usr/bin/gnuplot a = 2 ## trocando este valor voce pode obter a parĀ“bola em outros pontos. f(x) = 2*sin(3*x + 4) ; rho = 0.00001; Q_rho_f(x) = (f(x+rho)-f(x))/rho; # operador quociente Q2_rho_f(x) = (Q_rho_f(x+rho) - Q_rho_f(x))/rho # op. quoc. de segunda ordem a0 = f(a); a1 = Q_rho_f(a); a2 = Q2_rho_f(a); P(x) = a0 + a1*(x-a) + 0.5*a2*(x-a)**2; set xrange [0:3] set yrange [-3:3] plot f(x),0, P(x) pause -2