Matlab help

  • MATLAB
  • Thread starter roam
  • Start date
  • #1
1,266
11

Main Question or Discussion Point

I want to find q501 for the difference equation:

qn=qn-1 - 0.001 q2n-2
n≥2,
q0=1,
q1=2

I thought it's easier to do this on Matlab & I tried to write a script file for this problem:

Code:
k=0.001
q(1)=1
q(2)=2
for n=1:502
    q_n=q_n-1-k*(q_n-2)^2
end
I don't know why my code isn't working/giving the right answer. Does anyone know what's the problem with this script?
 

Answers and Replies

  • #2
f95toli
Science Advisor
Gold Member
2,998
493
You are not indexing the array in the right way
Element n of a vector q is written q(n) in Matlab (with the first element having index 1, not zero as in some other programming languages).

Try
Code:
k=0.001
q(1)=1
q(2)=2
for n=3:502
    q(n)=q(n-1)-1-k*q(n-2)^2
end
 

Related Threads on Matlab help

  • Last Post
Replies
2
Views
3K
  • Last Post
Replies
2
Views
2K
  • Last Post
Replies
7
Views
3K
  • Last Post
Replies
3
Views
4K
  • Last Post
Replies
2
Views
2K
  • Last Post
Replies
3
Views
2K
  • Last Post
Replies
6
Views
7K
  • Last Post
Replies
9
Views
3K
  • Last Post
Replies
1
Views
2K
  • Last Post
Replies
6
Views
2K
Top