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

forcal

多线程动态绘制图形是FC2D的基本特点。
参考：http://www.forcal.net/sysm/forcal9/fchtm/fc2d.htm

5 |3 T0 I+ H3 F  M请大家说说目前的动态绘制函数图形的库都有哪些？特点是什么？以便作为参考。

! E+ _4 c5 p+ |! ~$ ]! H请对FC2D留下宝贵意见，以便以后改进，谢谢！

柯雷

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

猴王再世

顶！

forcal

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

liujinshan_135

加油加油加油加油加油

forcal

FC2D最近更新：Plot函数用法有变化，增加了更多选项和属性
, q# j( O! x1 G- F& H参考：http://www.forcal.net/sysm/forcal9/fchtm/fc2d.htm

dingpeng123

谢谢了，很有用啊！！！

forcal

绘制如下含参变量二重积分的隐函数图形：
8 R+ Z/ w( y2 x+ @
Forcal代码：
( P0 `4 |8 [! D0 E% ~; L

2. !using("fc2d","IMSL","XSLSF");5 `' {2 Z, c( \2 n
   
4. ff(x,y::a,b)=y*cos(x+y*y+a)-b*x;0 v& Z) `& g/ f* N
   
6. h(x)=5*x;1 }3 C$ S: ~7 Z\\" n1 U
   
8. g(x::b)=-2*x+b;
   
9. 6 _! e: P& B$ `3 {2 `$ O1 G9 x
10. f(bb::a,b)= b=bb, TWODQ[HFor("ff"),0,a-b+2,HFor("g"),HFor("h"),1e-6,0,6,0];
    
11.   B0 Y$ }. q2 s; t
12. fb(fa:bb:a)= a=fa, bb=1, pqrt[HFor("f"),&bb,1e-6],bb;
    
13. ! @6 C$ j- |) }+ [' i' C: b' W( p3 u
14. ClearImslErr(),% J$ Y. m; a/ W8 s' d
    
16. Plot{Iclear, Ix : -3,3, Igrid : 1,* y- K/ U+ v+ |% A
    
18.      Iufun : HFor("fb"), Adots,30
    
19. . _3 r) O0 i8 E0 D0 e0 M4 O7 T
20. };
    
21. 2 N2 C' K; V1 }7 x) `

!using("fc2d","IMSL","XSLSF");<font class=\\"jammer\\">5 `' {2 Z, c( \2 n</font><br /> ff(x,y::a,b)=y*cos(x+y*y+a)-b*x;<font class=\\"jammer\\">0 v& Z) `& g/ f* N</font><br /> h(x)=5*x;<font class=\\"jammer\\">1 }3 C$ S: ~7 Z\\" n1 U</font><br /> g(x::b)=-2*x+b;<br /> <span style=\\"display:none\\">6 _! e: P& B$ `3 {2 `$ O1 G9 x</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\\">  B0 Y$ }. q2 s; t</span> fb(fa:bb:a)= a=fa, bb=1, pqrt[HFor("f"),&bb,1e-6],bb;<br /> <span style=\\"display:none\\">! @6 C$ j- |) }+ [' i' C: b' W( p3 u</span> ClearImslErr(),<font class=\\"jammer\\">% J$ Y. m; a/ W8 s' d</font><br /> Plot{Iclear, Ix : -3,3, Igrid : 1,<font class=\\"jammer\\">* y- K/ U+ v+ |% A</font><br />      Iufun : HFor("fb"), Adots,30<br /> <span style=\\"display:none\\">. _3 r) O0 i8 E0 D0 e0 M4 O7 T</span> };<br /> <span style=\\"display:none\\">2 N2 C' K; V1 }7 x) `</span>

图形如下：

