Using any arbitrary-precision library will be much slower than floating-point numbers, especially for something like calculating the Mandelbrot set which does many, many repeated iterations. If you want speed I'd recommend rewriting this in C using a library like gmplib.
The issue in your code is that you are using a bcadd
instead of a bcmul
for your first line inside the loop.