matlab - Matlab 微分

Question

前方差分法を使用して導関数を解くには、matlab で for ループを記述する必要があります。導出する関数は、0.25 のステップを使用して 0 から 10 までの 10+15x+20x^2 です。使ってみました

h=.25; 
x=[0:h:10]; 
y = 10+15*x+20*x.^2; 
y(1) = 45; size(x)  
for i=2:47,     
       y(i) = y(i-1) + h*(15+40*x); 
end

Answer 1

私は、まず、このようにします。

h=.25; 
x=[0:h:10]; 
y = 10+15*x+20*x.^2; 
diff(y)./diff(x)

または、代わりに、

syms x;
y = 20.*x.^2 + 15.*x + 10;
dy = diff(y,1);
h=.25; 
xx=[0:h:10];
res = subs(dy,xx);