TS L6

function ster_6
close all;
tspan = [0.0  10.0];
x0 = [1.5  0 0 7];

u = 5;

options_ode = odeset('RelTol',1e-6,'AbsTol',1e-6);
[tsol,xsol] = ode45(@(t,x) peniciline(t,x,u), tspan, x0, options_ode);

x2 = xsol(:,2);
x4 = xsol(:,4);
wskaznik_jakosci = x2(end)*x4(end)

function dx = peniciline(t,x,u)
dx = zeros(4,1);
h1 = 0.11*(x(3)/(0.006*x(1) + x(3)));
h2 = 0.0055*(x(3)/(0.0001 + x(3)*(1 + 10*x(3))));
dx(1) = h1*x(1) - u*(x(1)/(500*x(4)));
dx(2) = h2*x(1) - 0.01*x(2) - u*(x(2)/(500*x(4)));
dx(3) = -h1*x(1)/0.47 - h2*x(1)/1.2 - x(1)*((0.029*x(3))/(0.0001 + x(3)));
dx(4) = u/500;
end