function [U,u]=kalkulator(y);
%% input
%% vector y - absolute values of coefficients of Fourier series 
%%      (e.g voltage signal), number of coefficients is not known.
%% outputs:
%%    U - signal value vs. time, 120 values per time cycle of base frequency
%%    u - RMS value of the above signal  (square root of sum of squares)
%%
% obliczenia dla przebiegow odkszta³conych
% y    wielomian o nieznanej d³ugoœci, pomierzone wartoœci skuteczne napiêcia 
%                              o czêstotliwoœci podstawowej oraz jego harmoniczne 
% U    wektor funkcji czasu napiêcia odkszta³conego,  ?  Uk*sin (k*?*t).
%                                ma mieæ 120 punktów i obejmuje 2 okresy napiêcia 
% u    liczba, wartoœæ skuteczna z amplitud sk³adowych harmonicznych 
%>>[U,u]=kalkulator([1 0 1/3 0 1/5 0 1/7]); figure, plot(U) %wywolanie 

u=sqrt(y*y')  %jesli wektor pionowy, to y'*y 

n=length(y)
koniec=2*2*pi; krok=koniec/120
psi=krok:krok:koniec % 120 punktow (bez zera na pocz¹tku)
for i=1:n
   UU(i,:)=y(i)*sin(i*psi);
end
U=sum(UU)
plot(psi,UU') % veryfication only  %niepotrzebne, to tylko sprawdzenie,
whos          % veryfication only  %niepotrzebne, to tylko sprawdzenie,
