683/CH3/EX3.3/MS_3.sce


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17

// sum 3-3
clc;
clear;
m=25;
v=3;
E=210*10^3;
KE=0.5*m*v^2;
d=30;
L=2000;
A=%pi*d^2/4;
U=A*L/(2*E);
del=4*10^-5*A;
W=A*del;
sigi=sqrt(KE*10^3/(W+U));

  // printing data in scilab o/p window
  printf("del is %f N/mm^2     ",sigi);