%Keler box Blasius

x1=0; 
x2=10;
h=0.01;
N=(x2-x1)/h+1;

a=ones(N,1);
b=zeros(N,1);
c=zeros(N,1);
for i=1:N
   b(i)=0+i/N;
end;

err=1;
nr_iter=0;
while abs(err)>=1e-6
[A,B,C,r]=init_matrix(a,b,c,N,h);
delta=elim_bloc(A,B,C,r,N);
err=max(max(delta));
[a,b,c]=reeval(delta,a,b,c,N);
nr_iter=nr_iter+1;
end;
plot(x1:h:x2,b);
fprintf('\n iteratii: %d   f"(0)=%f \n',nr_iter,c(1))






 



