clear all
global R tau Q T0 p0 m V0
R=287.058;
tau=200.0;
Q=1;
T0=293.15;
p0=100000;
V0=0.001;
m=(p0*V0)/(R*T0);
simulation_opt=simset('Solver','ode5','FixedStep',1e-3);
Tsim=tau;
sim('Isotherm_mod',Tsim,simulation_opt);
figure(2);
plot(ans.Druck);
xlabel 'Zeit/s';
ylabel 'Druck/Pa';
title 'Isotherme Zustandsänderung'
function y=Gleichung_Volumen(u)
global V0 Q m R T0
t=u(1);
volumen=u(2);
y=volumen*(Q/(m*R*T0))+V0;
end
function y=Gleichung_Druck(u)
global V0 p0
t=u(1);
y=(p0*V0)/u;
end