%UKŁAD RL/ TEMAT2
R=input('R=')
L1=input('L1=')
L2=input('L2=')
L3=input('L3=')

tau1=L1/R;
tau2=L2/R;
tau3=L3/R;

t1=0:tau1/10:10.*tau1;
t2=0:tau2/10:10.*tau2;
t3=0:tau3/10:10.*tau3;

it1=1/L1.*exp(-t1/tau1);
it2=1/L2.*exp(-t2/tau2);
it3=1/L3.*exp(-t3/tau3);

ur1= (R/L1).*exp(-t1/tau1);
ur2= (R/L2).*exp(-t2/tau2);
ur3= (R/L3).*exp(-t3/tau3);

ult1= ((-R/L1).*exp(-t1/tau1));
ult2= ((-R/L2).*exp(-t2/tau2));
ult3= ((-R/L3).*exp(-t3/tau3));

figure;
hold on;

subplot(3,1,1);
plot(t1,it1,'--',t2,it2,'-',t3,it3,'-.');
xlabel('t[s]');
ylabel('i[A]');
title('prad')
legend('i1(dla L1=0.001)','i2(dla L2=0.002)','i3(dla L3=0.003)');

subplot(3,1,2);
plot(t1,ur1,'--',t2,ur2,'-',t3,ur3,'-.');
xlabel('t[s]');
ylabel('Ur[V]');
legend(' Rezystor [1]',' Rezystor [2]',' Rezystor [3]');

subplot(3,1,3);
plot(t1,ult1,'--',t2,ult2,'-',t3,ult3,'-.');
xlabel('t[s]');
ylabel('Ul[V]');
title('prad')
legend('Cewka [1]','Cewka [2]','Cewka [3]');