FORCAL图形库FC2D，多线程动态绘制图形

forcal

多线程动态绘制图形是FC2D的基本特点。
$ E% W( Z& ^7 q. S% n' {' d0 B参考：http://www.forcal.net/sysm/forcal9/fchtm/fc2d.htm
! m# i9 f. M" |% U3 s
请大家说说目前的动态绘制函数图形的库都有哪些？特点是什么？以便作为参考。

请对FC2D留下宝贵意见，以便以后改进，谢谢！

forcal

谢谢楼上两位！
可下载OpenFC体验一下：http://www.forcal.net/xiazai/forcal9/openfc32w.rar

forcal

FC2D最近更新：Plot函数用法有变化，增加了更多选项和属性
参考：http://www.forcal.net/sysm/forcal9/fchtm/fc2d.htm

forcal

绘制如下含参变量二重积分的隐函数图形：
1 i  a8 c+ U& q* \9 b
Forcal代码：
! K# n6 Y; r2 x* R4 Z

2. !using("fc2d","IMSL","XSLSF");
   
3. 3 x& S$ D) p. n
4. ff(x,y::a,b)=y*cos(x+y*y+a)-b*x;
   
5. 5 p\\" V, e3 |, h  f
6. h(x)=5*x;/ M, L5 I! G1 }6 d! Y) `5 C7 v% R
   
8. g(x::b)=-2*x+b;
   
9.   u0 @! E( k1 D5 f
10. f(bb::a,b)= b=bb, TWODQ[HFor("ff"),0,a-b+2,HFor("g"),HFor("h"),1e-6,0,6,0];
    
11. + k2 C  c, }/ s- q* n
12. fb(fa:bb:a)= a=fa, bb=1, pqrt[HFor("f"),&bb,1e-6],bb;4 Q\\" V# e& a$ ^* ?
    
14. ClearImslErr(),
    
15. % j6 i  p4 n$ z5 w4 e3 w( P0 ?
16. Plot{Iclear, Ix : -3,3, Igrid : 1,' c) U2 \/ v4 @# F) j$ O; o
    
18.      Iufun : HFor("fb"), Adots,30. K$ \4 P9 F8 S* J1 b2 j  V7 f6 \0 v
    
20. };
    
21. 0 S\\" T9 b/ q! q  E' o

!using("fc2d","IMSL","XSLSF");<br /> <span style=\\"display:none\\">3 x& S$ D) p. n</span> ff(x,y::a,b)=y*cos(x+y*y+a)-b*x;<br /> <span style=\\"display:none\\">5 p\\" V, e3 |, h  f</span> h(x)=5*x;<font class=\\"jammer\\">/ M, L5 I! G1 }6 d! Y) `5 C7 v% R</font><br /> g(x::b)=-2*x+b;<br /> <span style=\\"display:none\\">  u0 @! E( k1 D5 f</span> f(bb::a,b)= b=bb, TWODQ[HFor("ff"),0,a-b+2,HFor("g"),HFor("h"),1e-6,0,6,0];<br /> <span style=\\"display:none\\">+ k2 C  c, }/ s- q* n</span> fb(fa:bb:a)= a=fa, bb=1, pqrt[HFor("f"),&bb,1e-6],bb;<font class=\\"jammer\\">4 Q\\" V# e& a$ ^* ?</font><br /> ClearImslErr(),<br /> <span style=\\"display:none\\">% j6 i  p4 n$ z5 w4 e3 w( P0 ?</span> Plot{Iclear, Ix : -3,3, Igrid : 1,<font class=\\"jammer\\">' c) U2 \/ v4 @# F) j$ O; o</font><br />      Iufun : HFor("fb"), Adots,30<font class=\\"jammer\\">. K$ \4 P9 F8 S* J1 b2 j  V7 f6 \0 v</font><br /> };<br /> <span style=\\"display:none\\">0 S\\" T9 b/ q! q  E' o</span>

