rk4 - k3 = h * func(t(i,1)+.5*h,w(:,i)+.5*k2); k4 = h *

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
function [t,w] = rk4(func,a,b,h,initval) % % Runge Kutta 4th order % n = (b-a)/h; t(1,1) = a; w = initval; w for i = 1:n k1 = h * func(t(i,1),w(:,i)); k2 = h * func(t(i,1)+.5*h,w(:,i)+.5*k1);
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: k3 = h * func(t(i,1)+.5*h,w(:,i)+.5*k2); k4 = h * func(t(i,1)+h,w(:,i)+k3); w(:,i+1) = w(:,i) + (k1 + 2*k2 + 2*k3 + k4)/6; t(i+1,1) = t(i,1) + h; end;...
View Full Document

Ask a homework question - tutors are online