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,