% Function file for passive cable branched. Branch 2
function [ VVb2,hb2,hb22]= volt_cabranpab2(VVp1,lp1,lb21,lambda1);
% function to create values of V for various values of x on child branch 1
% function to create values of h for various values of N 
% function to generate values of h2
% calling sequence
% [ Vb2,h,h2] = volt_cabranpab2(N,Vp1,lp, lb21, lambdab21,ra1,
% rab21)

% Define variables 
% N1 --number of discretization points wrt space
% Vb2 -- voltage across different points
% Vb2o-- voltage at initial point Iinj/ ra
% lb2 --length of dendrite in mm
% lambdab21 -- space constant in mm
% Xb2 --length nondimensionalised
% x -- points along the dendrite
% I1 -- injected current in milliamperes
% I11b21 -- non dimensionalised current 
% rmb21 -- memb res of branch 1 in ohm mm
% t1 -- time constant in msec
% Tb2 -- nondimensionalised T ie t1 / taub1
% kk -- index variable
% Record of revisions 
% Date                        Programmer           Description of changes
% Feb 24/ 07                 Asha G                   original code
% October 25,08              Asha G                 correction for N(1)
% declare global values 
%pack ( 'volt_cabranpab2')
global  Rm Ri  tau I N 
%Vb2 = zeros( 1,N);
VVb2= zeros(1,N);
x = linspace ( 1,1+(lb21/lp1),N);
X =  x/ lambda1;
Lb2 = lb21 / lambda1;

VVb2 = linspace(-70,-70,N); 

%VVb2(1) = VVp1(N);
%VVb2( N) = (4/3) * VVb2( N-1) - ( 1/3) * VVb2( N-2); 

hb2 = Lb2 / ( N-1);
hb22 = hb2 * hb2; 

end % volt_ cabranpab2