Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- N = 80;
- a = vector(N); a[1]=1;
- for( k = 1, N-1, { a[k+1] = k*eulerphi(a[k])+a[eulerphi(k)];});
- a;
- b = vector(N);
- for( k = 1, N, { b[k] = log(a[k]);});
- x = vector(N); for( k = 1, N, { x[k] = k;});
- plothraw(x,b)
Add Comment
Please, Sign In to add comment