# Example 18.1      
 T1  =   85.0 N m 
 pmax   =   1400.0 kPa 
 ri   =   29.0 mm 
 ro   =   50.0 mm 
 f  =   0.06  
 N1 = T1 / (pi() * pmax * ri * f * (ro^2 - ri^2))  
 Nact   =   8.0  
 FOS = Nact / N1  
 ppmax = pmax / FOS ; kPa 
 Np = T1 / (pi() * ppmax * ri * f * (ro^2 - ri^2))   
 F1 = 2 *  T1 / (f * (ro + ri) * Np) ; N 

# Example 18.2      
BEGIN      
 F1  =   1000 N 
 H25 = 4 * F1 ; N 
 V25 = F1 ; N 
 H45 = 4 * F1 ; N 
 H34 = 4 * F1 ; N 
 H63 = 10.53 * F1 ; N 
 V63 = 2.11 * F1 ; N 
 H13 = 6.53 * F1 ; N 
 V13 = 2.11 * F1 ; N 
 H62 = 7.07 * F1 ; N 
 V62 = 1.41 * F1 ; N 
 H12 = 3.07 * F1 ; N 
 V12 = 2.41 * F1 ; N 
 H61 = 3.46 * F1 ; N 
 V61 = 0.7 * F1 ; N 
 R  =   250.0 mm 
 T0 = 3.52 * F1 * R ; N m 
 q  =   90.0 deg 
 w0  =   80.0 mm 
 A0 = w0 * 2 * R * sin(q/2)  ; mm^2 
 P3 = H63 / A0 ; Pa 
 pmax   =   0.4 N/mm^2 
 err = pmax - P3  ; Pa 
 Fb = sqrt(H61^2 + V61^2)  ; N 
END
solve(err, F1)

# Example 18.3
MKS(mm, deg, kPa, rpm)      
BEGIN      
 x1   =   150.0 mm 
 x2   =   200.0 mm 
 f = atan(x2 / x1) ; deg 
 q1 = f - pi() rad / 4  ; deg 
 q2 = q1 + pi() rad / 2  ; deg 
 d = sqrt(x1^2 + x2^2)  ; mm 
 sinq_max   =   1.0  
 pmax   =   445.6 kPa 
 b  =   50.0 mm 
 r  =   150.0 mm 
 d  =   250.0 mm 
 z = Number(2 * (q2 - q1)) + (sin(q1)^2 - sin(q2)^2)
 Mn = -pmax * b * r * d * z /  (4 * sinq_max)  ; N m 
 f  =   0.3  
 Mf = f * pmax * b * r * (r * (cos(q1) - cos(q2)) \
     + (d/4) * ( cos(q2)^2 - cos(q1)^2)) / sinq_max  ; N m 
 T0 = -r^2 * f * b * pmax * (cos(q1) - cos(q2)) / \
    sinq_max ; N m 
 n  =   300.0 rpm 
 v = r * n ; m/s 
 pmax_v  =   2.1 Mpa m / s 
 ppmax = pmax_v / v ; kPa 
 c  =   500.0 mm 
 F1 = (Mn + Mf) / c ; N 
 pmax2   =   391.7 kPa 
 err = -F1 *c + (pmax2 / pmax) * (Mn - Mf)  ; N m 
 pmax3   =   445.6 kPa 
 Tp =  (pmax2 / pmax3)  * T0 ; N m 
 Wdot = n * Tp  ; kW 
END   

# Example 18.4     
BEGIN     
 pmax   =   0.5 MPa
 r  =   250.0 mm
 b  =   80.0 mm
 P1 = pmax * r * b ; kN 
 f  =   0.2 
 phi  =   240.0 deg
 P2 = P1 / exp(f * phi)  ; kN 
 T0 = (P1 - P2) * r  ; N m 
 a  =   150.0 mm
 s0  =   35.0 mm
 c  =   700.0 mm
 F1 = (P2 * a - P1 * s0) / c ; N 
 sSL   =   64.9 mm
 FSL  = (P2 * a - P1 * sSL) / c ; N  
END      

# Homework Problems:

TBD