Mercurial > repos > public > sbplib

comparison +sbp/D2.m @ 256:87504db14c6d operator_remake

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Bug fix in D2.m, added 12th order.
author Martin Almquist <martin.almquist@it.uu.se>
date Wed, 07 Sep 2016 16:41:31 +0200
parents df3cc9c5dffc
children 6009f2712d13
comparison
equal deleted inserted replaced
255:df3cc9c5dffc 256:87504db14c6d
24 [H, HI, D1, D2, e_1, e_m, M, Q, S_1, S_m] = sbp.ordinary8(m,h); 24 [H, HI, D1, D2, e_1, e_m, M, Q, S_1, S_m] = sbp.ordinary8(m,h);
25 obj.borrowing.M.S = 0.0015; 25 obj.borrowing.M.S = 0.0015;
26 elseif order == 10 26 elseif order == 10
27 [H, HI, D1, D2, e_1, e_m, M, Q, S_1, S_m] = sbp.ordinary10(m,h); 27 [H, HI, D1, D2, e_1, e_m, M, Q, S_1, S_m] = sbp.ordinary10(m,h);
28 obj.borrowing.M.S = 0.0351; 28 obj.borrowing.M.S = 0.0351;
29 elseif order == 12
30 [H, HI, D1, D2, e_1, e_m, M, Q, S_1, S_m] = sbp.ordinary12(m,h);
31 obj.borrowing.M.S = []; % Need to compute!
29 else 32 else
30 error('Invalid operator order %d.',order); 33 error('Invalid operator order %d.',order);
31 end 34 end
32 35
33 obj.h = h; 36 obj.h = h;