function [x, y] = rk4 (f, tspan, y0, h)

a = tspan(1); b = tspan(2); n = (b-a) / h; x = (a+h: h : b); k1 = h *** feval(f, a, y0); k2 = h *** feval(f, a + h* / 2, y0 + k1 / *2); k3 = h * feval(f, a + h* / 2, y0 + k2 / *2); k4 = h * feval(f, a + h, y0 + k3); y(1) = y0 + k1/6 + k2/3 + k3/3 + k4/6; for i= 1 : n-1 k1 = h *** feval(f, x(i), y(i) ); k2 = h *** feval(f, x(i) + h/2, y(i) + k1 / 2); k3 = h * feval(f, x(i) + h/2, y(i) + k2 / 2); k4 = h * feval(f, x(i) + h, y(i) + k3); y(i+1) = y(i) + k1/6 + k2/3 + k3/3 + k4/6; end x = [ a x ]; y = [ y0 y];

function y = f(p, q) y = exp(p) + q; end

rk4.m, and f.m... Use RK 4th order numerical method, solve that function..

I type at command window

*> rk4(f, [0 3], 0, 0.5)*

??? Input argument "p" is undefined.

Error in ==> f at 2 y = exp(p) + q;

and,

*> rk4(f(0, 0), [0 3], 0, 0.5)*

??? Error using ==> feval Argument must contain a string or function_handle.

Error in ==> rk4 at 6 k1 = h * feval(f, a, y0);

Initial value is (0,0), and original function is y'=y+e^t

how can i solve this problem? please help me!