帽子の形状を生成する曲面の式を #case(0) に登録します。
//**********************************************************************************
#macro Surface_func(Ut,Vt,Func_no)
#switch(Func_no)
#case(0)
#local Rp=(pow(.2*4*Ut,2)+pow(.2*4*Vt,2));
#local F=< 4*Ut, -4/sqrt(Rp)*sin(pi*sqrt(Rp)), 4*Vt >;
#break
#end F
#end
//**********************************************************************************