Hi Bruno, Thanks for the note. As you mentioned the definition of u and v was missing. I would be very curious if anyone has an elegant way of eliminating the interior surfaces that only the outer facetet surface remains.
Thanks a lot, Peter
Here is the complete code:
clear all clc close all a=8; b=8; c=50; d=200; nf_uv_grid = 101; %alpha=0/180*pi; %beta=0/180*pi; gamma=0/180*pi;