load("glib")$ def taylor(N) { glib_open(); glib_window(-5,-2,5,2); glib_clear(); F = 0; for (I=0; I<=N; I++) { F=F+(-1)^I*x^(2*I+1)/fac(2*I+1); } print("sin(x) の Taylor 展開 :",0); print(2*N+1,0); print(" 次までは "); glib_line(-5,0,5,0); glib_line(0,-5,0,5); print(F); for (K=-5; K<=5; K = K+0.03) { glib_putpixel(K,subst(F,x,K)); } } print("Type in, for example, taylor(2); taylor(4);")$ end$