//Chapter 16,Example 6,page 560
//Determine the wave generated
clear
clc
// With refrence to table 16.1
C1 = 0.125*10^-6 // F
C2 = 1*10^-9 // F
R1 = 360 // ohm
R2 = 544 // ohm
V0 = 100 // kV
theta = sqrt(C1*C2*R1*R2)
neta = 1/(1+(1+R1/R2)*C2/C1)
alpha = R2*C1/(2*theta*neta)
printf("\n Theta = %f micro S",theta*10^6)
printf("\n Neta = %f",neta)
printf("\n Alpha = %f ",alpha)
// Coresponding to alpha the following can be deduced from Fig 16.12
T2 = 10.1*theta*10^6
T1 = T2/45
imp = T1/T2 // generated lighting impulse
// From equations 16.41 and 16.42
a1 = (alpha-sqrt(alpha^2-1))*10^-6/(theta) 
a2 = (alpha+sqrt(alpha^2-1))*10^-6/theta 
printf("\n T1 = %f microS", T1)
printf("\n T2 = %f microS", T2)
printf("\n Generated lighting impulse = %e wave", imp)
printf("\n alpha1 = %f microS", a1)
printf("\n alpha2 = %f microS", a2)
// According to equation 16.40
et = neta*(alpha*V0)/sqrt(alpha^2-1)
printf("\n e(t) = %f * (e^%ft - f^%ft)",et,-a1,-a2) // Equation of the wave form generated by the impulese

//Answers may vary due to round off error