Cでもそれなりに簡単に，数値計算系のプログラムを書く事はできますが
それを，表示するのがめんどくさいです．
csv吐いて，gnuplotに食わすのはつまらないのでMatlab互換のGNU Octaveを使ってみます．
リンク †
http://ayapin.film.s.dendai.ac.jp/~matuda/PlotUtils/Octave/octave.html
http://www.al.cs.kobe-u.ac.jp/~inamoto/unix-tools/useful/octave/index.html#sq
書きなぐり †
気に入ったら，丁寧に書くかも
- sin
octave:1> x=[0:0.1:2*pi]'; octave:2> sinx=[x,sin(x)]; octave:3> gplot sinx
行列は [1,2,3,4,5]みたいにつくる．[1:5]でも良い．
転置は[1;2;3;4;5]とか[1,2,3,4,5]'(シングルクウォートを後ろに付ける)gplot [x,sin(x)]
が出来ない辺りがショボイ．
gplotにはx0 y0 x1 y1 : :
を渡せば良いと思われる．
- circle
横に伸びるだけじゃつまらないrad = [0:0.1:2*pi]'; x = cos(rad); y = sin(rad); circle = [x,y]; gplot circle
- 関数
関数も作れますfunction plot_(args) x=[1:length(args)]'; line=[x,args']; gplot line endfunction
plot_ ([1,2,3])
こんな風に使えます．
- 行列の計算も楽チン
- スプライン補完
function ret=spline_(pnts) v0=(pnts(3,[1,2])-pnts(1,[1,2]))/2; v1=(pnts(4,[1,2])-pnts(2,[1,2]))/2; ft=[ 2*pnts(2,[1,2]) - 2*pnts(3,[1,2]) + v0 + v1; -3*pnts(2,[1,2]) + 3*pnts(3,[1,2]) - 2*v0 - v1; v0; pnts(2,[1,2])]; ii=1; for t= [0:0.05:1] ts=[t**3, t**2, t, 1]; tmp=ts*ft; ret(ii,1)=tmp(1); ret(ii,2)=tmp(2); ii++; endfor endfunction function ret=spline(pnts) sz=size(pnts); ret=[]; for ii=[1:sz(1)-3] ret= concat(ret', spline_(pnts([ii:ii+3],[1,2]))')'; endfor endfunction function plot_spline(pnts) pnts=concat(pnts(1,[1,2])', pnts')'; pnts=concat(pnts', pnts(size(pnts)(1), [1,2])')'; sp=spline(pnts); sz=size(sp); xs=sp([1:sz(1)],1); ys=sp([1:sz(1)],2); plot(xs,ys); endfunctionoctave:2> spline([1 1;3 1; 3 3; 1 3])
とかすると，スプライン補完します．
octave:3> plot_spline([1 1;3 1; 3 3; 1 3])
プロットもします．
ただ，C++でもこれ位なら演算子のオーバーロードだけで書けます．うーんoctaveの使い道がまだ分からない．
もっと重い処理をすると，Octaveの中の人の工夫が見えてくるかな
- スプライン補完
