Hello,(adsbygoogle = window.adsbygoogle || []).push({});

I'm having some trouble with my octave coding and would appreciate any input on where the issue lies.

The coding is as follows:

<Moderator's note: code tags added. Please use them when posting code.>Code (Matlab M):

age = [0:1:100]; %this is the age matrix, represented by a

time = [0:1:100]; %this is the time matrix, represented by t

population = zeros(101,101);

for a = 1:101;

population(a,1) = 0.5; %this is the initial population

B(a) = (1./(2.*a)); %this is the birthing probability at age (a)

endfor

for t = 1:100

for a = 2:101

population(1,t+1) = cumtrapz(B(a)); %this is the birth rate

population(a,t+1) = population(a-1,t); %this is the interation process

endfor

endfor

surf(age,time,population)

The issue I'm having lies in the population matrix computation; the cumtrapz integration is not recognizing the step process and all outputs from t = 2:101 in the population matrix are zeros. I would again appreciate any help to understand where the issue lies.

Thank You,

Jacob M.

**Physics Forums | Science Articles, Homework Help, Discussion**

Join Physics Forums Today!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# MATLAB Octave integral computation help

Tags:

Have something to add?

Draft saved
Draft deleted

Loading...

Similar Threads for Octave integral computation |
---|

MATLAB Double Numerical Integration |

MATLAB Integrating Trisurf |

Mathematica Complex output from a real integral |

Maple Computing Numerical Integrals with Maple |

Mathematica Cannot do the integral of the Hyper-geometric function? |

**Physics Forums | Science Articles, Homework Help, Discussion**