結構めんどくさい うまく動いてない途中で三重対角行列を係数行列とする連立１次方程式の解を求める必要があったが、TDMA（三重対角行列アルゴリズム）という解法がある。 function spline_interpolation clear all; close all; format compact; div = 100; …

N=3ですでにサインカーブがうまく補間されている function lagrange_interpolation clear all; close all; div = 100; x = linspace(0,pi,100); y = sin(x); plot(x,y); hold on; a = x(1); b = x(end); fx = (x-b)/(a-b)*sin(a) + (x-a)/(b-a)*sin(b); plot…

http://www.takuichi.net/hobby/edu/em/mom/numerical/index-j.htm を計算する リーマン積分（区分求積法） function riemann clear all; close all; divlist = [10e2,10e3,10e4,10e5,10e6,10e7]; errorlist = []; for cycle=1:6 div = divlist(cycle); x = …

スクリーンキャプチャした （１）isimを走らせるときはテストベンチを描いてもいいし、テストベンチがない場合でも、GUI/CUIでClockやConstantを与えることができる→ユニットテストが手軽にできて便利→isimがprintfとかデバッガ的な位置づけになる （２）Clo…

スクリーンキャプチャしといた （１）トップレベルに対してのみ回路図は生成できる （２）ダブルクリックでドリルダウン（サブモジュールの中に入る）できる （３）Ctrl-PでPDF出力できる

function conway2d clear all; close all; CAsize = 100; [X,Y] = meshgrid(1,CAsize+1); figure(); hold on; plot(X,Y,'k'); plot(Y,X,'k'); I = zeros(CAsize+1); surface(I); colormap(gray); axis off; I = rand(CAsize+1)>0.85; I(1,:) = 0; I(CAsize,:…