In[1]:= fiBL[x_,z_] = (1-x/h)(1-z/h) x z Out[1]= (1 - -) (1 - -) h h In[2]:= fiBR[x_,z_] = (x/h)(1-z/h) z x (1 - -) h Out[2]= --------- h In[3]:= fiTR[x_,z_] = (x/h)(z/h) x z Out[3]= --- 2 h In[4]:= fiTL[x_,z_] = (1-x/h)(z/h) x (1 - -) z h Out[4]= --------- h In[5]:= psiL[x_,z_] = 1 - x/h x Out[5]= 1 - - h In[8]:= psiB[x_,z_] = 1/2(-1-z/h) z -1 - - h Out[8]= ------ 2 In[9]:= psiT[x_,z_] = z/h z Out[9]= - h In[11]:= Integrate[fiBL[x,z] fiBL[x,z], {x, 0, h}, {z, 0 , h}] 2 h Out[11]= -- 9 In[12]:= Integrate[fiBR[x,z] fiBL[x,z], {x, 0, h}, {z, 0 , h}] 2 h Out[12]= -- 18 In[13]:= Integrate[fiTR[x,z] fiBL[x,z], {x, 0, h}, {z, 0 , h}] 2 h Out[13]= -- 36 In[14]:= Integrate[fiTL[x,z] fiBL[x,z], {x, 0, h}, {z, 0 , h}] 2 h Out[14]= -- 18 In[15]:= Integrate[psiL[x,z] fiBL[x,z], {x, 0, h}, {z, 0 , h}] 2 h Out[15]= -- 6 In[16]:= Integrate[psiR[x,z] fiBL[x,z], {x, 0, h}, {z, 0 , h}] 2 h Out[16]= -- 12 In[17]:= Integrate[fiBR[x,z] fiBR[x,z], {x, 0, h}, {z, 0 , h}] 2 h Out[17]= -- 9 In[18]:= Integrate[fiTR[x,z] fiBR[x,z], {x, 0, h}, {z, 0 , h}] 2 h Out[18]= -- 18 In[19]:= Integrate[fiTL[x,z] fiBR[x,z], {x, 0, h}, {z, 0 , h}] 2 h Out[19]= -- 36 In[20]:= Integrate[psiL[x,z] fiBR[x,z], {x, 0, h}, {z, 0 , h}] 2 h Out[20]= -- 12 In[21]:= Integrate[psiR[x,z] fiBR[x,z], {x, 0, h}, {z, 0 , h}] 2 h Out[21]= -- 6 In[22]:= Integrate[fiTR[x,z] fiTR[x,z], {x, 0, h}, {z, 0 , h}] 2 h Out[22]= -- 9 In[23]:= Integrate[fiTL[x,z] fiTR[x,z], {x, 0, h}, {z, 0 , h}] 2 h Out[23]= -- 18 In[24]:= Integrate[psiL[x,z] fiTR[x,z], {x, 0, h}, {z, 0 , h}] 2 h Out[24]= -- 12 In[25]:= Integrate[psiR[x,z] fiTR[x,z], {x, 0, h}, {z, 0 , h}] 2 h Out[25]= -- 6 In[26]:= Integrate[fiTL[x,z] fiTL[x,z], {x, 0, h}, {z, 0 , h}] 2 h Out[26]= -- 9 In[27]:= Integrate[psiL[x,z] fiTL[x,z], {x, 0, h}, {z, 0 , h}] 2 h Out[27]= -- 6 In[28]:= Integrate[psiR[x,z] fiTL[x,z], {x, 0, h}, {z, 0 , h}] 2 h Out[28]= -- 12 *** for a_55 **************** In[29]:= Integrate[psiB[x,z] fiBL[x,z], {x, 0, h}, {z, 0 , h}] 2 -h Out[29]= --- 6 In[30]:= Integrate[psiB[x,z] psiB[x,z], {x, 0, h}, {z, 0 , h}] 2 7 h Out[30]= ---- 12 **** for a_56 ************* In[31]:= Integrate[psiB[x,z] fiBR[x,z], {x, 0, h}, {z, 0 , h}] 2 -h Out[31]= --- 6 **** for a_57 ************* In[32]:= Integrate[psiB[x,z] fiTR[x,z], {x, 0, h}, {z, 0 , h}] 2 -5 h Out[32]= ----- 24 **** for a_58 ************* In[33]:= Integrate[psiB[x,z] fiTL[x,z], {x, 0, h}, {z, 0 , h}] 2 -5 h Out[33]= ----- 24 ***** for a_{5,11} *********** In[35]:= Integrate[psiT[x,z] fiBL[x,z], {x, 0, h}, {z, 0 , h}] 2 h Out[35]= -- 12 In[36]:= Integrate[psiT[x,z] psiB[x,z], {x, 0, h}, {z, 0 , h}] 2 -5 h Out[36]= ----- 12 ***** for a_{6,11} *********** In[37]:= Integrate[psiT[x,z] fiBR[x,z], {x, 0, h}, {z, 0 , h}] 2 h Out[37]= -- 12 ***** for a_{7,11} *********** In[38]:= Integrate[psiT[x,z] fiTR[x,z], {x, 0, h}, {z, 0 , h}] 2 h Out[38]= -- 6 ***** for a_{8,11} *********** In[39]:= Integrate[psiT[x,z] fiTL[x,z], {x, 0, h}, {z, 0 , h}] 2 h Out[39]= -- 6 ***** for a_{99} *********** In[40]:= Integrate[psiL[x,z] psiL[x,z], {x, 0, h}, {z, 0 , h}] 2 h Out[40]= -- 3 ***** for a_{9,10} *********** In[42]:= Integrate[psiR[x,z] psiL[x,z], {x, 0, h}, {z, 0 , h}] 2 h Out[42]= -- 6 **** for a_{10,10} *********** In[43]:= Integrate[psiR[x,z] psiR[x,z], {x, 0, h}, {z, 0 , h}] 2 h Out[43]= -- 3 *** for a_{11,11} *********** In[44]:= Integrate[psiT[x,z] psiT[x,z], {x, 0, h}, {z, 0 , h}] 2 h Out[44]= -- 3 In[45]:= dxfiBL[x_,z_] = D[fiBL[x,z],x] z 1 - - h Out[45]= -(-----) h In[46]:= dzfiBL[x_,z_] = D[fiBL[x,z],z] x 1 - - h Out[46]= -(-----) h In[47]:= dxfiBR[x_,z_] = D[fiBR[x,z],x] z 1 - - h Out[47]= ----- h In[48]:= dzfiBR[x_,z_] = D[fiBR[x,z],z] x Out[48]= -(--) 2 h In[49]:= dxfiTR[x_,z_] = D[fiTR[x,z],x] z Out[49]= -- 2 h In[50]:= dzfiTR[x_,z_] = D[fiTR[x,z],z] x Out[50]= -- 2 h In[51]:= dxfiTL[x_,z_] = D[fiTL[x,z],x] z Out[51]= -(--) 2 h In[52]:= dzfiTL[x_,z_] = D[fiTL[x,z],z] x 1 - - h Out[52]= ----- h In[53]:= dxpsiL[x_,z_] = D[psiL[x,z],x] 1 Out[53]= -(-) h In[54]:= dxpsiR[x_,z_] = D[psiR[x,z],x] 1 Out[54]= - h In[55]:= dzpsiB[x_,z_] = D[psiB[x,z],z] -1 Out[55]= --- 2 h In[56]:= dzpsiT[x_,z_] = D[psiT[x,z],z] 1 Out[56]= - h *****************for a_11 *************** In[57]:= Integrate[dxfiBL[x,z] dxfiBL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[57]= - 3 In[58]:= Integrate[dzfiBL[x,z] dzfiBL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[58]= - 3 *****************for a_12 *************** In[59]:= Integrate[dxfiBR[x,z] dxfiBL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[59]= -(-) 3 In[60]:= Integrate[dzfiBR[x,z] dzfiBL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[60]= - 6 *****************for a_13 *************** In[61]:= Integrate[dxfiTR[x,z] dxfiBL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[61]= -(-) 6 In[62]:= Integrate[dzfiTR[x,z] dzfiBL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[62]= -(-) 6 *****************for a_14 *************** In[63]:= Integrate[dxfiTL[x,z] dxfiBL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[63]= - 6 In[64]:= Integrate[dzfiTL[x,z] dzfiBL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[64]= -(-) 3 *****************for a_15 *************** In[65]:= Integrate[dzfiBL[x,z] dxfiBL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[65]= - 4 In[66]:= Integrate[dxfiBL[x,z] dzfiBL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[66]= - 4 In[67]:= Integrate[dzpsiB[x,z] dxfiBL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[67]= - 4 *****************for a_16 *************** In[68]:= Integrate[dzfiBR[x,z] dxfiBL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[68]= - 4 In[69]:= Integrate[dxfiBR[x,z] dzfiBL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[69]= -(-) 4 *****************for a_17 *************** In[70]:= Integrate[dzfiTR[x,z] dxfiBL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[70]= -(-) 4 In[71]:= Integrate[dxfiTR[x,z] dzfiBL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[71]= -(-) 4 *****************for a_18 *************** In[72]:= Integrate[dzfiTL[x,z] dxfiBL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[72]= -(-) 4 In[73]:= Integrate[dxfiTL[x,z] dzfiBL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[73]= - 4 ************************* for a_19 ******************* In[74]:= Integrate[dxpsiL[x,z] dxfiBL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[74]= - 2 ************************* for a_1,10 ******************* In[75]:= Integrate[dxpsiR[x,z] dxfiBL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[75]= -(-) 2 ************************* for a_1,11 ******************* In[76]:= Integrate[dzpsiT[x,z] dxfiBL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[76]= -(-) 2 ************************ for a_22 ******************* In[77]:= Integrate[dxfiBR[x,z] dxfiBR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[77]= - 3 In[78]:= Integrate[dzfiBR[x,z] dzfiBR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[78]= - 3 ************************ For a_23 **************************** In[79]:= Integrate[dxfiTR[x,z] dxfiBR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[79]= - 6 In[80]:= Integrate[dzfiTR[x,z] dzfiBR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[80]= -(-) 3 ************************ For a_24 **************************** In[81]:= Integrate[dxfiTL[x,z] dxfiBR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[81]= -(-) 6 In[82]:= Integrate[dzfiTL[x,z] dzfiBR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[82]= -(-) 6 ************************ For a_25 **************************** In[83]:= Integrate[dzfiBL[x,z] dxfiBR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[83]= -(-) 4 In[84]:= Integrate[dxfiBL[x,z] dzfiBR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[84]= - 4 In[85]:= Integrate[dzpsiB[x,z] dxfiBR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[85]= -(-) 4 ************************ For a_26 **************************** In[86]:= Integrate[dzfiBR[x,z] dxfiBR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[86]= -(-) 4 In[87]:= Integrate[dxfiBR[x,z] dzfiBR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[87]= -(-) 4 ************************ For a_27 **************************** In[88]:= Integrate[dzfiTR[x,z] dxfiBR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[88]= - 4 In[89]:= Integrate[dxfiTR[x,z] dzfiBR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[89]= -(-) 4 *********************** For a_28 **************************** In[90]:= Integrate[dzfiTL[x,z] dxfiBR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[90]= - 4 In[91]:= Integrate[dxfiTL[x,z] dzfiBR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[91]= - 4 ********************** For a_2,9 **************************** In[94]:= Integrate[dxpsiL[x,z] dxfiBR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[94]= -(-) 2 *********************** For a_2,10 **************************** In[92]:= Integrate[dxpsiR[x,z] dxfiBR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[92]= - 2 ********************** For a_2,11 **************************** In[95]:= Integrate[dzpsiT[x,z] dxfiBR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[95]= - 2 ************************ For a_33 **************************** In[96]:= Integrate[dxfiTR[x,z] dxfiTR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[96]= - 3 In[97]:= Integrate[dzfiTR[x,z] dzfiTR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[97]= - 3 ************************ For a_34 **************************** In[98]:= Integrate[dxfiTL[x,z] dxfiTR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[98]= -(-) 3 In[99]:= Integrate[dzfiTL[x,z] dzfiTR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[99]= - 6 ************************ For a_35 **************************** In[100]:= Integrate[dzfiBL[x,z] dxfiTR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[100]= -(-) 4 In[101]:= Integrate[dxfiBL[x,z] dzfiTR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[101]= -(-) 4 In[102]:= Integrate[dzpsiB[x,z] dxfiTR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[102]= -(-) 4 ************************ For a_36 **************************** In[103]:= Integrate[dzfiBR[x,z] dxfiTR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[103]= -(-) 4 In[104]:= Integrate[dxfiBR[x,z] dzfiTR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[104]= - 4 ************************ For a_37 **************************** In[105]:= Integrate[dzfiTR[x,z] dxfiTR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[105]= - 4 In[106]:= Integrate[dxfiTR[x,z] dzfiTR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[106]= - 4 ************************ For a_38 **************************** In[107]:= Integrate[dzfiTL[x,z] dxfiTR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[107]= - 4 In[108]:= Integrate[dxfiTL[x,z] dzfiTR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[108]= -(-) 4 ************************ For a_39 **************************** In[109]:= Integrate[dxpsiL[x,z] dxfiTR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[109]= -(-) 2 ************************ For a_3,10 **************************** In[110]:= Integrate[dxpsiR[x,z] dxfiTR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[110]= - 2 ************************ For a_3,11 **************************** In[111]:= Integrate[dzpsiT[x,z] dxfiTR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[111]= - 2 ************************ For a_44 **************************** In[112]:= Integrate[dxfiTL[x,z] dxfiTL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[112]= - 3 In[113]:= Integrate[dzfiTL[x,z] dzfiTL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[113]= - 3 ************************ For a_45 **************************** In[114]:= Integrate[dzfiBL[x,z] dxfiTL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[114]= - 4 In[115]:= Integrate[dxfiBL[x,z] dzfiTL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[115]= -(-) 4 In[116]:= Integrate[dzpsiB[x,z] dxfiTL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[116]= - 4 ************************ For a_46 **************************** In[117]:= Integrate[dzfiBR[x,z] dxfiTL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[117]= - 4 In[118]:= Integrate[dxfiBR[x,z] dzfiTL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[118]= - 4 *********************** For a_47 **************************** In[119]:= Integrate[dzfiTR[x,z] dxfiTL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[119]= -(-) 4 In[120]:= Integrate[dxfiTR[x,z] dzfiTL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[120]= - 4 *********************** For a_48 **************************** In[121]:= Integrate[dzfiTL[x,z] dxfiTL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[121]= -(-) 4 In[122]:= Integrate[dxfiTL[x,z] dzfiTL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[122]= -(-) 4 *********************** For a_49 **************************** In[123]:= Integrate[dxpsiL[x,z] dxfiTL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[123]= - 2 *********************** For a_4,10 **************************** In[124]:= Integrate[dxpsiR[x,z] dxfiTL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[124]= -(-) 2 *********************** For a_4,11 **************************** In[125]:= Integrate[dzpsiT[x,z] dxfiTL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[125]= -(-) 2 *********************** For a_55 **************************** In[126]:= Integrate[dzpsiB[x,z] dzfiBL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[126]= - 4 In[127]:= Integrate[dzpsiB[x,z] dzpsiB[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[127]= - 4 *********************** For a_56 **************************** In[128]:= Integrate[dzpsiB[x,z] dzfiBR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[128]= - 4 *********************** For a_57 **************************** In[129]:= Integrate[dzpsiB[x,z] dzfiTR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[129]= -(-) 4 *********************** For a_58 **************************** In[130]:= Integrate[dzpsiB[x,z] dzfiTL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[130]= -(-) 4 *********************** For a_59 **************************** In[131]:= Integrate[dxpsiL[x,z] dzfiBL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[131]= - 2 In[132]:= Integrate[dxpsiL[x,z] dzpsiB[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[132]= - 2 *********************** For a_5,10 **************************** In[133]:= Integrate[dxpsiR[x,z] dzfiBL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[133]= -(-) 2 In[134]:= Integrate[dxpsiR[x,z] dzpsiB[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[134]= -(-) 2 *********************** For a_5,11 **************************** In[135]:= Integrate[dzpsiT[x,z] dzfiBL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[135]= -(-) 2 In[136]:= Integrate[dzpsiT[x,z] dzpsiB[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[136]= -(-) 2 *********************** For a_66 **************************** In[137]:= Integrate[dzpsiB[x,z] dzfiBR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[137]= - 4 ********************** For a_6,9 **************************** In[138]:= Integrate[dxpsiL[x,z] dzfiBR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[138]= - 2 In[139]:= Integrate[dxpsiL[x,z] dzpsiB[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[139]= - 2 ********************** For a_6,10 **************************** In[140]:= Integrate[dxpsiR[x,z] dzfiBR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[140]= -(-) 2 In[141]:= Integrate[dxpsiL[x,z] dzpsiB[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[141]= - 2 ********************** For a_6,11 **************************** In[142]:= Integrate[dzpsiT[x,z] dzfiBR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[142]= -(-) 2 In[143]:= Integrate[dzpsiT[x,z] dzpsiB[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[143]= -(-) 2 ********************** For a_7,9 **************************** In[144]:= Integrate[dxpsiL[x,z] dzfiTR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[144]= -(-) 2 ********************** For a_7,10 **************************** In[145]:= Integrate[dxpsiR[x,z] dzfiTR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[145]= - 2 ********************* For a_7,11 **************************** In[146]:= Integrate[dzpsiT[x,z] dzfiTR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[146]= - 2 ********************** For a_8,9 **************************** In[147]:= Integrate[dxpsiL[x,z] dzfiTL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[147]= -(-) 2 ********************** For a_8,10 **************************** In[148]:= Integrate[dxpsiR[x,z] dzfiTL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[148]= - 2 ********************** For a_8,11 **************************** In[149]:= Integrate[dzpsiT[x,z] dzfiTL[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[149]= - 2 ********************** For a_99 **************************** In[150]:= Integrate[dxpsiL[x,z] dxpsiL[x,z] , {x, 0, h}, {z, 0 , h}] Out[150]= 1 ********************** For a_9,10 **************************** In[151]:= Integrate[dxpsiR[x,z] dxpsiL[x,z] , {x, 0, h}, {z, 0 , h}] Out[151]= -1 ********************** For a_9,11 **************************** In[153]:= Integrate[dzpsiT[x,z] dxpsiL[x,z] , {x, 0, h}, {z, 0 , h}] Out[153]= -1 ********************** For a_10,10 **************************** In[154]:= Integrate[dxpsiR[x,z] dxpsiR[x,z] , {x, 0, h}, {z, 0 , h}] Out[154]= 1 ********************** For a_10,11 **************************** In[155]:= Integrate[dzpsiB[x,z] dxpsiR[x,z] , {x, 0, h}, {z, 0 , h}] 1 Out[155]= -(-) 2 ********************** For a_11,11 **************************** In[156]:= Integrate[dzpsiT[x,z] dzpsiT[x,z] , {x, 0, h}, {z, 0 , h}] Out[156]= 1