{"html":"<iframe src=\"https://hatenablog-parts.com/embed?url=https%3A%2F%2Fseinzumtode.hatenadiary.jp%2Fentry%2F2025%2F04%2F12%2F142817\" title=\"LU\u5206\u89e3\u3067\u306e\u9023\u7acb\u65b9\u7a0b\u5f0f\u306e\u89e3\u653e - 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>","height":"190","blog_url":"https://seinzumtode.hatenadiary.jp/","provider_name":"Hatena Blog","blog_title":"how to code something","description":"clear; close all;clc; A = [2 2 1; 3 6 2; 4 -2 5]; b = [4 2 8]'; [L, U] = lu_nopivot(A); y = forward_subs(L,b) [U2,y2] = forward_elim(U,y); %\u4e00\u5ea6\u5bfe\u89d2\u6210\u5206\u306e\u4fc2\u6570\u3092\u30bc\u30ed\u306b\u3059\u308b\u5fc5\u8981\u3042\u308a x = backward_subs(U2,y2); function [L,U] = lu_nopivot(A) n=size(A,1); L=eye(n); for k=1:n L(k+1:n,k) = A(k+1:n,k)/A(k,k); for l=k+1:n A(l,:)\u2026","title":"LU\u5206\u89e3\u3067\u306e\u9023\u7acb\u65b9\u7a0b\u5f0f\u306e\u89e3\u653e","provider_url":"https://hatena.blog","type":"rich","image_url":null,"url":"https://seinzumtode.hatenadiary.jp/entry/2025/04/12/142817","width":"100%","categories":[],"author_name":"seinzumtode","author_url":"https://blog.hatena.ne.jp/seinzumtode/","published":"2025-04-12 14:28:17","version":"1.0"}