clc,clear
printf('Example 3.30\n\n')

V=400, I=50
I_a_g=250 ,I_a_m =300  //armature current for generator and motor
I_sh_g=2.5 ,I_sh_m =2.4 //field current for generator and motor
R_a=0.1 //armature resistance

arm_cu_loss_g = R_a*I_a_g^2   //armature copper loss for generator
arm_cu_loss_m = R_a*I_a_m^2   //armature copper loss for motor
power_drawn=V*I
IFW_losses = power_drawn - (arm_cu_loss_g + arm_cu_loss_m)  //Iron , friction and windage losses
IFW_losses_each= IFW_losses /2  // Iron , friction and windage losses for each machine

//for motor
field_cu_loss_m= V*I_sh_m //field copper loss for motor
total_loss_m= field_cu_loss_m + IFW_losses_each  +  arm_cu_loss_m
motor_input=V * I_a_m
motor_output= motor_input - total_loss_m
eta_m = 100*(motor_output/motor_input)    //motor efficiency
printf('Efficiency of motor is %.2f percent\n',eta_m)

//for generator
field_cu_loss_g= V*I_sh_g //field copper loss for generator
total_loss_g = field_cu_loss_g + arm_cu_loss_g + IFW_losses_each
generator_output=V*I_a_g
generator_input = generator_output + total_loss_g
eta_g = 100*(generator_output/generator_input)//generator efficiency
printf('Efficiency of generator is %.2f percent\n',eta_g)