Matlab - NO PUEDO CORRER UNA FUNCION DE MEODOS NUMERICOS..AYUDA

function [k, p, err, P]=puntofijo(p0,tol, max1)

P(1)=p0;

for k=2:max1

    P(k)=feval(@fun1,P(k-1));

    err=abs(P(k)-P(k-1));

    relerr=err/(abs(P(k))+eps);

    p=P(k);

    if err<tol ||  relerr<tol , break, end

end

if k==max1

    disp('se ha excedido el número maximo de iteraciones')

end

P=P';

function g=fun1(x)

g=sqrt(2*x+3);

>> p0=4; tol=0.01; max1=100;  [k, p, err, P]=puntofijo(p0,tol, max1)

k =

     5

p =

    3.0114

err =

    0.0229

P =

    4.0000

    3.3166

    3.1037

    3.0344

    3.0114

function [k, p, err, P]=puntofijooctave(g,p0,tol, max1)

  fid=fopen('funcion_octave.m','w+');

    fprintf(fid,'function y=funcion_octave(x)\n');

    fprintf(fid,'y=%s;',g);

  fclose(fid);

P(1)=p0;

for k=2:max1

    P(k)=feval(@funcion_octave,P(k-1));

    err=abs(P(k)-P(k-1));

    relerr=err/(abs(P(k))+eps);

    p=P(k);

    if err<tol ||  relerr<tol , break, end

end

if k==max1

    disp('se ha excedido el número maximo de iteraciones')

end

P=P';

>> g='sqrt(2*x+3)'; p0=4; tol=0.01; max1=100;  [k, p, err, P]=puntofijooctave(g,p0,tol, max1)

k = 5

p = 3.0114

err = 0.022945

P =

   4.0000

   3.3166

   3.1037

   3.0344

   3.0114