Mercurial > repos > public > sbplib

comparison +sbp/ordinary12.m @ 257:5714fda0db85 operator_remake

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Bug fix 12th order 2nd der.
author Martin Almquist <martin.almquist@it.uu.se>
date Wed, 07 Sep 2016 16:46:20 +0200
parents df3cc9c5dffc
children
comparison
equal deleted inserted replaced
256:87504db14c6d 257:5714fda0db85
280 e_1(1)= 1; 280 e_1(1)= 1;
281 e_m = zeros(m,1); 281 e_m = zeros(m,1);
282 e_m(end)= 1; 282 e_m(end)= 1;
283 283
284 S_1 = (e_1'*D1)'; 284 S_1 = (e_1'*D1)';
285 S_m = (e_m'*D1); 285 S_m = (e_m'*D1)';
286 286
287 Q = H*D1-(-e_1*e_1' + e_m*e_m'); 287 Q = H*D1-(-e_1*e_1' + e_m*e_m');
288 M = -(H*D2-(-e_1*S_1' + e_m*S_m')); 288 M = -(H*D2-(-e_1*S_1' + e_m*S_m'));
289 end 289 end