%Specify Cylindrical Coordinate Components
rhonot = 0.5;
delrho = 0.5;
delpsi = pi/8;

%Boundary Conditions
bc = [1 1 65;2 1 65;3 1 65;4 1 65;5 1 65;6 1 65;7 1 65;8 1 65;9 1 65;10 1 65;11 1 65;12 1 65;13 1 65;14 1 65;...
    2 1 0;2 2 0;2 8 0;2 9 0;2 10 0;2 15 0;3 1 0;3 2 0;3 3 0;3 7 0;3 8 0;3 9 0;3 10 0;3 11 0;...
    3 14 0;3 15 0;4 1 0;4 2 0;4 3 0;4 5 0;4 7 0;4 8 0;4 9 0;4 10 0;4 11 0;4 14 0;4 15 0;
    5 1 0;5 2 0;5 3 0;5 4 0;5 6 0;5 7 0;5 8 0;5 9 0;5 10 0;5 11 0;5 12 0;5 13 0;5 14 0;5 15 0];