Index exceeds array bounds.

clc; clear;
dt = 0.5; 
Glucose = zeros(1,361);
Glucose(1,1) = 6000; 
Glucose_Released = 110;
Glucose_Used = zeros(1,361);
Usage_Fraction = zeros(1,361);
t = 0:dt:180; % We need to generate a time vector for our plot
for k = 2:length(t)
    Glucose_Used(k) = Glucose(k) * Usage_Fraction(k); % Based on Euler?sformula
    Glucose(k-1) = Glucose(k) + [Glucose_Released(k) - Glucose_Used(k)]*dt;
end

Please help!!

[Sudhaunsh Deshpande]

@Adam Danz Here it is

[Adam Danz]

Question edited to format code.

@Sudhaunsh Deshpande, looks like David hit the nail on the head. Let us know if you have any other questions.

[Sudhaunsh Deshpande]

It did kinda work, the values are now zero which means that my guy is dead, but hey atleast the damn thing runs!

Answer 1

[David Fletcher]

You have defined Glucose_Released as a scaler value

Glucose_Released = 110;

You are then trying to index it, which will cause an error

Glucose_Released(k)

[Sudhaunsh Deshpande]

The thing is Glucose_Released is the rate of glucose being released.

So the only option is:

Glucose_Released a matrix of (361X1) with the value 110

is that what you mean?

[David Fletcher]

It is your calculation - I am merely pointing out that trying to index a scaler is the cause of your error. If the rate of glucose release is invariant, do you need to index it at all - or does Glucose_Released really represent the amount of glucose released over time? I don't know - as I say I have no idea what you are doing, you asked for the source of the error: I provided it.

[Sudhaunsh Deshpande]

Okay that actually makes sense

thank you :)

[ellierose]

so how did you modify the script in the end to make it work?

[Sudhaunsh Deshpande]

Pardon my text editing:)

for i=1: length (t) Glucose Used(i) = Glucose(i) * Usage_Fraction(i); Glucose(i+1) = Glucose(i) + (G lucose_Released(i) Glucose_Used(i)) *dt; Insulin breakdown(i) = Insulin(i)/k; Insulin(i+1) Insulin(1) [Insulin_Secreted(i) Insulin_breakdown (1)] *dt;

[Walter Roberson]

@Sudhaunsh Deshpande

You have

Glucose(i+1) = Glucose(i) + (G lucose_Released(i) Glucose_Used(i)) *dt

We can assume an accidental space,

Glucose(i+1) = Glucose(i) + (Glucose_Released(i) Glucose_Used(i)) *dt

but there needs to be something between Glucose_Released(i) and Glucose_Used(i)