// sum 3-16
clc;
clear;
D=22;
d=20;
r=1;
K=2.2;
sigmax=130;
sigmax=sigmax/K;
Z=%pi*d^3/32;
M=sigmax*Z*10^-3;

  // printing data in scilab o/p window
  printf("M is %0.3f Nm    ",M);