dataCell=f1('tartalyrendszer.txt','text');
data=dataCell{1};
r1=data(1);
r2=data(2);
tspan=data(3:4);
y0=data(5:6);
fv=@(t,y) [-r1*y(1);r1*y(1)-r2*y(2)];
[T,Y]=ode45(fv,tspan,y0);
plot(T,Y(:,1),'bo-.','Linewidth',1.5)
hold on
plot(T,Y(:,2),'r+--','Linewidth',1.5)
plot([0,T(10)],[Y(10,1),Y(10,1)],'g','Linewidth',1.5)
plot([0,T(10)],[Y(10,2),Y(10,2)],'g','Linewidth',1.5)
plot([0,T(15)],[Y(15,2),Y(15,2)],'g','Linewidth',1.5)
plot([0,T(15)],[Y(15,1),Y(15,1)],'g','Linewidth',1.5)
plot([0,T(20)],[Y(20,1),Y(20,1)],'g','Linewidth',1.5)
kulonbseg=abs(Y(:,1)-Y(:,2));
metszespont=[T(min(kulonbseg)==kulonbseg),Y(min(kulonbseg)==kulonbseg)];
plot([metszespont(1) metszespont(1)],[1 metszespont(2)],'k--','Linewidth',1.5)
title('Tartály rendszer','FontSize',24)
set(gca,'Fontsize',12);
legend({'A tartály','B tartály'},'FontSize',12)
xlabel(sprintf('idõ [óra], metszéspont: %3.2f',metszespont(1)),'Fontsize',12,'fontweight','bold');
tartomany=T<=metszespont(1);
T1=trapz(T(tartomany),Y(tartomany,1));
T2=trapz(T(tartomany),Y(tartomany,2));
kozosTerulet=abs(T1-T2);
ylabel(sprintf('vízmennyiség [liter], terület: %3.2f',kozosTerulet),'Fontsize',12,'fontweight','bold');
fid=fopen('eredmeny.txt','w');
fprintf(fid,'Tartály rendszer\nMetszéspont: %3.5f\nKözös terület: %3.5f\n',metszespont(1),kozosTerulet);
fclose(fid);