seinzumtode https://blog.hatena.ne.jp/seinzumtode/ how to code something https://seinzumtode.hatenadiary.jp/ clear; close all; c = 1.0; dt = 0.05; dx = 0.1; xmax = 5; xmin = -1; xs = xmin:dx:xmax; q = zeros(1,length(xs)); q(xs<0) = 1; plot(xs,q,'b'); ns = [0,16,32,48]; nmax = max(ns); qs = zeros(length(ns),length(q)); qs(1,:) = q; qs_index = 2; for k=1:nmax qold = q; dq = zeros(1,length(xs)); c_a = abs(c);… 190 <iframe src="https://hatenablog-parts.com/embed?url=https%3A%2F%2Fseinzumtode.hatenadiary.jp%2Fentry%2F2025%2F08%2F02%2F232633" title="近似LDU分解による移流方程式の計算 - how to code something" class="embed-card embed-blogcard" scrolling="no" frameborder="0" style="display: block; width: 100%; height: 190px; max-width: 500px; margin: 10px 0px;"></iframe> https://cdn-ak.f.st-hatena.com/images/fotolife/s/seinzumtode/20250802/20250802232231.png Hatena Blog https://hatena.blog 2025-08-02 23:26:33 rich https://seinzumtode.hatenadiary.jp/entry/2025/08/02/232633 1.0 100%