How to set starting condition for Keller box method?

the initial condition are: f'(0)=1; f(0)=0; f'(inf)=0; p(inf)=0 G(inf)=f(inf) H(inf)=k1 where k1=0.2

so in matlab i set like this _______________________________________________________________________________________ % To generate the initial value for velocity, temperature and

    etau15 = 1 / eta(np,1);
    etau152 = 2 / eta(np,1);
    etanpq2 = eta(np,1) / 2;
    etanpq3 = eta(np,1) / 3;

    for j = 1:np
    deta(j,k) = deleta;    
      etab = eta(j,1) / eta(np,1);
      etab2 = etab * etab;
      etab3 = etab2 * etab;

        f(j,1) = eta(j,1) - eta(j,1)^2/eta(np,1) + eta(j,1)^3/3*(eta(np,1))^2;
        u(j,1) = (1-etab)^2;
        v(j,1) = -2*(eta(np,1)-eta(j,1))/(eta(np,1))^2;

        p(j,1) = (1-etab)^2;
        g(j,1) = - etau15;
        s(j,1) = -etau15;

    end
____________________________________________________________________________________

The answer is corect for value f,u,v but p,g and s are wrong.. anyone have suggestion of this???

[vijayakumar]

i am doing the same method but i was struggle how to set the Convective boundary condtions in my problem can you expalin how you define the code for mentioned boundary condition

Answer 1

[Shahid Rafiq]

what are ur odes? kindly give full code then it will be possible to check !

[Poornima T]

yes,kindly give us full details of what you have coded