clc; v=230; // source voltage ra=0.1; // resistance of armature ia=100; // armature current n=1600; // speed of dc shunt motor wl=300; // friction and windage losses lo=1200; // no load core loss lc=2500; // copper losses Ls=0.01; // stray losses as a fraction of output Ea=v-ia*ra; // counter EMF pe=Ea*ia; // electromagnetic power wo=wl+lo; // no load rotational losses po=pe-wo; // shaft power + stray load losses psh=po/(1+Ls); Tsh=(psh*60)/(2*%pi*n); printf('Shaft torque is %f Nm\n',Tsh); pi=pe+lc; // power input to motor nm=(psh/pi)*100; printf('Motor efficiency is %f percent',nm);