mandel.js 1.2 KB

1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253
  1. /*
  2. * Mandelbrot example:
  3. *
  4. * $ ./duk mandel.js
  5. * [...]
  6. */
  7. function mandel() {
  8. var w = 76, h = 28, iter = 100;
  9. var i, j, k, c;
  10. var x0, y0, xx, yy, xx2, yy2;
  11. var line;
  12. for (i = 0; i < h; i++) {
  13. y0 = (i / h) * 2.5 - 1.25;
  14. for (j = 0, line = []; j < w; j++) {
  15. x0 = (j / w) * 3.0 - 2.0;
  16. for (k = 0, xx = 0, yy = 0, c = '#'; k < iter; k++) {
  17. /* z -> z^2 + c
  18. * -> (xx+i*yy)^2 + (x0+i*y0)
  19. * -> xx*xx+i*2*xx*yy-yy*yy + x0 + i*y0
  20. * -> (xx*xx - yy*yy + x0) + i*(2*xx*yy + y0)
  21. */
  22. xx2 = xx*xx; yy2 = yy*yy;
  23. if (xx2 + yy2 < 4.0) {
  24. yy = 2*xx*yy + y0;
  25. xx = xx2 - yy2 + x0;
  26. } else {
  27. /* xx^2 + yy^2 >= 4.0 */
  28. if (k < 3) { c = '.'; }
  29. else if (k < 5) { c = ','; }
  30. else if (k < 10) { c = '-'; }
  31. else { c = '='; }
  32. break;
  33. }
  34. }
  35. line.push(c);
  36. }
  37. print(line.join(''));
  38. }
  39. }
  40. try {
  41. mandel();
  42. } catch (e) {
  43. print(e.stack || e);
  44. }