seinzumtode https://blog.hatena.ne.jp/seinzumtode/ how to code something https://seinzumtode.hatenadiary.jp/ を考える。 clear; close all; Ns=[8,10,15,20,30]; x0=0; xN=20; figure(); hold on; for idx=1:length(Ns) N = Ns(idx); x=linspace(x0,xN,N); y=zeros(1,N); y(1)=0.1; h=(xN-x0)/N; f = @(y) y-y^2; for k=2:N y(k) = y(k-1) + h*(f(y(k-1))); end plot(1:N,y,'DisplayName',sprintf('N=%d',N)); xlim([0,10]); end [x2,y2]=… 190 <iframe src="https://hatenablog-parts.com/embed?url=https%3A%2F%2Fseinzumtode.hatenadiary.jp%2Fentry%2F2020%2F05%2F18%2F111737" title="微分方程式の初期値問題（オイラー法での数値振動） - 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://chart.apis.google.com/chart?cht=tx&chl=%20%5Cdfrac%7Bdy%7D%7Bdx%7D%3Dy-y%5E3 Hatena Blog https://hatena.blog 2020-05-18 11:17:37 rich https://seinzumtode.hatenadiary.jp/entry/2020/05/18/111737 1.0 100%