require 'plot numeric trig viewmat' NB. Usage: NB. * arua3d 3 NB. * fun_01 plot_aqua 20 NB. ----------class aquarium------------- NB. aqua2x=: 3 : '({@>|.TMP),L:0;("1),.(# TMP)#<{@>TMP=:i: y' NB. aqua2 2 3 aqua2=: 3 : '|.{(i:{.y);i:{:y' NB. real for log // x ,y, xy// not include 0 aqua2_x_real=: 3 : '|.{(i:{.y);i.{:y' aqua2_y_real=: 3 : '|.{(i.{.y);i:{:y' aqua2_xy_real=:3 : '|.{(i.{.y);i.{:y' NB. aqua2xr=: 3 : '({@>|.TMP),L:0;("1),.(# TMP)#<{@>TMP=:i. y' NB. -----class plot--------------------- plot_aqua=: 1 : 0 TMP=:;("1) u L:0 aqua2 y 'SR SC'=: <. -: $ TMP pd 'reset' pd (i:SR);(i:SC);TMP pd 'show' ) plot_aqua2_adv=: 1 : 0 NB.Usage: fun_04 test_adv 0 30 NB. M is 0 1 2 3/ aqua2 and real x y xy 'M Y0'=: y select. M case. 0 do. TMP=:;("1) u (L:0) aqua2 Y0 case. 1 do. TMP=: ;("1) u (L:0) aqua2_x_real Y0 case. 2 do. TMP=: ;("1) u (L:0) aqua2_y_real Y0 case. 3 do. TMP=: ;("1) u (L:0) aqua2_xy_real Y0 end. 'SR SC'=: <. -: $ TMP pd 'reset' pd (i:SR);(i:SC);TMP pd 'show' ) NB. -----class function--------------------- fun_01=: 3 : '-/ y ^2' NB. x^2 - y^2 fun_02=: 3 : '+/ y ^2' NB. x^2 + y^2 NB. ----------------------------------- fun_03=: 3 : '+/(2*{.y),(^.{. y), -({:y)^2' NB. f(x,y)=2x-(y^2)+logx NB. complex number fun_04=: 3 : '+/6 3 _5 *(*/ y),^&2 y ' NB. f=3x^2 + 6xy-5y^2 NB. fun_04 plot_aqua 30 fun_05=: 3 : '1 - {: y' NB. 1-y fun_06=: 3 : '-{:y' NB. -y fun_07=: 3 : '(1x1 * {. y)-{:y'