I am getting error:Warning: Failure at t=0.000000e+00. Unable to ... (2025)

clc

clear

global datei Temp c0

datei='input_data.m';

run(datei);

zspan=[0 V/A_hc];

%intergration of the system

options=odeset('RelTol',1e-12,'AbsTol',1e-12);%,'NonNegative',2:7);

for n=1:size(c,1)

switch dataset

case 1

name=[experiment ' ' num2str(y0_nh3_scr(n,2)) ' ppm NH3 ' num2str(y0_nh3_scr(n,1)) ' vol NO ' num2str(y0_nh3_scr(n,5)) '% O2' ];

case 2

name=[experiment ' ' num2str(y0_nh3_ox(n,2)) ' ppm NH3 ' num2str(y0_nh3_ox(n,5)) ' vol O2'];

case 3

name=[experiment ' ' num2str(y0_hcho_ox(n,3)) ' ppm HCHO ' num2str(y0_hcho_ox(n,5)) ' vol O2'];

end

results=zeros(length(T),size(c,2)+2);

for m=1:length(T)

c0=c(n,1:4).*T_N./T(m);

c_o2=c0_o2(n).*T_N./T(m);

p_o2=p0_o2(n);

y0=[c0 0];

Temp=T(m);

tic

disp(['Calculating Step ' num2str(m) ' of ' num2str(length(T))])

[t,y]=ode15s(@massbalance,zspan,y0,options);

toc

results(m,:)=[T(m)-273.15 y(end,:)];

end

end

results(:,2:5)=results(:,2:5).*Vm.*1e6.*(results(:,1)+273.15)./T_N;

switch dataset

case 1

results(:,2:3)=(1-results(:,2:3)./y0_nh3_scr(n,1:2)).*100;

case 2

results(:,3)=(1-results(:,3)./y0_nh3_ox(n,2)).*100;

case 3

results(:,4)=(1-results(:,4)./y0_hcho_ox(n,3)).*100;

end

name=[name ' ' char(num2str(fix(clock)))];

save(name,'results','-double','-ascii')

plot(results(:,1),results(:,2),'b',results(:,1),results(:,3),'r',results(:,1),results(:,6),'g')

xlabel('T / ^oC')

ylabel('y / ppm')

legend('NO','NH3','N2O')

function dy_dz = massbalance(takt,y)

global datei Temp c0

run(datei)

F=F*Temp/T_N;

u=F/A_hc;

if c0(1)==0

kin=kin_ox;

else

kin=kin_scr;

end

r1=kin(1,1).*exp(-kin(1,2)./R.*(1/Temp-1/kin(1,3))).*y(2).*(1-y(5));

r2=kin(2,1).*exp(-(kin(2,2)-y(5).*kin(2,3))./(R.*Temp)).*y(5);

r3=kin(3,1).*exp(-kin(3,2)./(R.*Temp)).*y(1).*y(5);

r4=kin(4,1).*exp(-kin(4,2)./(R.*Temp)).*y(5);

r5=kin(5,1).*exp(-kin(5,2)./(R.*Temp)).*y(4);

r6=kin(6,1).*exp(-kin(6,2.)/(R.*Temp)).*y(3);

%% mass balances

dy_dz=[1/u*(-r3)

1/u*(-r1+r2)

1/u*(-r6)

1/u*(-r5)

(r1-r2-r3-r4)/(u*gam)];

end