图形如下：
$ b; @' d1 {" J( A7 M0 {% V
所有的点(a,b)如下：
-3. , 5.677964685436554e-002
-2.793103448275862 , 0.113926626758415
-2.586206896551724 , 0.1794836885445007
-2.379310344827586 , 0.2582958827136312
7 Y3 b- t; w+ x3 I" ~/ s-2.172413793103448 , 0.3678246226201311
-1.965517241379311 , 0.4854053136248959
7 F! Y8 T- i/ ~) Y2 J-1.758620689655173 , 0.3109675572745501
) b8 J0 |& x; b8 o-1.551724137931035 , 0.4482758622854263
-1.344827586206897 , 0.6551724398984934
4 k; j9 I& S# S- ?6 D# s$ z# P-1.137931034482759 , 0.8620689655216192
-0.931034482758621 , 1.06896551724139
-0.7241379310344831 , 0.9928357579204139
5 A- y# f: p# I7 i7 M-0.5172413793103452 , 1.482758620687618
-0.3103448275862073 , 1.689655172413478
8 ^4 W  u2 j# Z& G  U" W-0.1034482758620694 , 1.810803012030456
0.1034482758620686 , 1.762128999758652
0.3103448275862065 , 1.895805254542486
) ^0 n" ~+ s8 J0.5172413793103444 , 2.022877410009786
# B0 D* l2 z1 t- ?( y0.7241379310344823 , 2.153694202098552
; F/ p1 G0 U, J! h0.9310344827586202 , 2.320359396161155
+ q+ Q4 f2 w, }. X: S6 s8 F1.137931034482758 , 2.513284341824756
1.344827586206896 , 2.7264453989435
8 t* E( U( |6 v2 R. `$ [, O1.551724137931034 , 2.925615642424742
% y2 P" m: C' w  h( l& I1.758620689655172 , 3.096868663763255
1.96551724137931 , 3.264359697730042
2.172413793103448 , 3.443077817907455
2.379310344827586 , 4.4368980745608
2.586206896551724 , 4.635103786282911
2.793103448275862 , 2.815981856569838e-004
0 k$ b+ x; j2 G) i1 _3. , 2.885091408164247e-004

/ S- w) X+ f! j; t

forcal

图怎么看不清呢？再上传一次试试。
# ?8 u! S! J0 O+ I楼上的问题图：
1 Z- Z9 Z" d. ~$ L. b) z2 Q
5 I9 q! f+ o+ Q- p函数图形：
" B+ B3 r9 H7 k, L  G% j- n3 |

forcal

8＃的公式决定的隐函数是个多值函数。
所有a=b-2的数对(a,b)都是方程的解。还有一些其他解。
3 P" |$ R' W3 p# N% ?% m, ^* p所以用隐函数绘图会更合适，难点在于计算量非常大。

以下是Forcal的隐函数绘图代码：a取-3～3，b取-2～5，绘图点数取50

2. !using("fc2d","IMSL","XSLSF");
   
3. 6 J0 g5 J; ?8 a$ q
4. ff(x,y::a,b)=y*cos(x+y*y+a)-b*x;( z1 i/ U1 v8 k+ |4 p
   
6. h(x)=5*x;
   
7. ' Q& @) d0 u9 X, N  o
8. g(x::b)=-2*x+b;
   
9. 1 C' G3 r0 }! t$ H* [\\" D
10. f(aa,bb::a,b)= a=aa,b=bb, TWODQ[HFor("ff"),0,a-b+2,HFor("g"),HFor("h"),1e-6,0,6,0];/ u9 f! o# E6 ~\\" N, w. M/ f$ R
    
12. ClearImslErr(),
    
13. 5 q+ k\\" r/ d7 V4 y
14. Plot{Iclear, Ix : -3,3, Iy : -2,5, Iydynamic : 0, Igrid : 1,\\" O7 O1 `  ]# f5 W/ ~% E8 X
    
16.      Iifun : HFor("f"), Adots,50
    
17. / K\\" ?' u) t$ d% o7 R: d1 P
18. };
    
19. 0 r3 _, \$ U% X: c, u6 c( W8 C: [4 T

!using("fc2d","IMSL","XSLSF");<br /> <span style=\\"display:none\\">6 J0 g5 J; ?8 a$ q</span> ff(x,y::a,b)=y*cos(x+y*y+a)-b*x;<font class=\\"jammer\\">( z1 i/ U1 v8 k+ |4 p</font><br /> h(x)=5*x;<br /> <span style=\\"display:none\\">' Q& @) d0 u9 X, N  o</span> g(x::b)=-2*x+b;<br /> <span style=\\"display:none\\">1 C' G3 r0 }! t$ H* [\\" D</span> f(aa,bb::a,b)= a=aa,b=bb, TWODQ[HFor("ff"),0,a-b+2,HFor("g"),HFor("h"),1e-6,0,6,0];<font class=\\"jammer\\">/ u9 f! o# E6 ~\\" N, w. M/ f$ R</font><br /> ClearImslErr(),<br /> <span style=\\"display:none\\">5 q+ k\\" r/ d7 V4 y</span> Plot{Iclear, Ix : -3,3, Iy : -2,5, Iydynamic : 0, Igrid : 1,<font class=\\"jammer\\">\\" O7 O1 `  ]# f5 W/ ~% E8 X</font><br />      Iifun : HFor("f"), Adots,50<br /> <span style=\\"display:none\\">/ K\\" ?' u) t$ d% o7 R: d1 P</span> };<br /> <span style=\\"display:none\\">0 r3 _, \$ U% X: c, u6 c( W8 C: [4 T</span>

图形如下：