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

forcal

多线程动态绘制图形是FC2D的基本特点。
0 k0 h+ Q0 D1 o6 o+ y参考：http://www.forcal.net/sysm/forcal9/fchtm/fc2d.htm
. B( w% Y5 w8 B7 [( Z+ v
& H/ o- r% }* s& U请大家说说目前的动态绘制函数图形的库都有哪些？特点是什么？以便作为参考。
1 Y2 o1 A9 B+ |9 w6 a8 b+ Y. ?
请对FC2D留下宝贵意见，以便以后改进，谢谢！

柯雷

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

猴王再世

顶！

forcal

谢谢楼上两位！
5 E% N- A8 e: A8 }可下载OpenFC体验一下：http://www.forcal.net/xiazai/forcal9/openfc32w.rar

liujinshan_135

加油加油加油加油加油

forcal

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

dingpeng123

谢谢了，很有用啊！！！

forcal

绘制如下含参变量二重积分的隐函数图形：
( D7 m- Y. e2 X  ?& W6 Q
Forcal代码：
5 k. _1 |8 w# k& N) l, u9 X

2. !using("fc2d","IMSL","XSLSF");+ _9 |, ~, g$ A( ^
   
4. ff(x,y::a,b)=y*cos(x+y*y+a)-b*x;$ W, l1 |; o$ n4 }; R7 w+ E$ E
   
6. h(x)=5*x;
   
7. 7 c\\" P  O: C0 i3 H8 x
8. g(x::b)=-2*x+b;# M( P1 @. e. ]1 p& S6 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. 7 `2 O; |: Q  A) _# h
12. fb(fa:bb:a)= a=fa, bb=1, pqrt[HFor("f"),&bb,1e-6],bb;. b' Z9 t5 f) }0 ^( y* A$ p
    
14. ClearImslErr(),
    
15. 8 h# F* s2 P1 F
16. Plot{Iclear, Ix : -3,3, Igrid : 1,
    
17. / a8 \' L8 Q- N5 _
18.      Iufun : HFor("fb"), Adots,30
    
19. $ y8 D\\" [& _' o+ n4 p1 o0 A4 q* z
20. };
    
21. % I/ ?& R  Q8 |/ K% T' F# ?

!using("fc2d","IMSL","XSLSF");<font class=\\"jammer\\">+ _9 |, ~, g$ A( ^</font><br /> ff(x,y::a,b)=y*cos(x+y*y+a)-b*x;<font class=\\"jammer\\">$ W, l1 |; o$ n4 }; R7 w+ E$ E</font><br /> h(x)=5*x;<br /> <span style=\\"display:none\\">7 c\\" P  O: C0 i3 H8 x</span> g(x::b)=-2*x+b;<font class=\\"jammer\\"># M( P1 @. e. ]1 p& S6 F. `& _</font><br /> 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\\">7 `2 O; |: Q  A) _# h</span> fb(fa:bb:a)= a=fa, bb=1, pqrt[HFor("f"),&bb,1e-6],bb;<font class=\\"jammer\\">. b' Z9 t5 f) }0 ^( y* A$ p</font><br /> ClearImslErr(),<br /> <span style=\\"display:none\\">8 h# F* s2 P1 F</span> Plot{Iclear, Ix : -3,3, Igrid : 1,<br /> <span style=\\"display:none\\">/ a8 \' L8 Q- N5 _</span>      Iufun : HFor("fb"), Adots,30<br /> <span style=\\"display:none\\">$ y8 D\\" [& _' o+ n4 p1 o0 A4 q* z</span> };<br /> <span style=\\"display:none\\">% I/ ?& R  Q8 |/ K% T' F# ?</span>

图形如下：
( v$ S/ R4 |) R% P
所有的点(a,b)如下：
-3. , 5.677964685436554e-002
; i& m2 ]) f4 Z4 p* R-2.793103448275862 , 0.113926626758415
-2.586206896551724 , 0.1794836885445007
4 B: P+ b) D2 @4 E4 A-2.379310344827586 , 0.2582958827136312
-2.172413793103448 , 0.3678246226201311
' W( p; ]! ~4 E- F& f. [+ M% {-1.965517241379311 , 0.4854053136248959
, [1 M- C7 s; T# {3 }-1.758620689655173 , 0.3109675572745501
- S! U% Q" d8 e! A7 S-1.551724137931035 , 0.4482758622854263
! p+ k/ G8 g! b4 O- }-1.344827586206897 , 0.6551724398984934
. G3 d. a% V4 d8 V) c! o& j8 e-1.137931034482759 , 0.8620689655216192
' ?) d6 B) t* B# m5 x-0.931034482758621 , 1.06896551724139
-0.7241379310344831 , 0.9928357579204139
-0.5172413793103452 , 1.482758620687618
-0.3103448275862073 , 1.689655172413478
$ w8 D, U0 k; a-0.1034482758620694 , 1.810803012030456
  u: {6 K" a/ n- p0.1034482758620686 , 1.762128999758652
0.3103448275862065 , 1.895805254542486
0.5172413793103444 , 2.022877410009786
0.7241379310344823 , 2.153694202098552
0.9310344827586202 , 2.320359396161155
1.137931034482758 , 2.513284341824756
1.344827586206896 , 2.7264453989435
1.551724137931034 , 2.925615642424742
1.758620689655172 , 3.096868663763255
1.96551724137931 , 3.264359697730042
" k% X$ `  V- m8 ?5 Q2.172413793103448 , 3.443077817907455
2.379310344827586 , 4.4368980745608
2.586206896551724 , 4.635103786282911
0 s. W( N. I1 v; \2.793103448275862 , 2.815981856569838e-004
3. , 2.885091408164247e-004
% C1 C' z: n7 z" F0 M6 j

forcal

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

函数图形：

( x2 \1 [5 M' k" }

forcal

8＃的公式决定的隐函数是个多值函数。
* B2 w1 U5 q) q3 g" h1 `; C所有a=b-2的数对(a,b)都是方程的解。还有一些其他解。
所以用隐函数绘图会更合适，难点在于计算量非常大。
9 C6 H0 P6 Y$ g) c
+ {% H, T1 `9 P( @9 Y以下是Forcal的隐函数绘图代码：a取-3～3，b取-2～5，绘图点数取50
# E9 K( D  S2 g3 j9 |3 V' \! \& }+ o

2. !using("fc2d","IMSL","XSLSF");
   
3. ) Y) L+ e- }& B
4. ff(x,y::a,b)=y*cos(x+y*y+a)-b*x;
   
5. $ E2 |) `6 D- s1 t0 x+ b5 h
6. h(x)=5*x;  `; y% B7 M/ N8 W$ |4 G
   
8. g(x::b)=-2*x+b;% i5 _7 i& f1 A. ?' ]: N: }. N& R$ l
   
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];
    
11. % V$ Q# I5 l4 z# @# S4 W* }
12. ClearImslErr(),  H& y4 x9 R% f) Q; a, U8 n5 m
    
14. Plot{Iclear, Ix : -3,3, Iy : -2,5, Iydynamic : 0, Igrid : 1,
    
15. 6 n$ ]4 E; X0 V& W; j1 P' g$ K
16.      Iifun : HFor("f"), Adots,50
    
17. 1 Q9 B! \' {* d$ u# `
18. };/ j! l- d$ r; _0 n- U3 ~8 i* X
    

!using("fc2d","IMSL","XSLSF");<br /> <span style=\\"display:none\\">) Y) L+ e- }& B</span> ff(x,y::a,b)=y*cos(x+y*y+a)-b*x;<br /> <span style=\\"display:none\\">$ E2 |) `6 D- s1 t0 x+ b5 h</span> h(x)=5*x;<font class=\\"jammer\\">  `; y% B7 M/ N8 W$ |4 G</font><br /> g(x::b)=-2*x+b;<font class=\\"jammer\\">% i5 _7 i& f1 A. ?' ]: N: }. N& R$ l</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];<br /> <span style=\\"display:none\\">% V$ Q# I5 l4 z# @# S4 W* }</span> ClearImslErr(),<font class=\\"jammer\\">  H& y4 x9 R% f) Q; a, U8 n5 m</font><br /> Plot{Iclear, Ix : -3,3, Iy : -2,5, Iydynamic : 0, Igrid : 1,<br /> <span style=\\"display:none\\">6 n$ ]4 E; X0 V& W; j1 P' g$ K</span>      Iifun : HFor("f"), Adots,50<br /> <span style=\\"display:none\\">1 Q9 B! \' {* d$ u# `</span> };<font class=\\"jammer\\">/ j! l- d$ r; _0 n- U3 ~8 i* X</font><br />

图形如下：
) J" |, _' q  j& i; ]; w4 F
) p) ]3 S+ K+ ~; E0 j* u2 M0 ^- j; q* @
7 r- e8 h7 M) Z6 J! L1 n