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

forcal

多线程动态绘制图形是FC2D的基本特点。
参考：http://www.forcal.net/sysm/forcal9/fchtm/fc2d.htm
( e& N' x+ w. C# m5 q4 B3 W
% w6 R0 s) e- J" R2 [6 D请大家说说目前的动态绘制函数图形的库都有哪些？特点是什么？以便作为参考。

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

柯雷

这个我也不会啊啊，不过顶之

猴王再世

顶！

forcal

谢谢楼上两位！
: A3 {( b$ E4 @; k  J+ e可下载OpenFC体验一下：http://www.forcal.net/xiazai/forcal9/openfc32w.rar

liujinshan_135

加油加油加油加油加油

forcal

FC2D最近更新：Plot函数用法有变化，增加了更多选项和属性
' z0 n: d0 [$ _% [# C参考：http://www.forcal.net/sysm/forcal9/fchtm/fc2d.htm

dingpeng123

谢谢了，很有用啊！！！

forcal

绘制如下含参变量二重积分的隐函数图形：
" J$ B' Z1 U9 H
Forcal代码：
1 ^' |5 `" h8 M) Y  |0 [2 A

2. !using("fc2d","IMSL","XSLSF");
   
3. ! o9 h' M/ W8 _3 r1 E
4. ff(x,y::a,b)=y*cos(x+y*y+a)-b*x;. D  j8 f5 K# i; n# B. N  X1 X
   
6. h(x)=5*x;
   
7. + B4 @4 q5 \8 A9 b
8. g(x::b)=-2*x+b;4 n! Q1 t( s) ?; {7 ?' o6 Y9 U& y
   
10. f(bb::a,b)= b=bb, TWODQ[HFor("ff"),0,a-b+2,HFor("g"),HFor("h"),1e-6,0,6,0];' d5 Y  K- f+ b( }0 R  V
    
12. fb(fa:bb:a)= a=fa, bb=1, pqrt[HFor("f"),&bb,1e-6],bb;
    
13. # x) Y) ?/ W& B% n8 l+ C
14. ClearImslErr(),
    
15. / m+ Q) g/ w3 H5 h9 Q  w( M4 j& e
16. Plot{Iclear, Ix : -3,3, Igrid : 1,; i& t% R. a& Z- Z& y2 J* V# w+ u& F
    
18.      Iufun : HFor("fb"), Adots,300 U  b% F7 o3 d& f
    
20. };1 s) x$ d& H0 P, h& U
    

!using("fc2d","IMSL","XSLSF");<br /> <span style=\\"display:none\\">! o9 h' M/ W8 _3 r1 E</span> ff(x,y::a,b)=y*cos(x+y*y+a)-b*x;<font class=\\"jammer\\">. D  j8 f5 K# i; n# B. N  X1 X</font><br /> h(x)=5*x;<br /> <span style=\\"display:none\\">+ B4 @4 q5 \8 A9 b</span> g(x::b)=-2*x+b;<font class=\\"jammer\\">4 n! Q1 t( s) ?; {7 ?' o6 Y9 U& y</font><br /> f(bb::a,b)= b=bb, TWODQ[HFor("ff"),0,a-b+2,HFor("g"),HFor("h"),1e-6,0,6,0];<font class=\\"jammer\\">' d5 Y  K- f+ b( }0 R  V</font><br /> fb(fa:bb:a)= a=fa, bb=1, pqrt[HFor("f"),&bb,1e-6],bb;<br /> <span style=\\"display:none\\"># x) Y) ?/ W& B% n8 l+ C</span> ClearImslErr(),<br /> <span style=\\"display:none\\">/ m+ Q) g/ w3 H5 h9 Q  w( M4 j& e</span> Plot{Iclear, Ix : -3,3, Igrid : 1,<font class=\\"jammer\\">; i& t% R. a& Z- Z& y2 J* V# w+ u& F</font><br />      Iufun : HFor("fb"), Adots,30<font class=\\"jammer\\">0 U  b% F7 o3 d& f</font><br /> };<font class=\\"jammer\\">1 s) x$ d& H0 P, h& U</font><br />

图形如下：

0 |; X& V/ R; ?所有的点(a,b)如下：
. |3 p- X, m2 s; ^4 M5 _-3. , 5.677964685436554e-002
- h7 S6 e/ _5 Q-2.793103448275862 , 0.113926626758415
-2.586206896551724 , 0.1794836885445007
) R  s' L  U# B( W6 D7 O. u' O-2.379310344827586 , 0.2582958827136312
-2.172413793103448 , 0.3678246226201311
3 c+ r7 y5 m5 B( g' X-1.965517241379311 , 0.4854053136248959
$ n4 S. {5 n5 J-1.758620689655173 , 0.3109675572745501
, U% O8 M0 p4 j-1.551724137931035 , 0.4482758622854263
-1.344827586206897 , 0.6551724398984934
-1.137931034482759 , 0.8620689655216192
-0.931034482758621 , 1.06896551724139
. f! o- G3 z& x) y; t3 i-0.7241379310344831 , 0.9928357579204139
-0.5172413793103452 , 1.482758620687618
4 a4 F! Q, ^0 y-0.3103448275862073 , 1.689655172413478
$ w2 [" t# }/ y4 z% O& L9 m- m% W-0.1034482758620694 , 1.810803012030456
0.1034482758620686 , 1.762128999758652
0.3103448275862065 , 1.895805254542486
9 `3 S/ _$ J6 _7 \+ z* d9 w1 c0.5172413793103444 , 2.022877410009786
) `1 a: r; u, m0.7241379310344823 , 2.153694202098552
7 M; O& N+ C; Q+ v0.9310344827586202 , 2.320359396161155
1.137931034482758 , 2.513284341824756
1.344827586206896 , 2.7264453989435
# i4 @  d( O8 S% Y- V4 v  y1.551724137931034 , 2.925615642424742
% a: a. d( Y+ E7 y1.758620689655172 , 3.096868663763255
1.96551724137931 , 3.264359697730042
2.172413793103448 , 3.443077817907455
$ r3 q! V4 ]/ s1 ~& d2 s: W2 j2.379310344827586 , 4.4368980745608
2.586206896551724 , 4.635103786282911
0 c6 H8 z& [: T, |' ?% e2.793103448275862 , 2.815981856569838e-004
3. , 2.885091408164247e-004

forcal

图怎么看不清呢？再上传一次试试。
+ O5 W8 }  b* X. U$ @7 t楼上的问题图：

  T# z- L9 X' N函数图形：

forcal

8＃的公式决定的隐函数是个多值函数。
- H* X/ k- N6 W! g2 x/ J( x" g, E所有a=b-2的数对(a,b)都是方程的解。还有一些其他解。
; Q' @/ F4 J  B所以用隐函数绘图会更合适，难点在于计算量非常大。

1 k1 A* y) R3 H. s% @1 a以下是Forcal的隐函数绘图代码：a取-3～3，b取-2～5，绘图点数取50

2. !using("fc2d","IMSL","XSLSF");4 F4 _6 L  w& h5 A5 P, T- e
   
4. ff(x,y::a,b)=y*cos(x+y*y+a)-b*x;/ i  I! d: B+ G+ x
   
6. h(x)=5*x;
   
7. ' I* }% V9 Z5 P# g/ M6 j* o* G$ _7 |
8. g(x::b)=-2*x+b;1 u' T! G1 l1 u  e9 V# J' _& \
   
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];( t  d1 G& U1 a# z
    
12. ClearImslErr(),
    
13. . _$ w+ J$ n8 b* K& ?% V: {
14. Plot{Iclear, Ix : -3,3, Iy : -2,5, Iydynamic : 0, Igrid : 1,' g& u1 H! O& @' G: C
    
16.      Iifun : HFor("f"), Adots,50
    
17. 4 v/ B( |8 p9 h8 u3 ]7 T7 @  C# S! I
18. };
    
19. 3 I* S# S/ I8 F$ m7 N0 T0 l

!using("fc2d","IMSL","XSLSF");<font class=\\"jammer\\">4 F4 _6 L  w& h5 A5 P, T- e</font><br /> ff(x,y::a,b)=y*cos(x+y*y+a)-b*x;<font class=\\"jammer\\">/ i  I! d: B+ G+ x</font><br /> h(x)=5*x;<br /> <span style=\\"display:none\\">' I* }% V9 Z5 P# g/ M6 j* o* G$ _7 |</span> g(x::b)=-2*x+b;<font class=\\"jammer\\">1 u' T! G1 l1 u  e9 V# J' _& \</font><br /> 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\\">( t  d1 G& U1 a# z</font><br /> ClearImslErr(),<br /> <span style=\\"display:none\\">. _$ w+ J$ n8 b* K& ?% V: {</span> Plot{Iclear, Ix : -3,3, Iy : -2,5, Iydynamic : 0, Igrid : 1,<font class=\\"jammer\\">' g& u1 H! O& @' G: C</font><br />      Iifun : HFor("f"), Adots,50<br /> <span style=\\"display:none\\">4 v/ B( |8 p9 h8 u3 ]7 T7 @  C# S! I</span> };<br /> <span style=\\"display:none\\">3 I* S# S/ I8 F$ m7 N0 T0 l</span>

图形如下：
) d9 l) Z  I, U, Y
2 m' N' @: I' D9 o$ Z