CHEMOSTAT

Graf P(t), S(t)

Slika 4-6a  	4-6b

Clear[v,k,y,a1,k1,a2,k2,p,s,q];
v=1;
k=1;
y=0.99;
omega = 0.1;
s0=1;
xend=40;
eq=NDSolve[{
s'[x]==-v s[x]/(k+s[x])p[x]/y+ omega s0 - omega s[x],
p'[x]==v s[x]/(k+s[x])p[x]-omega p[x], 
p[0]==0.02,s[0]==2},{p[x],s[x]},{x,0,xend}];
p[x_]=N[p[x]]/.eq;
s[x_]=N[s[x]]/.eq;
g1=Plot[s[x],{x,0,xend},PlotRange->{0,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," S(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Red,Thickness->0.005},{Green,Thickness->0.005}}]
g2=Plot[p[x],{x,0,xend}, PlotRange->{0,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," P(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Blue,Thickness->0.005},{Green,Thickness->0.005}}]

n[x_]=p[x]+q[x];
g4=Plot[n[x],{x,0,xend}, PlotRange->{0,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," P(t)+Q(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Blue,Thickness->0.005},{Green,Thickness->0.005}}]

BIOFILM

Graf P(t), S(t), Q(t)

Slika 4-7a  	4-7b	4-7c

Clear[v,k,y,a1,k1,a2,k2,p,s,q];
v=1;
k=1;
y=0.99;
omega = 0.1;
s0=1;
vq=0.2;
kq=1;
a1=0.005;
k1=1;
a2=1;
k2=1;
xend=40;
eq=NDSolve[{
s'[x]==-v s[x]/(k+s[x])p[x]/y -vq s[x]/(kq+s[x])q[x]/y+ omega s0 - omega s[x],
p'[x]==v s[x]/(k+s[x])p[x]-omega p[x] -a1/(k1+s[x]) p[x]+ a2 s[x]/(k2+s[x]) q[x], 
q'[x]==vq s[x]/(kq+s[x])q[x]+a1/(k1+s[x]) p[x]- a2 s[x]/(k2+s[x]) q[x], 
p[0]==0.02,s[0]==2,q[0]==0.01},{p[x],s[x],q[x]},{x,0,xend}];
p[x_]=N[p[x]]/.eq;
s[x_]=N[s[x]]/.eq;
q[x_]=N[q[x]]/.eq;
g1=Plot[s[x],{x,0,xend},PlotRange->{0,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," S(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Red,Thickness->0.005},{Green,Thickness->0.005}}]
g2=Plot[p[x],{x,0,xend}, PlotRange->{0,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," P(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Blue,Thickness->0.005},{Green,Thickness->0.005}}]
g3=Plot[q[x],{x,0,xend}, PlotRange->{0,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," Q(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Blue,Thickness->0.005},{Green,Thickness->0.005}}]

n[x_]=p[x]+q[x];
g4=Plot[n[x],{x,0,xend}, PlotRange->{0,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," P(t)+Q(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Blue,Thickness->0.005},{Green,Thickness->0.005}}]

konstantni alpha i beta

Clear[v,k,y,a1,k1,a2,k2,p,s,q];
v=1;
k=1;
y=0.99;
omega = 0.1;
s0=1;
vq=0.2;
kq=1;
a1=0.05;
k1=1;
a2=1;
k2=1;
xend=40;
eq=NDSolve[{
s'[x]==-v s[x]/(k+s[x])p[x]/y -vq s[x]/(kq+s[x])q[x]/y+ omega s0 - omega s[x],
p'[x]==v s[x]/(k+s[x])p[x]-omega p[x] -a1 p[x]+ a2  q[x], 
q'[x]==vq s[x]/(kq+s[x])q[x]+a1 p[x]- a2  q[x], 
p[0]==0.02,s[0]==2,q[0]==0.01},{p[x],s[x],q[x]},{x,0,xend}];
p[x_]=N[p[x]]/.eq;
s[x_]=N[s[x]]/.eq;
q[x_]=N[q[x]]/.eq;
g1=Plot[s[x],{x,0,xend},PlotRange->{0,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," S(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Red,Thickness->0.005},{Green,Thickness->0.005}}]
g2=Plot[p[x],{x,0,xend}, PlotRange->{0,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," P(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Blue,Thickness->0.005},{Green,Thickness->0.005}}]
g3=Plot[q[x],{x,0,xend}, PlotRange->{0,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," Q(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Blue,Thickness->0.005},{Green,Thickness->0.005}}]

n[x_]=p[x]+q[x];
g4=Plot[n[x],{x,0,xend}, PlotRange->{0,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," P(t)+Q(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Blue,Thickness->0.005},{Green,Thickness->0.005}}]

HRANIDBENI LANAC s ograničenjem bez mirovanja

Graf P1(t), P2(t), P3(t)

Slika 4-8a  	4-8b	4-8c

Clear[v,k,y,a1,k1,a2,k2,p1,s1,p2,s2,p3,s3];
v1=1;
k1=1;
y1=1;
v2=1;
k2=1;
y2=1;
v3=1;
k3=1;
y3=1;
m12=1;
m23=1;
omega = 0.1;
s0=1;
vq=0.2;
kq=1;
a1=0.005;
k1=1;
a2=1;
k2=1;
xend=40;
eq=NDSolve[{
s1'[x]==-v1 s1[x]/(k1+s1[x])p1[x]/y1 -m12 p1[x],
p1'[x]==v1 s1[x]/(k1+s1[x])p1[x], 
s2'[x]==-v2 s2[x]/(k2+s2[x])p2[x]/y2 +m12 p1[x]-m23 p2[x],
p2'[x]==v2 s2[x]/(k2+s2[x])p2[x], 
s3'[x]==-v3 s3[x]/(k3+s3[x])p3[x]/y3 +m23 p2[x],
p3'[x]==v3 s3[x]/(k3+s3[x])p3[x], 
p1[0]==0.02,s1[0]==1,s2[0]==0,p2[0]==0.01,s3[0]==0,p3[0]==0.01},{s1[x],p1[x],s2[x],p2[x],s3[x],p3[x]},{x,0,xend}];
p1[x_]=N[p1[x]]/.eq;
p2[x_]=N[p2[x]]/.eq;
p3[x_]=N[p3[x]]/.eq;
g1=Plot[p1[x],{x,0,xend},PlotRange->{0,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," (Subscript[P, 1])(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Blue,Thickness->0.005},{Green,Thickness->0.005}}]
g2=Plot[p2[x],{x,0,xend}, PlotRange->{0,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," (Subscript[P, 2])(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Blue,Thickness->0.005},{Green,Thickness->0.005}}]
g3=Plot[p3[x],{x,0,xend}, PlotRange->{0,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," (Subscript[P, 3])(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Blue,Thickness->0.005},{Green,Thickness->0.005}}]

s1[x_]=N[s1[x]]/.eq;
s2[x_]=N[s2[x]]/.eq;
s3[x_]=N[s3[x]]/.eq;
g1=Plot[s1[x],{x,0,xend},PlotRange->{-0.5,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," (Subscript[S, 1])(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Red,Thickness->0.005},{Green,Thickness->0.005}}]
g2=Plot[s2[x],{x,0,xend}, PlotRange->{-0.5,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," (Subscript[S, 2])(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Red,Thickness->0.005},{Green,Thickness->0.005}}]
g3=Plot[s3[x],{x,0,xend}, PlotRange->{-0.5,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," (Subscript[S, 3])(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Red,Thickness->0.005},{Green,Thickness->0.005}}]

HRANIDBENI LANAC s ograničenjem bez mirovanja modifikacija potrošnje

Graf P1(t), P2(t), P3(t)

Slika 4-8a  	4-8b	4-8c

Clear[v,k,y,a1,k1,a2,k2,p1,s1,p2,s2,p3,s3];
v1=1;
k1=1;
y1=1;
v2=1;
k2=1;
y2=1;
v3=1;
k3=1;
y3=1;
m12=1;
m23=1;
omega = 0.1;
s0=1;
vq=0.2;
kq=1;
a1=0.005;
k1=1;
a2=1;
k2=1;
xend=50;
eq=NDSolve[{
s1'[x]==-v1 s1[x]/(k1+s1[x])p1[x]/y1 -s1[x]/(1+s1[x]) m12 p1[x],
p1'[x]==v1 s1[x]/(k1+s1[x])p1[x], 
s2'[x]==-v2 s2[x]/(k2+s2[x])p2[x]/y2 +s1[x]/(1+s1[x])m12 p1[x]-s2[x]/(1+s2[x])m23 p2[x],
p2'[x]==v2 s2[x]/(k2+s2[x])p2[x], 
s3'[x]==-v3 s3[x]/(k3+s3[x])p3[x]/y3 +s2[x]/(1+s2[x])m23 p2[x],
p3'[x]==v3 s3[x]/(k3+s3[x])p3[x], 
p1[0]==0.02,s1[0]==1,s2[0]==0,p2[0]==0.01,s3[0]==0,p3[0]==0.01},{s1[x],p1[x],s2[x],p2[x],s3[x],p3[x]},{x,0,xend}];
p1[x_]=N[p1[x]]/.eq;
p2[x_]=N[p2[x]]/.eq;
p3[x_]=N[p3[x]]/.eq;
g1=Plot[p1[x],{x,0,xend},PlotRange->{0,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," (Subscript[P, 1])(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Blue,Thickness->0.005},{Green,Thickness->0.005}}]
g2=Plot[p2[x],{x,0,xend}, PlotRange->{0,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," (Subscript[P, 2])(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Blue,Thickness->0.005},{Green,Thickness->0.005}}]
g3=Plot[p3[x],{x,0,xend}, PlotRange->{0,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," (Subscript[P, 3])(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Blue,Thickness->0.005},{Green,Thickness->0.005}}]

s1[x_]=N[s1[x]]/.eq;
s2[x_]=N[s2[x]]/.eq;
s3[x_]=N[s3[x]]/.eq;
g1=Plot[s1[x],{x,0,xend},PlotRange->{-0.5,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," (Subscript[S, 1])(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Red,Thickness->0.005},{Green,Thickness->0.005}}]
g2=Plot[s2[x],{x,0,xend}, PlotRange->{-0.5,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," (Subscript[S, 2])(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Red,Thickness->0.005},{Green,Thickness->0.005}}]
g3=Plot[s3[x],{x,0,xend}, PlotRange->{-0.5,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," (Subscript[S, 3])(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Red,Thickness->0.005},{Green,Thickness->0.005}}]

HRANIDBENI LANAC s ograničenjem s mirovanja

Graf P1(t), P2(t), P3(t)

Slika 4-9a  	4-9b	4-9c

Clear[v,k,y,a1,k1,a2,k2,p,q,s,p1,s1,p2,s2,p3,s3,q1,q2,q3,t];
v1=1;
k1=1;
y1=1;
v2=1;
k2=1;
y2=1;
v3=1;
k3=1;
y3=1;
m12=1;
m23=1;
a1=1;
k1=1;
a11=1;
k11=1;
a2=1;
k2=1;
a22=1;
k22=1;
a3=1;
k3=1;
a33=1;
k33=1;
v=1;
k=1;
y=0.99;
a1=1;
k1=1;
a2=1;
k2=1;
xend=40;

eq=NDSolve[{
s1'[x]==-v1 s1[x]/(k1+s1[x])p1[x]/y1 -m12 p1[x],
p1'[x]==v1 s1[x]/(k1+s1[x])p1[x]-a1/(k1+s1[x]) p1[x]+ a11 s1[x]/(k11+s1[x]) q1[x], 
q1'[x]==a1/(k1+s1[x]) p1[x]- a11 s1[x]/(k11+s1[x]) q1[x],
p1[0]==0.02,s1[0]==3,q1[0]==0.001},{p1[x],s1[x],q1[x]},{x,0,xend}];
t[x_]=N[p1[x]]/.eq;
p1[x_]=t[x][[1]];
s1[x_]=N[s1[x]]/.eq;
g1=Plot[p1[x],{x,0,xend},PlotRange->{0,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," (Subscript[P, 1])(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Blue,Thickness->0.005},{Green,Thickness->0.005}}]
eq2=NDSolve[{
s2'[x]==-v2 s2[x]/(k2+s2[x])p2[x]/y2 +m12 p1[x] - m23 p2[x],
p2'[x]==v2 s2[x]/(k2+s2[x])p2[x]-a2/(k2+s2[x]) p2[x]+ a22 s2[x]/(k22+s2[x]) q2[x], 
q2'[x]==a2/(k2+s2[x]) p2[x]- a22 s2[x]/(k22+s2[x]) q2[x],
s2[0]==0,p2[0]==0.01,q2[0]==0.001},{s2[x],p2[x],q2[x]},{x,0,xend}];
t2[x_]=N[p2[x]]/.eq2;
p2[x_]=t2[x][[1]];
s2[x_]=N[s2[x]]/.eq2;

g2=Plot[p2[x],{x,0,xend}, PlotRange->{0,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," (Subscript[P, 2])(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Blue,Thickness->0.005},{Green,Thickness->0.005}}]
eq=NDSolve[{
s3'[x]==-v3 s3[x]/(k3+s3[x])p3[x]/y3 +m23 p2[x],
p3'[x]==v3 s3[x]/(k3+s3[x])p3[x]-a3/(k3+s3[x]) p3[x]+ a33 s3[x]/(k33+s3[x]) q3[x], 
q3'[x]==a3/(k3+s3[x]) p3[x]- a33 s3[x]/(k33+s3[x]) q3[x],
s3[0]==0,p3[0]==0.01,q3[0]==0.001},{s3[x],p3[x],q3[x]},{x,0,xend}];
p3[x_]=N[p3[x]]/.eq;
s3[x_]=N[s3[x]]/.eq;
g3=Plot[p3[x],{x,0,xend}, PlotRange->{0,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," (Subscript[P, 3])(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Blue,Thickness->0.005},{Green,Thickness->0.005}}]

h1=Plot[s1[x],{x,0,xend},PlotRange->{-0.5,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," (Subscript[S, 1])(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Red,Thickness->0.005},{Green,Thickness->0.005}}]
h2=Plot[s2[x],{x,0,xend}, PlotRange->{-0.5,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," (Subscript[S, 2])(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Red,Thickness->0.005},{Green,Thickness->0.005}}]
h3=Plot[s3[x],{x,0,xend}, PlotRange->{-0.5,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," (Subscript[S, 3])(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Red,Thickness->0.005},{Green,Thickness->0.005}}]

HRANIDBENI LANAC s ograničenjem s mirovanja

Graf P1(t), P2(t), P3(t)

Slika 4-10a  	4-10b	4-10c

Clear[v,k,y,a1,k1,a2,k2,p,q,s,p1,s1,p2,s2,p3,s3,q1,q2,q3,t];
v1=1;
k1=1;
y1=1;
v2=1;
k2=1;
y2=1;
v3=1;
k3=1;
y3=1;
m12=1;
m23=1;
a1=1;
k1=1;
a11=1;
k11=1;
a2=1;
k2=1;
a22=1;
k22=1;
a3=1;
k3=1;
a33=1;
k33=1;
v=1;
k=1;
y=0.99;
a1=1;
k1=1;
a2=1;
k2=1;
xend=40;

eq=NDSolve[{
s1'[x]==-v1 s1[x]/(k1+s1[x])p1[x]/y1 -m12 s1[x]/(1+s1[x])p1[x],
p1'[x]==v1 s1[x]/(k1+s1[x])p1[x]-a1/(k1+s1[x]) p1[x]+ a11 s1[x]/(k11+s1[x]) q1[x], 
q1'[x]==a1/(k1+s1[x]) p1[x]- a11 s1[x]/(k11+s1[x]) q1[x],
p1[0]==0.02,s1[0]==6,q1[0]==0.001},{p1[x],s1[x],q1[x]},{x,0,xend}];
t[x_]=N[p1[x]]/.eq;
p1[x_]=t[x][[1]];
t[x_]=N[s1[x]]/.eq;
s1[x_]=t[x][[1]];
g1=Plot[p1[x],{x,0,xend},PlotRange->{0,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," (Subscript[P, 1])(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Blue,Thickness->0.005},{Green,Thickness->0.005}}]
eq2=NDSolve[{
s2'[x]==-v2 s2[x]/(k2+s2[x])p2[x]/y2 +m12 s1[x]/(1+s1[x])p1[x] - m23 s2[x]/(1+s2[x])p2[x],
p2'[x]==v2 s2[x]/(k2+s2[x])p2[x]-a2/(k2+s2[x]) p2[x]+ a22 s2[x]/(k22+s2[x]) q2[x], 
q2'[x]==a2/(k2+s2[x]) p2[x]- a22 s2[x]/(k22+s2[x]) q2[x],
s2[0]==0,p2[0]==0.01,q2[0]==0.001},{s2[x],p2[x],q2[x]},{x,0,xend}];
t2[x_]=N[p2[x]]/.eq2;
p2[x_]=t2[x][[1]];
t2[x_]=N[s2[x]]/.eq2;
s2[x_]=t2[x][[1]];

g2=Plot[p2[x],{x,0,xend}, PlotRange->{0,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," (Subscript[P, 2])(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Blue,Thickness->0.005},{Green,Thickness->0.005}}]
eq=NDSolve[{
s3'[x]==-v3 s3[x]/(k3+s3[x])p3[x]/y3 +m23 s2[x]/(1+s2[x])p2[x],
p3'[x]==v3 s3[x]/(k3+s3[x])p3[x]-a3/(k3+s3[x]) p3[x]+ a33 s3[x]/(k33+s3[x]) q3[x], 
q3'[x]==a3/(k3+s3[x]) p3[x]- a33 s3[x]/(k33+s3[x]) q3[x],
s3[0]==0,p3[0]==0.01,q3[0]==0.001},{s3[x],p3[x],q3[x]},{x,0,xend}];
p3[x_]=N[p3[x]]/.eq;
s3[x_]=N[s3[x]]/.eq;
g3=Plot[p3[x],{x,0,xend}, PlotRange->{0,2.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," (Subscript[P, 3])(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Blue,Thickness->0.005},{Green,Thickness->0.005}}]

h1=Plot[s1[x],{x,0,xend},PlotRange->{-0.5,3.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," (Subscript[S, 1])(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Red,Thickness->0.005},{Green,Thickness->0.005}}]
h2=Plot[s2[x],{x,0,xend}, PlotRange->{-0.5,3.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," (Subscript[S, 2])(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Red,Thickness->0.005},{Green,Thickness->0.005}}]
h3=Plot[s3[x],{x,0,xend}, PlotRange->{-0.5,3.1},AspectRatio->0.4,FrameLabel->{"Vrijeme"," (Subscript[S, 3])(t)"},Frame->{True,True,False,False},FrameStyle->{Directive[12,Thickness[0.004]],Directive[12,Thickness[0.004]]},PlotStyle->{{Red,Thickness->0.005},{Green,Thickness->0.005}}]