integration_routine

integration_routine - k2=eqns(n, x0+(3/4)*h, y+(3/4)*k1*h);

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

View Full Document Right Arrow Icon
function [y] = integration_routine(n, y, x0, xf, dx) while xf>x0 if (xf-x0) < dx; dx=(xf-x0); else end h=dx; k1=eqns(n, x0, y);
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: k2=eqns(n, x0+(3/4)*h, y+(3/4)*k1*h); y=y+((1/3)*k1+(2/3)*k2)*h; x0=x0+dx; end...
View Full Document

This note was uploaded on 07/12/2011 for the course ECON 011 taught by Professor Bob during the Spring '11 term at Tufts.

Ask a homework question - tutors are online