function f = Eg(n,a)  % a = alpha
format long g;
Ezero = 15000; % Einheit in Bar
Ezero_inv = (1/Ezero)
m = 10
pzero = 0
m_inv = (-1/m)
n_inv = (1/n)
m_num = (m+1)*m_inv
n_num = (n+1)*n_inv



 p = 0:1:5;
 
 d = length(p);

for i = 2:1:d

   
 Zaehler(i) = (1-a)*(1+(m*p(i))/Ezero)^m_inv + a*(pzero /p(i))^n_inv
Nenner(i) = ((Ezero_inv)*(1-a)*(1+((m*p(i))/Ezero))^m_num )+((a/(n*pzero))*(pzero /p(i))^n_num)

y(i) = Zaehler(i) /Nenner(i)

plot(p,y(i));
 end
 


% hold on
%plot(n(900:end),[ew(900:end), elms(900:end),enlms(900:end)]);
%legend('Wiener filter denoised sinusoid',...
 %   'LMS denoised sinusoid', 'NLMS denoised sinusoid');
%xlabel('Time index (n)');
%ylabel('Amplitude');
% hold off