Problem 255. Spiraling out of control....

Created by Matt Fig

Solve Later

{"relationships":[{"relationshipType":"http://schemas.mathworks.com/matlab/code/2013/relationships/document","targetMode":"","relationshipId":"rId1","target":"/matlab/document.xml"},{"relationshipType":"http://schemas.mathworks.com/matlab/code/2013/relationships/output","targetMode":"","relationshipId":"rId2","target":"/matlab/output.xml"}],"parts":[{"partUri":"/matlab/document.xml","relationship":[],"contentType":"application/vnd.mathworks.matlab.code.document+xml","content":"<?xml version=\"1.0\" encoding=\"UTF-8\"?>\n<w:document xmlns:w=\"http://schemas.openxmlformats.org/wordprocessingml/2006/main\"><w:body><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>Imagine a square with corners at the origin and at (1,1). An interesting spiral can be created by tracing the path of a moving point that starts at the origin and intersects the leftmost opposite edge at a point (F,1), where F is the fraction of the length of the intersected edge. The point then continues on in this way, always heading for the leftmost opposite side of its path and intersecting the segment at a point which is fraction F of its length. For example, if F = 1/2, the first few points are:</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"code\"/></w:pPr><w:r><w:t><![CDATA[P = [0 0;0.5 1;1 0.5;0.5 0;0.25 0.5;.75 .75];\nplot(P(:,1),P(:,2),'*-') \naxis equal]]></w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>If F = 1/3, then the first few positions of the particle are:</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"code\"/></w:pPr><w:r><w:t><![CDATA[P = [0 0;1/3 1;1 2/3;2/3 0;1/9 1/3;5/9 8/9];]]></w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>Write a function that takes two inputs and returns two outputs. The first input argument, N, is the Nth location of the particle. The second input argument is the fraction F, where F is on (0,1). The first returned argument is a two element row vector that has the coordinates of the Nth location of the particle as [x,y]. The second returned argument is the limiting value of the particles location as N--&gt;infinity. For example,</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"code\"/></w:pPr><w:r><w:t><![CDATA[[P,L] = my_spiral(3,1/2)]]></w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>returns</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"code\"/></w:pPr><w:r><w:t><![CDATA[P = [1 0.5]]]></w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>and</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"code\"/></w:pPr><w:r><w:t><![CDATA[L = [?x? ?y?]]]></w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>Where both values of L are doubles accurate to 100*eps.</w:t></w:r></w:p></w:body></w:document>"},{"partUri":"/matlab/output.xml","contentType":"text/xml","content":"<?xml version=\"1.0\" encoding=\"UTF-8\" standalone=\"no\" ?><embeddedOutputs><metaData><evaluationState>manual</evaluationState><layoutState>code</layoutState><outputStatus>ready</outputStatus></metaData><outputArray type=\"array\"/><regionArray type=\"array\"/></embeddedOutputs>"}]}

Imagine a square with corners at the origin and at (1,1). An interesting spiral can be created by tracing the path of a moving point that starts at the origin and intersects the leftmost opposite edge at a point (F,1), where F is the fraction of the length of the intersected edge. The point then continues on in this way, always heading for the leftmost opposite side of its path and intersecting the segment at a point which is fraction F of its length. For example, if F = 1/2, the first few points are:<pre class="language-matlab">P = [0 0;0.5 1;1 0.5;0.5 0;0.25 0.5;.75 .75];
plot(P(:,1),P(:,2),'*-') 
axis equal
</pre>If F = 1/3, then the first few positions of the particle are:<pre class="language-matlab">P = [0 0;1/3 1;1 2/3;2/3 0;1/9 1/3;5/9 8/9];
</pre>Write a function that takes two inputs and returns two outputs. The first input argument, N, is the Nth location of the particle. The second input argument is the fraction F, where F is on (0,1). The first returned argument is a two element row vector that has the coordinates of the Nth location of the particle as [x,y]. The second returned argument is the limiting value of the particles location as N-->infinity. For example,<pre class="language-matlab">[P,L] = my_spiral(3,1/2)
</pre>returns<pre class="language-matlab">P = [1 0.5] 
</pre>and<pre class="language-matlab">L = [?x? ?y?] 
</pre>Where both values of L are doubles accurate to 100*eps.

Imagine a square with corners at the origin and at (1,1).  An interesting spiral can be created by tracing the path of a moving point that starts at the origin and intersects the leftmost opposite edge at a point (F,1), where F is the fraction of the length of the intersected edge.  The point then continues on in this way, always heading for the leftmost opposite side of its path and intersecting the segment at a point which is fraction F of its length.  For example, if F = 1/2, the first few points are:

P = [0 0;0.5 1;1 0.5;0.5 0;0.25 0.5;.75 .75];
  plot(P(:,1),P(:,2),'*-') 
  axis equal

If F = 1/3, then the first few positions of the particle are:

P = [0 0;1/3 1;1 2/3;2/3 0;1/9 1/3;5/9 8/9];

Write a function that takes two inputs and returns two outputs. The first input argument, N, is the Nth location of the particle. The second input argument is the fraction F, where F is on (0,1). The first returned argument is a two element row vector that has the coordinates of the Nth location of the particle as [x,y]. The second returned argument is the limiting value of the particles location as N-->infinity. For example,

[P,L] = my_spiral(3,1/2)

returns

P = [1 0.5]

and

L = [?x? ?y?]

Where both values of L are doubles accurate to 100*eps.

Solve

Solution Stats

96 Solutions
8 Solvers

Last Solution submitted on Jul 25, 2021

Last 200 Solutions

Problem Comments

2 Comments

2 Comments

Stephen on 29 Feb 2012

I can solve the simpler test cases, but how do you handle the eps stuff without too much error or crashing the program?

Leonard Gipson on 13 Dec 2021

good question

Solution Comments

Solution 54511

3 Comments

3 Comments

Stephen on 1 Mar 2012

any suggestions for handling my error?

Matt Fig on 2 Mar 2012

I see what you are trying to do, but there is an analytical solution for L. It seems you only have a problem with L, not p (as far as the test cases).

Stephen on 15 Mar 2012

I'll take another look

Solution 29518

1 Comment

1 Comment

Matt Fig on 5 Feb 2012

Octogenarian, Shmoctogenarian! Hope that one didn't keep you up too late, Roger. :-)

Solution 29382

1 Comment

1 Comment

@bmtran (Bryant Tran) on 5 Feb 2012

I feel like there should be an analytical solution for L, but it's too late at night to do any non-computational numerical analysis of this sequence.

Problem Recent Solvers8

Problem Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!