...probably some differences in handling of NAN / INF, etc. The algorithm makes your $d and $e values overflow rather quickly.

By reducing the number of iterations from 99 to something like 9 you should get a coarse approximation of the desired shape -- on both platforms.