3次元プロット

実験屋らしい使い方も書いておこう.
これは原子によるレーザーの散乱レート
\gamma_p = \frac \gamma 2 \frac s {1+s+(2\delta/\gamma)^2}
をレーザー強度および離調の2変数関数としてプロットしたものである.飽和パラメタsはレーザー強度を飽和強度Isat = 1.67 mW/cm^2(Rb D2線)で正規化したもので,0から10まで動かしている.また離調δは0から20 MHzまで動かしている(Rb D2線の自然幅γ=6.07 MHz).

load 'plot numeric'
gamma =. 6.07
n =. 20
plot (steps 0 10,n) (([ ; ]) , <@((o. 2*gamma)&*)@(([ % >:@(+ *:@((gamma & (%~))@+:)))"0/)) (steps 0 20,n)