tspan = [0 20];y0 = [0 0.01];[z,y] = ode45(@odefcn, tspan, y0);plot(z,y(:,1),'-o',z,y(:,2),'-.')function dydt = odefcn(z,y)dydt1 = y(2);dydt2 = NS*z.*y(1);dydt=[dydt1;dydt2];end function M = NS(z)z = [2 3 5 7 10 15 20 ];r =[3.5 3.7 4 6 7.2 8 9];n=length(z);% Calculation of differentiation from the above datas
for i=1:n-1 M=zeros(); M(i)=(r(i+1)-r(i))./(z(i+1)-z(i)); end end
MATLAB: I am having trouble multiplying the main ode function with an external function NS which is supposed to be multiplied on the RHS. Thanks for the great help.
differentiation data pointsode45
Best Answer