need to have graph for EA (0.9,0.8,0.7,0.6,0.5) on same axis or same graph(in linearly increasing shape)
Show older comments
clear all;close all;clc; k=8.617385*10^(-5); T0=298; EA=0.9; Ta=[25,30,35,40,45,50,55,60,65,70,75,80,85,90,95]; l=length(Ta); n=(EA/k)*((1/T0)-(1/(298))); T=5; for i=1:l
n=(EA/k)*((1/T0)-(1/((i*T)+293)));
AF(i)=exp(n);
end
disp(AF(i))
k=8.617385*10^(-5);
T0=298;
EA=0.8;
Ta=[25,30,35,40,45,50,55,60,65,70,75,80,85,90,95];
l=length(Ta);
n=(EA/k)*((1/T0)-(1/(298)));
T=5;
for j=1:l
n=(EA/k)*((1/T0)-(1/((j*T)+293)));
AF(j)=exp(n);
end
disp(AF(j))
k=8.617385*10^(-5);
T0=298;
EA=0.7;
Ta=[25,30,35,40,45,50,55,60,65,70,75,80,85,90,95];
l=length(Ta);
n=(EA/k)*((1/T0)-(1/(298)));
T=5;
for a=1:l
n=(EA/k)*((1/T0)-(1/((a*T)+293)));
AF(a)=exp(n);
end
disp(AF(a))
k=8.617385*10^(-5);
T0=298;
EA=0.6;
Ta=[25,30,35,40,45,50,55,60,65,70,75,80,85,90,95];
l=length(Ta);
n=(EA/k)*((1/T0)-(1/(298)));
T=5;
for b=1:l
n=(EA/k)*((1/T0)-(1/((b*T)+293)));
AF(b)=exp(n);
end
disp(AF(b))
k=8.617385*10^(-5);
T0=298;
EA=0.5;
Ta=[25,30,35,40,45,50,55,60,65,70,75,80,85,90,95];
l=length(Ta);
n=(EA/k)*((1/T0)-(1/(298)));
T=5;
for c=1:l
n=(EA/k)*((1/T0)-(1/((c*T)+293)));
AF(c)=exp(n);
end
disp (AF(c))
plot(Ta,AF(i),':bs','LineWidth',2);
hold on
plot(Ta,AF(j),'--mo','LineWidth',2);
plot(Ta,AF(a),'-.r*','LineWidth',2);
plot(Ta,AF(b),'-.gd','LineWidth',2);
plot(Ta,AF(c),'-c+','LineWidth',2);
grid on;
hold off
xlim([0 100]);
ylim([0 800]);
Answers (0)
Categories
Find more on Startup and Shutdown in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!