1 files changed, 50 insertions, 0 deletions
diff --git a/1445/CH10/EX10.16/ch10_ex_16.sce b/1445/CH10/EX10.16/ch10_ex_16.sce
new file mode 100644
index 000000000..f2fd39150
--- /dev/null
+++ b/1445/CH10/EX10.16/ch10_ex_16.sce
@@ -0,0 +1,50 @@
+//CHAPTER 10- THREE-PHASE INDUCTION MACHINES
+//Example 16
+
+disp("CHAPTER 10");
+disp("EXAMPLE 16");
+
+//VARIABLE INITIALIZATION
+ratio1=1.5;                    //ratio of T_est and T_efl
+ratio2=2.5;                    //ratio of T_em and T_efl
+
+//SOLUTION
+s=1;
+
+//solution (a)
+//directly solving the quadratic equation
+a=1;
+b=-3.333;
+c=1;
+D=(b)^2-(4*a*c);               //discriminant
+sm1=(-b+sqrt(D))/(2*a);
+sm2=(-b-sqrt(D))/(2*a);
+if(sm1<=0 & sm2<=0) then
+disp("The value of the slip at maximum torque is not valid");
+else if(sm1>0 & sm1<1) 
+disp(sprintf("The slip at maximum torque is %f",sm1));
+else if(sm2>0 & sm2<1) 
+disp(sprintf("The slip at maximum torque is %f",sm2));
+end; 
+
+//solution (b)
+//directly solving the quadratic equation
+a=1;
+b=-1.665;
+c=0.111;
+D=(b)^2-(4*a*c);          
+ans1=(-b+sqrt(D))/(2*a); 
+ans2=(-b-sqrt(D))/(2*a);
+if(ans1>0 & ans1<1)
+disp(sprintf("The full load slip is %f",ans1));
+sfl=ans1;
+else if(ans2>0 & ans2<1)
+disp(sprintf("The full load slip is %f",ans2));
+sfl=ans2;
+end;
+
+//solution (c)
+I=sqrt(ratio1/sfl);
+disp(sprintf("The rotor current at the starting in terms of full load current is %f A",I));
+
+//END