Comparing (37.16) with (37.2), we have
co = 0 |
(37.23) |
Comparing (37.17) and (37.18) with (37.2), we have
pl = cl -1 ql = 0 pr = cr qr = 0 |
(37.24) (37.25) (37.26) (37.27) |
Now the Matlab functions can be written as
Function pdefun
Function icfun
Function bcfun