function konkka(A,B,r,ui,uo)

% funktio konkka(A,B,r,ui,uo) laskee virtalähteen rippelin
% suotokondensaattorin navoissa
% Parametrit:
% A=suotokondensaattorin koko uf
% B=kuormittava virta mA
% r=virtalähteen sisäinen vast  ohm
% ui=vaihtojännitteen huippuarvo V
% uo=kondensaattorin DC-jännite alkutilanteessa V

	c=A/1e6;
	i=B/1e3;   
	n=50;
	dt=0.02/n;
	t=0;
    tout=[0];
    uin=[0];
    uout=[0];

for k=1:1000,
	t=t+dt;
	tout=[tout t];
	u= abs(ui*sin(2*pi*50*t));
	uin=[uin u];
	uout=[uout uo];
	uch=u-uo;
	if uch>0,
		uc=uo-(i*dt)/c+dt*uch/(r*c);
		if uc>u,
			uo=u;
		else
			uo=uc;
		end
	else
		uc=uo-i*dt/c;
		uo=uc;
	end
end

	plot(tout,uin,tout,uout);
