why is showing Inf and Nan value how to remove that

a(1)=0;

b(1)=1;

c(1)=0;

k3=28;

k1=8/3;

k2=10;

s1=0;

s2=0;

s3=0;

T=0.1;

t=0:0.01:0.1;

tau = t ./ T;

for i=1:20

a(i+1)=(1/(i+1)).*(b(i)-a(i))*k2;

b(i+1)=(1/(i+1)).*(T*a(i).*(k3-c(i))-b(i));

c(i+1)=(1/(i+1))*(T*a(i).*b(i)-k1*c(i));

s1=s1+(a(i).*tau.^i);

s2=s2+(b(i).*tau.^i);

s3=s3+(c(i).*t.au^i);

disp([s1;s2;s3])

end

plot(t, s1, t, s2, t, s3)

disp is showing INF NAn

[Benjamin Thompson]

Your sample code does not run. Can you correct it? This appears to be the line.

s3=s3+(c(i).*t.au^i);

[Jan]

@shiv gaur: You have been instructed to apply a proper code formatting. Demonstrating, that you ignore hints of the community, can discourages readers to post answers.

As soon as you apply a proper formatting, you can even run the code in the browser and then you would see, that it stops with an error.

[shiv gaur]

a(1)=0;

b(1)=1;

c(1)=0;

k3=28;

k1=8/3;

k2=10;

s1=0;

s2=0;

s3=0;

T=0.1;

t=0:0.01:0.1;

tau = t ./ T;

for i=1:20

a(i+1)=(1/(i+1)).*(b(i)-a(i))*k2;

b(i+1)=(1/(i+1)).*(T*a(i).*(k3-c(i))-b(i));

c(i+1)=(1/(i+1))*(T*a(i).*b(i)-k1*c(i));

s1=s1+(a(i).*tau.^i);

s2=s2+(b(i).*tau.^i);

s3=s3+(c(i).*tau.^i);

disp([s1;s2;s3])

end

plot(t, s1, t, s2, t, s3)

this is the format program display is showing like that

NaN Inf Inf Inf Inf Inf Inf Inf Inf Inf Inf

NaN -Inf -Inf -Inf -Inf -Inf -Inf -Inf -Inf -Inf -Inf

NaN Inf Inf Inf Inf Inf Inf Inf Inf Inf Inf

NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN

NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN

>>

[Mathieu NOE]

hello

@Benjamin Thompson

there is a small typo in the posted code - this should be :

s3=s3+(c(i).*tau.^i); % here

like s1 and s2 above

@shiv gaur

seems a,b,and c are diverging inside the for loop, their values are skyrocketing towards very big numbers up to the point it's Inf or NaN

check again your recursion equations and data initilaization

[Mathieu NOE]

hello again

simply testing the code for 10 iterations shows the diverging trend...

what was the expected result ?

[Jan]

@shiv gaur: You have been asked repeatedly to use a proper code formatting. Please read and consider this: https://www.mathworks.com/matlabcentral/answers/help/rtc#rtc_summary . Thanks.

Answer 1

[Jan]

Your function is growing exponentially. It is the mathematically correct behaviour, that it reaches Inf. "Fixing" this would mean, to compute something else, but this is not a meaningful solution.