Circuit Theory/Phasors/Examples/Example 7/mfile7a

format short
R = 10;
L = .01;
w = 377;
Vm = 120*2^.5;
omega = 2*pi/3;
Im = Vm/(R^2+(L*w)^2)^.5;
alpha = omega - atan(L*w/R);
syms t;
I = Im * cos(w*t+alpha);
Vr = R*I;
VL = L*diff(I,t);
fprintf('I is')
vpa(I,3)
fprintf('Vr is')
vpa(Vr,3)
fprintf('VL is')
vpa(VL,3)
fprintf('checking .. should equal 0')
subs(Vr,t,0) + subs(VL,t,0) - Vm*cos(omega)