level 13
LNSZDZG
楼主
在画板论坛里学习了xklppp老师的《几何画板下的编程思考:javascript》,感触颇深,不敢独享,(没征求xklppp老师的意见)先发在这里。
UF代码:
Julia && Mandelbrot {
;这是根据xklppp老师的程序写的一个综合程序
; javascript 程序可以参考如下网页
; http://www.inrm3d.cn/viewthread.php?tid=4831&extra=page%3D1
init:
; 程序初始化
complex z = #pixel
complex c = @seed
; 分形类型选择
if @type==1
z = @start
c = #pixel
endif
;变换模块1 "c" "1/c" "c-1/c" "c+1/c" "c^2" "1/c^2" "c-1/c^2" "c+1/c^2" "√c" "1/√c" "c-1/√c" "c+1/√c"
if @tran==1
c = 1/c
elseif @tran==2
c = c-1/c
elseif @tran==3
c = c+1/c
elseif @tran==4
c = c*c
elseif @tran==5
c = 1/c^2
elseif @tran==6
c = c-1/c^2
elseif @tran==7
c = c+1/c^2
elseif @tran==8
c = sqrt(c)
elseif @tran==9
c = 1/sqrt(c)
elseif @tran==10
c =c-1/sqrt(c)
elseif @tran==11
c = c+1/sqrt(c)
elseif @tran==12
c = log(c^3)/3
endif
;变换模块2
if @tran2==1
z = 1/z
elseif @tran2==2
z = z-1/z
elseif @tran2==3
z = z+1/z
elseif @tran2==4
z = z*z
elseif @tran2==5
z = 1/z^2
elseif @tran2==6
z = z-1/z^2
elseif @tran2==7
z = z+1/z^2
elseif @tran2==8
z = sqrt(z)
elseif @tran2==9
z = 1/sqrt(z)
elseif @tran2==10
z =z-1/sqrt(z)
elseif @tran2==11
z = z+1/sqrt(z)
endif
;变换模块3
c = @fn5(@fn4(@fn3(@fn1(c))))
;变换模块4
z = @fn8(@fn7(@fn6(@fn2(z))))
;循环模块
loop:
z = z^2 +c
z = @fn9(z)
;变量说明模块
bailout:
|z| < @bailout
default:
title = "Julia && Mandelbrot"
maxiter = 200
periodicity = 0
center = (0.0,0.0)
magn = 1
param seed
caption = "seed"
default = (-0.93,0.305)
endparam
param type
caption="Set type"
default=1
enum="Julia" "Mandelbrot"
endparam
param tran
caption="trans of c"
default=0
enum="c" "1/c" "c-1/c" "c+1/c" "c^2" "1/c^2" "c-1/c^2" "c+1/c^2" "√c" "1/√c" "c-1/√c" "c+1/√c" "log(c^3)/3"
endparam
param tran2
caption="trans of z"
default=0
enum="z" "1/z" "z-1/z" "z+1/z" "z^2" "1/z^2" "z-1/z^2" "z+1/z^2" "√z" "1/√z" "z-1/√z" "z+1/√z"
endparam
func fn1
caption="function1 of c"
default= ident()
endfunc
func fn3
caption="function2 of c"
default= ident()
endfunc
func fn4
caption="function3 of c"
default= ident()
endfunc
func fn5
caption="function4 of c"
default= ident()
endfunc
func fn2
caption="function1 of z"
default= ident()
endfunc
func fn6
caption="function2 of z"
default= ident()
endfunc
func fn7
caption="function3 of z"
default= ident()
endfunc
func fn8
caption="function4 of z"
default= ident()
endfunc
func fn9
caption="function of z"
default= ident()
endfunc
param start
caption = "start"
default = (0,0.00)
endparam
param bailout
caption = "Bailout"
default = 10E5
endparam
}
MC工作表过会分享。
2014年11月06日 08点11分
1
UF代码:
Julia && Mandelbrot {
;这是根据xklppp老师的程序写的一个综合程序
; javascript 程序可以参考如下网页
; http://www.inrm3d.cn/viewthread.php?tid=4831&extra=page%3D1
init:
; 程序初始化
complex z = #pixel
complex c = @seed
; 分形类型选择
if @type==1
z = @start
c = #pixel
endif
;变换模块1 "c" "1/c" "c-1/c" "c+1/c" "c^2" "1/c^2" "c-1/c^2" "c+1/c^2" "√c" "1/√c" "c-1/√c" "c+1/√c"
if @tran==1
c = 1/c
elseif @tran==2
c = c-1/c
elseif @tran==3
c = c+1/c
elseif @tran==4
c = c*c
elseif @tran==5
c = 1/c^2
elseif @tran==6
c = c-1/c^2
elseif @tran==7
c = c+1/c^2
elseif @tran==8
c = sqrt(c)
elseif @tran==9
c = 1/sqrt(c)
elseif @tran==10
c =c-1/sqrt(c)
elseif @tran==11
c = c+1/sqrt(c)
elseif @tran==12
c = log(c^3)/3
endif
;变换模块2
if @tran2==1
z = 1/z
elseif @tran2==2
z = z-1/z
elseif @tran2==3
z = z+1/z
elseif @tran2==4
z = z*z
elseif @tran2==5
z = 1/z^2
elseif @tran2==6
z = z-1/z^2
elseif @tran2==7
z = z+1/z^2
elseif @tran2==8
z = sqrt(z)
elseif @tran2==9
z = 1/sqrt(z)
elseif @tran2==10
z =z-1/sqrt(z)
elseif @tran2==11
z = z+1/sqrt(z)
endif
;变换模块3
c = @fn5(@fn4(@fn3(@fn1(c))))
;变换模块4
z = @fn8(@fn7(@fn6(@fn2(z))))
;循环模块
loop:
z = z^2 +c
z = @fn9(z)
;变量说明模块
bailout:
|z| < @bailout
default:
title = "Julia && Mandelbrot"
maxiter = 200
periodicity = 0
center = (0.0,0.0)
magn = 1
param seed
caption = "seed"
default = (-0.93,0.305)
endparam
param type
caption="Set type"
default=1
enum="Julia" "Mandelbrot"
endparam
param tran
caption="trans of c"
default=0
enum="c" "1/c" "c-1/c" "c+1/c" "c^2" "1/c^2" "c-1/c^2" "c+1/c^2" "√c" "1/√c" "c-1/√c" "c+1/√c" "log(c^3)/3"
endparam
param tran2
caption="trans of z"
default=0
enum="z" "1/z" "z-1/z" "z+1/z" "z^2" "1/z^2" "z-1/z^2" "z+1/z^2" "√z" "1/√z" "z-1/√z" "z+1/√z"
endparam
func fn1
caption="function1 of c"
default= ident()
endfunc
func fn3
caption="function2 of c"
default= ident()
endfunc
func fn4
caption="function3 of c"
default= ident()
endfunc
func fn5
caption="function4 of c"
default= ident()
endfunc
func fn2
caption="function1 of z"
default= ident()
endfunc
func fn6
caption="function2 of z"
default= ident()
endfunc
func fn7
caption="function3 of z"
default= ident()
endfunc
func fn8
caption="function4 of z"
default= ident()
endfunc
func fn9
caption="function of z"
default= ident()
endfunc
param start
caption = "start"
default = (0,0.00)
endparam
param bailout
caption = "Bailout"
default = 10E5
endparam
}
MC工作表过会分享。

