(* ::Package:: *)

Print["initilisation.m is run"];


L=10;tol=0.0001;Vb1=-1.0;Vb2=1.0;
iini=0;ilast=L;
jini=0;jlast=L;
xini=0;xlast=1.0;yini=0;ylast=1.0;
Delta=(xlast-xini)/L;
For[i=iini,i<=ilast,i++,
For[j=jini,j<=jlast,j++,
V[i,j]=0.0;
];(*end for i*)
];(*end for i*)

For[i=iini,i<=ilast,i++,
x[i]=xini+Delta*i;
];(*end for i*)
For[j=jini,j<=jlast,j++,
y[j]=yini+Delta*j;
];(*end for i*)

j=jini;
For[i=iini,i<=ilast,i++,
V[i,j]=(2/(xlast-xini))*x[i]-1;
];(*end for i*)

j=jlast;
For[i=iini,i<=ilast,i++,
V[i,j]=(2/(xlast-xini))*x[i]-1;
];(*end for i*)

txyV=Table[{x[i],y[j],V[i,j]},{j,jini,jlast},{i,iini,ilast}];
txyV=Flatten[txyV,1];
