function [iQ,iD] = park(iq,id,ang)
% This function computes the park transformation to change
% the stationary dq stator variables to synchronously rotating 
% dq stator variables 
% Inputs:  
%         iq = Stationary q-axis stator variable
%         id = Stationary q-axis stator variable
%         ang = Synchronously rotating angle between 0-2*pi (rad)
% Outputs:
%         iQ = Synchronously rotating q-axis stator variable
%         iD = Synchronously rotating d-axis stator variable

iQ = -id*sin(2*pi*ang)+iq*cos(2*pi*ang);
iD = id*cos(2*pi*ang)+iq*sin(2*pi*ang);