import sys from math import exp n=int(sys.argv[1]) # Orden av polynom xmax=float(sys.argv[2]) # Regner fra -xmax xmax npunkt=int(2*xmax/0.1) # ca. antall punkter med dx=0.1 if(npunkt<=2): npunkt=2 #sikrer minimum 2 punkter dx=2*xmax/npunkt c=[1.0]; i=1 while (i <=n): c.append(c[i-1]/i) # blir 1 over fakultet til i i=i+1 i=0 while (i <=npunkt): x=i*dx-xmax; p=c[n]; k=n-1 while(k>=0): p=c[k]+x*p k=k-1 fasit=exp(x); diff=fasit-p print x,p,fasit,diff i=i+1