/ W% u  n/ u; `7 G4 _9 `所有的点(a,b)如下：
8 C1 n9 y1 h9 o9 ^8 E/ ~" p2 W-3. , 5.677964685436554e-002
2 L0 J0 s7 M; r8 r. e3 H: V5 z9 {-2.793103448275862 , 0.113926626758415
5 Z, J, v7 A8 \' [% d-2.586206896551724 , 0.1794836885445007
-2.379310344827586 , 0.2582958827136312
0 r; U6 }0 p2 Z-2.172413793103448 , 0.3678246226201311
-1.965517241379311 , 0.4854053136248959
-1.758620689655173 , 0.3109675572745501
5 P) L/ Q5 T$ Y3 `3 v3 U% e# d-1.551724137931035 , 0.4482758622854263
-1.344827586206897 , 0.6551724398984934
-1.137931034482759 , 0.8620689655216192
5 f2 X* z* ?, f1 _1 l: Y-0.931034482758621 , 1.06896551724139
-0.7241379310344831 , 0.9928357579204139
-0.5172413793103452 , 1.482758620687618
-0.3103448275862073 , 1.689655172413478
-0.1034482758620694 , 1.810803012030456
0.1034482758620686 , 1.762128999758652
0.3103448275862065 , 1.895805254542486
5 Y; {8 C2 x) i0.5172413793103444 , 2.022877410009786
0.7241379310344823 , 2.153694202098552
& @* X5 j) R' v4 l0 }0.9310344827586202 , 2.320359396161155
$ s. f7 ~9 W* }0 m. l9 g+ P5 n1.137931034482758 , 2.513284341824756
" e: ?; m3 F: I8 T1.344827586206896 , 2.7264453989435
" N) j1 [  {: s1.551724137931034 , 2.925615642424742
; T% v. P7 J4 h  G7 o6 X% n1.758620689655172 , 3.096868663763255
1.96551724137931 , 3.264359697730042
2.172413793103448 , 3.443077817907455
2.379310344827586 , 4.4368980745608
  J9 R& }0 G- T4 m  \; E2.586206896551724 , 4.635103786282911
2.793103448275862 , 2.815981856569838e-004
, A9 c* f3 J$ Q; ]3. , 2.885091408164247e-004
1 N) o0 r, @3 N" n! o

forcal

图怎么看不清呢？再上传一次试试。
楼上的问题图：

0 u- m% @! U: p5 R( [: }函数图形：

$ c6 Q) G. I" \0 z: [4 F

forcal

8＃的公式决定的隐函数是个多值函数。
所有a=b-2的数对(a,b)都是方程的解。还有一些其他解。
" P( \1 p: V* l+ V% ?/ p所以用隐函数绘图会更合适，难点在于计算量非常大。
+ c3 k" r8 x8 ^; B* O& a
以下是Forcal的隐函数绘图代码：a取-3～3，b取-2～5，绘图点数取50
! k, p1 v, }1 f, A; G

2. !using("fc2d","IMSL","XSLSF");
   
3. 6 u* Y) k4 _+ }9 y
4. ff(x,y::a,b)=y*cos(x+y*y+a)-b*x;
   
5. 1 ?' H\\" h/ z9 w6 d
6. h(x)=5*x;
   
7. $ z: V- c3 L! M7 C; B' c' r
8. g(x::b)=-2*x+b;
   
9. % `1 C5 O, ^; Q8 I5 m4 V5 q
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];' n* ^4 q0 g% l4 D8 K
    
12. ClearImslErr(),
    
13. : b2 d% k6 l6 `& Q) S
14. Plot{Iclear, Ix : -3,3, Iy : -2,5, Iydynamic : 0, Igrid : 1,
    
15. 8 l6 c! S9 r, B# l$ z
16.      Iifun : HFor("f"), Adots,50* M2 m: x8 {! ^+ G9 x* n2 Y% A
    
18. };
    
19. . C: s( Y/ N5 x9 m# m( h

!using("fc2d","IMSL","XSLSF");<br /> <span style=\\"display:none\\">6 u* Y) k4 _+ }9 y</span> ff(x,y::a,b)=y*cos(x+y*y+a)-b*x;<br /> <span style=\\"display:none\\">1 ?' H\\" h/ z9 w6 d</span> h(x)=5*x;<br /> <span style=\\"display:none\\">$ z: V- c3 L! M7 C; B' c' r</span> g(x::b)=-2*x+b;<br /> <span style=\\"display:none\\">% `1 C5 O, ^; Q8 I5 m4 V5 q</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\\">' n* ^4 q0 g% l4 D8 K</font><br /> ClearImslErr(),<br /> <span style=\\"display:none\\">: b2 d% k6 l6 `& Q) S</span> Plot{Iclear, Ix : -3,3, Iy : -2,5, Iydynamic : 0, Igrid : 1,<br /> <span style=\\"display:none\\">8 l6 c! S9 r, B# l$ z</span>      Iifun : HFor("f"), Adots,50<font class=\\"jammer\\">* M2 m: x8 {! ^+ G9 x* n2 Y% A</font><br /> };<br /> <span style=\\"display:none\\">. C: s( Y/ N5 x9 m# m( h</span>

图形如下：

1 R1 C* y# u! _7 n2 n
/ ^) J9 V  E# l' l6 U! L. m$ ]