Problem 282. Circle and Quadratic

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 the quadratic curve with equation</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>y=y(x)=ax^2+bx+c</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>On the concave side of this curve there is a circle of radius R. The circle is as close the the extremum of the quadratic as possible without resulting in the curves crossing each other. Write a function which takes as inputs a,b,c, and R and returns the coordinates of the center of the circle.</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>For example, if</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"code\"/></w:pPr><w:r><w:t><![CDATA[a=1; b=0; c=10; R=pi;]]></w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>then the function returns</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"code\"/></w:pPr><w:r><w:t><![CDATA[T = circ_puzz(a,b,c,R)\nT =\n     0   20.1196044010894]]></w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>This can be visualized as follows:</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"code\"/></w:pPr><w:r><w:t><![CDATA[ P = @(x) a*x.^2 + b*x + c;  % Quadratic\n C = @(x) real(-sqrt(R^2-(x-T(1)).^2) + T(2));  % Lower half circle\n x = linspace(-R,R,10000);  % Range of plotted data.\n plot(x,C(x),'r',x,-C(x)+2*T(2),'r',x,P(x),T(1),T(2),'*k')\n ylim([0,30])\n axis equal]]></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>"}]}

Solve

Solution Stats

183 Solutions
14 Solvers

Last Solution submitted on May 30, 2023

Last 200 Solutions

Problem Comments

3 Comments

3 Comments

Gunther S on 14 Feb 2012

the visualisation has an error, it should be

C = @(x) real(-sqrt(R^2-(x-T(1)).^2) + T(2)); % Lower half circle

Matt Fig on 16 Feb 2012

Thanks Gunther. The function does produce the correct graphic, but only when T(1)=0 as in my example.

Rafael S.T. Vieira on 9 Nov 2022

The third test is a pain. Why such a small circle for high precision? It's already a single-point tangent circle. It is no fun to look for the specific method the author had in mind when we already have the answer. IMHO, problems should be open; they accept all kinds of solutions (leaving people to explore clever ideas). Or they should be closed, and a specific method is requested (so people don't wander to dead ends). Precision should be a requirement only for closed ones.

Solution Comments

Show comments

Problem Recent Solvers14

Problem Tags

Community Treasure Hunt

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

Start Hunting!

Problem 282. Circle and Quadratic

Solution Stats

Last 200 Solutions

Problem Comments

Solution Comments

Problem Recent Solvers14

Suggested Problems

More from this Author6

Problem Tags

Community Treasure Hunt

Problem 282. Circle and Quadratic

Solution Stats

Last 200 Solutions

Problem Comments

Solution Comments

Problem Recent Solvers14

Suggested Problems

More from this Author6

Problem Tags

Community Treasure Hunt

Players