// パラメータ定義
int maxiter = chi("iterations");
float scale = 1 / chf("scale");
float offset_x = chf("offset_x");
float offset_z = chf("offset_z");
float escape_threshold = chf("escape_threshold");
float amplifier = chf("amplifier");
float core_height = chf("core_height");

// メインルーチンの'mandelset'
function float mandelset(float x0, y0; int imax; float escape) {
    float x = 0;
    float y = 0;
    float xnew, ynew;
    int i;

    for (i = 0; i < imax; i++) {
        xnew = x * x - y * y + x0;
        ynew = 2 * x * y + y0;
        if (xnew * xnew + ynew * ynew > escape) {
            return float(i) / float(imax); // 正規化された反復深さ
        }
        x = xnew;
        y = ynew;
    }
    return -1.0; // マンデルブロ集合内のポイント
}

// マンデルブロ集合の計算
float value = mandelset(@P.x * scale - offset_x, @P.z * scale - offset_z, maxiter, escape_threshold);

if (value >= 0) {
    @height = value * amplifier;
    @height += sin(value * 10) * 0.1; // Bands
} 

else {
    @height = core_height;
}