Question: Is the for..do loop in my procedure causing a problem?


I trying to get a proc to work. In the 1st half of the document, I derive four sections for a pedal curve and plot them fine. In the second I try doing it with a procedure called Pedal. Am having problems. I think it might be something to do with the for.. do loop. 

restart

NULLNULL

Typesetting:-mrow(Typesetting:-mi("Trunc", italic = "true", font_style_name = "2D Input", mathvariant = "italic"), Typesetting:-mo("≔", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.2777778em", rspace = "0.2777778em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo("proc", bold = "true", font_style_name = "2D Input", mathvariant = "bold", fontweight = "bold", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mfenced(Typesetting:-mrow(Typesetting:-mi("F", italic = "true", font_style_name = "2D Input", mathvariant = "italic"), Typesetting:-mo("∷", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.2777778em", rspace = "0.2777778em"), Typesetting:-mfenced(Typesetting:-mrow(Typesetting:-mi("procedure", italic = "true", font_style_name = "2D Input", mathvariant = "italic"), Typesetting:-mo(",", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "true", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.3333333em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mi("`+`", italic = "true", font_style_name = "2D Input", mathvariant = "italic")), font_style_name = "2D Input", mathvariant = "normal", open = "{", close = "}"), Typesetting:-mo(",", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "true", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.3333333em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mi("odr", italic = "true", font_style_name = "2D Input", mathvariant = "italic"), Typesetting:-mo("∷", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.2777778em", rspace = "0.2777778em"), Typesetting:-mi("posint", italic = "true", font_style_name = "2D Input", mathvariant = "italic"), Typesetting:-mo("≔", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.2777778em", rspace = "0.2777778em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mn("2", font_style_name = "2D Input", mathvariant = "normal"), Typesetting:-mo(",", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "true", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.3333333em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mi("v", italic = "true", font_style_name = "2D Input", mathvariant = "italic"), Typesetting:-mo("∷", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.2777778em", rspace = "0.2777778em"), Typesetting:-mi("list", italic = "true", font_style_name = "2D Input", mathvariant = "italic"), Typesetting:-mfenced(Typesetting:-mrow(Typesetting:-mi("name", italic = "true", font_style_name = "2D Input", mathvariant = "italic")), font_style_name = "2D Input", mathvariant = "normal"), Typesetting:-mo("≔", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.2777778em", rspace = "0.2777778em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mfenced(Typesetting:-mrow(Typesetting:-mi("x", italic = "true", font_style_name = "2D Input", mathvariant = "italic"), Typesetting:-mo(",", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "true", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.3333333em"), Typesetting:-mi("y", italic = "true", font_style_name = "2D Input", mathvariant = "italic")), font_style_name = "2D Input", mathvariant = "normal", open = "[", close = "]")), font_style_name = "2D Input", mathvariant = "normal"), Typesetting:-mi(""), Typesetting:-mspace(height = "0.0ex", width = "0.0em", depth = "0.0ex", linebreak = "newline"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo("description", bold = "true", font_style_name = "2D Input", mathvariant = "bold", fontweight = "bold", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-ms(" Truncates an algebraic equation to required degree"), Typesetting:-mo(";", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "true", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.2777778em"), Typesetting:-mspace(height = "0.0ex", width = "0.0em", depth = "0.0ex", linebreak = "newline"), Typesetting:-mspace(height = "0.0ex", width = "0.0em", depth = "0.0ex", linebreak = "auto"), Typesetting:-mo("local", bold = "true", font_style_name = "2D Input", mathvariant = "bold", fontweight = "bold", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mi("f", italic = "true", font_style_name = "2D Input", mathvariant = "italic"), Typesetting:-mo("≔", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.2777778em", rspace = "0.2777778em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mi("`if`", italic = "true", font_style_name = "2D Input", mathvariant = "italic"), Typesetting:-mfenced(Typesetting:-mrow(Typesetting:-mi("F", italic = "true", font_style_name = "2D Input", mathvariant = "italic"), Typesetting:-mo("∷", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.2777778em", rspace = "0.2777778em"), Typesetting:-mi("procedure", italic = "true", font_style_name = "2D Input", mathvariant = "italic"), Typesetting:-mo(",", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "true", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.3333333em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mi("F", italic = "true", font_style_name = "2D Input", mathvariant = "italic"), Typesetting:-mfenced(Typesetting:-mrow(Typesetting:-mi("v", italic = "true", font_style_name = "2D Input", mathvariant = "italic"), Typesetting:-mfenced(Typesetting:-mrow(Typesetting:-mi("")), font_style_name = "2D Input", mathvariant = "normal", open = "[", close = "]")), font_style_name = "2D Input", mathvariant = "normal"), Typesetting:-mo(",", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "true", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.3333333em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mi("F", italic = "true", font_style_name = "2D Input", mathvariant = "italic")), font_style_name = "2D Input", mathvariant = "normal"), Typesetting:-mo(";", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "true", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.2777778em"), Typesetting:-mi("print", italic = "true", font_style_name = "2D Input", mathvariant = "italic"), Typesetting:-mfenced(Typesetting:-mrow(Typesetting:-mi("F", italic = "true", font_style_name = "2D Input", mathvariant = "italic")), font_style_name = "2D Input", mathvariant = "normal"), Typesetting:-mo(";", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "true", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.2777778em"), Typesetting:-mspace(height = "0.0ex", width = "0.0em", depth = "0.0ex", linebreak = "newline"), Typesetting:-mspace(height = "0.0ex", width = "0.0em", depth = "0.0ex", linebreak = "auto"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo("if", bold = "true", font_style_name = "2D Input", mathvariant = "bold", fontweight = "bold", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo("not", bold = "true", font_style_name = "2D Input", mathvariant = "bold", fontweight = "bold", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mi("f", italic = "true", font_style_name = "2D Input", mathvariant = "italic"), Typesetting:-mo("∷", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.2777778em", rspace = "0.2777778em"), Typesetting:-mi("`+`", italic = "true", font_style_name = "2D Input", mathvariant = "italic"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo("then", bold = "true", font_style_name = "2D Input", mathvariant = "bold", fontweight = "bold", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mspace(height = "0.0ex", width = "0.0em", depth = "0.0ex", linebreak = "newline"), Typesetting:-mspace(height = "0.0ex", width = "0.0em", depth = "0.0ex", linebreak = "auto"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo("error", bold = "true", font_style_name = "2D Input", mathvariant = "bold", fontweight = "bold", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-ms("Can't truncate 1-term expression"), Typesetting:-mspace(height = "0.0ex", width = "0.0em", depth = "0.0ex", linebreak = "newline"), Typesetting:-mspace(height = "0.0ex", width = "0.0em", depth = "0.0ex", linebreak = "auto"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo("else", bold = "true", font_style_name = "2D Input", mathvariant = "bold", fontweight = "bold", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mspace(height = "0.0ex", width = "0.0em", depth = "0.0ex", linebreak = "newline"), Typesetting:-mspace(height = "0.0ex", width = "0.0em", depth = "0.0ex", linebreak = "auto"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mi("select", italic = "true", font_style_name = "2D Input", mathvariant = "italic"), Typesetting:-mfenced(Typesetting:-mrow(Typesetting:-mi("q", italic = "true", font_style_name = "2D Input", mathvariant = "italic"), Typesetting:-mo("→", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.2777778em", rspace = "0.2777778em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mi("degree", italic = "true", font_style_name = "2D Input", mathvariant = "italic"), Typesetting:-mfenced(Typesetting:-mrow(Typesetting:-mi("q", italic = "true", font_style_name = "2D Input", mathvariant = "italic"), Typesetting:-mo(",", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "true", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.3333333em"), Typesetting:-mi("v", italic = "true", font_style_name = "2D Input", mathvariant = "italic")), font_style_name = "2D Input", mathvariant = "normal"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo("≤", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.2777778em", rspace = "0.2777778em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mi("odr", italic = "true", font_style_name = "2D Input", mathvariant = "italic"), Typesetting:-mo(",", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "true", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.3333333em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mi("f", italic = "true", font_style_name = "2D Input", mathvariant = "italic")), font_style_name = "2D Input", mathvariant = "normal"), Typesetting:-mi(""), Typesetting:-mspace(height = "0.0ex", width = "0.0em", depth = "0.0ex", linebreak = "newline"), Typesetting:-mspace(height = "0.0ex", width = "0.0em", depth = "0.0ex", linebreak = "auto"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo("fi", bold = "true", font_style_name = "2D Input", mathvariant = "bold", fontweight = "bold", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mspace(height = "0.0ex", width = "0.0em", depth = "0.0ex", linebreak = "newline"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo("end", bold = "true", font_style_name = "2D Input", mathvariant = "bold", fontweight = "bold", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo(" ", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo("proc", bold = "true", font_style_name = "2D Input", mathvariant = "bold", fontweight = "bold", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.0em", rspace = "0.0em"), Typesetting:-mo(":", font_style_name = "2D Input", mathvariant = "normal", fence = "false", separator = "false", stretchy = "false", symmetric = "false", largeop = "false", movablelimits = "false", accent = "false", lspace = "0.2777778em", rspace = "0.2777778em"))

NULLNULLNULL

NULL

a := 3; b := -1; c := -1

C := (x^2+y^2+12*x+9)^2-4*(2*x+3)^3

(x^2+y^2+12*x+9)^2-4*(2*x+3)^3

(1)

p21 := plots:-implicitplot(C, x = -3 .. 3, y = -5 .. 5, colour = "Salmon", gridrefine = 2, size = [300, 300])

 

sol1 := [solve(C, y)]

[(-x^2-12*x-9+2*(8*x^3+36*x^2+54*x+27)^(1/2))^(1/2), -(-x^2-12*x-9+2*(8*x^3+36*x^2+54*x+27)^(1/2))^(1/2), (-x^2-12*x-9-2*(8*x^3+36*x^2+54*x+27)^(1/2))^(1/2), -(-x^2-12*x-9-2*(8*x^3+36*x^2+54*x+27)^(1/2))^(1/2)]

(2)

NULL

f1 := expand(eval(C, [x = X+r, y = Y+s]))

X^4+4*X^3*r+2*X^2*Y^2+4*X^2*Y*s+6*X^2*r^2+2*X^2*s^2+4*X*Y^2*r+8*X*Y*r*s+4*X*r^3+4*X*r*s^2+Y^4+4*Y^3*s+2*Y^2*r^2+6*Y^2*s^2+4*Y*r^2*s+4*Y*s^3+r^4+2*r^2*s^2+s^4-8*X^3-24*X^2*r+24*X*Y^2+48*X*Y*s-24*X*r^2+24*X*s^2+24*Y^2*r+48*Y*r*s-8*r^3+24*r*s^2+18*X^2+36*X*r+18*Y^2+36*Y*s+18*r^2+18*s^2-27

(3)

ltg1 := simplify(expand(eval(Trunc(f1, 1, [Y, X]), [X = x-r, Y = y-s])))

-3*r^4+(4*x+16)*r^3+(-6*s^2+4*s*y-24*x-18)*r^2+((4*x-48)*s^2+48*s*y+36*x)*r-27-3*s^4+4*s^3*y+(24*x-18)*s^2+36*s*y

(4)

NULLThis*section*works*out*the*pedal*curve

This*section*works*out*the*pedal*curve

(5)

xp1 := 1; yp1 := 2; lprp1 := -coeff(ltg1, y)*x+coeff(ltg1, x)*y+K1; for i to nops(sol1) do s := eval(sol1[i], x = r); print("s  ", s); K1 := solve(eval(lprp1, [x = xp1, y = yp1]), K1); print("K1  ", K1); sol2 := solve([ltg1, lprp1], [x, y]); Pedal || i := [rhs(sol2[1, 1]), rhs(sol2[1, 2])]; p9 || i := plot([op(Pedal || i), r = -20 .. 20], colour = "MediumSeaGreen"); unassign('K1', 'sol2') end do

"K1  ", -4*(-r^2-12*r-9-2*(8*r^3+36*r^2+54*r+27)^(1/2))^(3/2)-4*r^2*(-r^2-12*r-9-2*(8*r^3+36*r^2+54*r+27)^(1/2))^(1/2)+16*(8*r^3+36*r^2+54*r+27)^(1/2)*r-48*r*(-r^2-12*r-9-2*(8*r^3+36*r^2+54*r+27)^(1/2))^(1/2)+192*r^2+96*(8*r^3+36*r^2+54*r+27)^(1/2)-36*(-r^2-12*r-9-2*(8*r^3+36*r^2+54*r+27)^(1/2))^(1/2)+576*r+432

(6)

seq(simplify(Pedal || i), i = 1 .. nops(sol1))

[((((2*r+3)^3)^(1/2)+6*r+9)*(-r^2-12*r-9+2*((2*r+3)^3)^(1/2))^(1/2)-2*(r-3)*(r^2+(5/2)*r-((2*r+3)^3)^(1/2)+3/2))/(4*r^2-6*r-18), ((-2*r^2-3*r-((2*r+3)^3)^(1/2))*(-r^2-12*r-9+2*((2*r+3)^3)^(1/2))^(1/2)+6*(r-3)*(r-(1/6)*((2*r+3)^3)^(1/2)+3/2))/(4*r^2-6*r-18)], [((-((2*r+3)^3)^(1/2)-6*r-9)*(-r^2-12*r-9+2*((2*r+3)^3)^(1/2))^(1/2)-2*(r-3)*(r^2+(5/2)*r-((2*r+3)^3)^(1/2)+3/2))/(4*r^2-6*r-18), ((2*r^2+((2*r+3)^3)^(1/2)+3*r)*(-r^2-12*r-9+2*((2*r+3)^3)^(1/2))^(1/2)+6*(r-3)*(r-(1/6)*((2*r+3)^3)^(1/2)+3/2))/(4*r^2-6*r-18)], [((6*r-((2*r+3)^3)^(1/2)+9)*(-r^2-2*((2*r+3)^3)^(1/2)-12*r-9)^(1/2)-2*(r-3)*(r^2+(5/2)*r+((2*r+3)^3)^(1/2)+3/2))/(4*r^2-6*r-18), ((-2*r^2+((2*r+3)^3)^(1/2)-3*r)*(-r^2-2*((2*r+3)^3)^(1/2)-12*r-9)^(1/2)+6*(r-3)*(r+(1/6)*((2*r+3)^3)^(1/2)+3/2))/(4*r^2-6*r-18)], [((((2*r+3)^3)^(1/2)-6*r-9)*(-r^2-2*((2*r+3)^3)^(1/2)-12*r-9)^(1/2)-2*(r-3)*(r^2+(5/2)*r+((2*r+3)^3)^(1/2)+3/2))/(4*r^2-6*r-18), ((2*r^2-((2*r+3)^3)^(1/2)+3*r)*(-r^2-2*((2*r+3)^3)^(1/2)-12*r-9)^(1/2)+6*(r-3)*(r+(1/6)*((2*r+3)^3)^(1/2)+3/2))/(4*r^2-6*r-18)]

(7)

plots:-display(p21, seq(p9 || i, i = 1 .. nops(sol1)), size = [600, 600], scaling = constrained, caption = " Curve and it's Pedal curve")

 

This I can't get to work. All 4 sections of the pedal curve are the same and should not be.

Parse:-ConvertTo1D, "`%1` is not a module or member", RonanRoutines

Pedal2 := Pedal(C, xp1, yp1)

"K||i", -8*r^3+(-4*(-r^2-12*r-9-2*(8*r^3+36*r^2+54*r+27)^(1/2))^(1/2)+48)*r^2-8*(-(-r^2-12*r-9-2*(8*r^3+36*r^2+54*r+27)^(1/2))^(1/2)-3)^2*r-4*(-r^2-12*r-9-2*(8*r^3+36*r^2+54*r+27)^(1/2))^(3/2)+48*r^2+576*r+432+96*(8*r^3+36*r^2+54*r+27)^(1/2)-36*(-r^2-12*r-9-2*(8*r^3+36*r^2+54*r+27)^(1/2))^(1/2)

(8)

seq(simplify(Pedal2[i]), i = 1 .. nops(sol1))

[((4*r^2+6*((2*r+3)^3)^(1/2)+30*r+36)*(-r^2-2*((2*r+3)^3)^(1/2)-12*r-9)^(1/2)+(-r^2-34*r-69)*((2*r+3)^3)^(1/2)-26*r^3-243*r^2-540*r-351)/((2*r+30)*((2*r+3)^3)^(1/2)+36*r^2+162*r+162), (((r+9)*((2*r+3)^3)^(1/2)+14*r^2+51*r+45)*(-r^2-2*((2*r+3)^3)^(1/2)-12*r-9)^(1/2)+(12*r+72)*((2*r+3)^3)^(1/2)+4*r^3+126*r^2+432*r+378)/((2*r+30)*((2*r+3)^3)^(1/2)+36*r^2+162*r+162)], [((4*r^2+6*((2*r+3)^3)^(1/2)+30*r+36)*(-r^2-2*((2*r+3)^3)^(1/2)-12*r-9)^(1/2)+(-r^2-34*r-69)*((2*r+3)^3)^(1/2)-26*r^3-243*r^2-540*r-351)/((2*r+30)*((2*r+3)^3)^(1/2)+36*r^2+162*r+162), (((r+9)*((2*r+3)^3)^(1/2)+14*r^2+51*r+45)*(-r^2-2*((2*r+3)^3)^(1/2)-12*r-9)^(1/2)+(12*r+72)*((2*r+3)^3)^(1/2)+4*r^3+126*r^2+432*r+378)/((2*r+30)*((2*r+3)^3)^(1/2)+36*r^2+162*r+162)], [((4*r^2+6*((2*r+3)^3)^(1/2)+30*r+36)*(-r^2-2*((2*r+3)^3)^(1/2)-12*r-9)^(1/2)+(-r^2-34*r-69)*((2*r+3)^3)^(1/2)-26*r^3-243*r^2-540*r-351)/((2*r+30)*((2*r+3)^3)^(1/2)+36*r^2+162*r+162), (((r+9)*((2*r+3)^3)^(1/2)+14*r^2+51*r+45)*(-r^2-2*((2*r+3)^3)^(1/2)-12*r-9)^(1/2)+(12*r+72)*((2*r+3)^3)^(1/2)+4*r^3+126*r^2+432*r+378)/((2*r+30)*((2*r+3)^3)^(1/2)+36*r^2+162*r+162)], [((4*r^2+6*((2*r+3)^3)^(1/2)+30*r+36)*(-r^2-2*((2*r+3)^3)^(1/2)-12*r-9)^(1/2)+(-r^2-34*r-69)*((2*r+3)^3)^(1/2)-26*r^3-243*r^2-540*r-351)/((2*r+30)*((2*r+3)^3)^(1/2)+36*r^2+162*r+162), (((r+9)*((2*r+3)^3)^(1/2)+14*r^2+51*r+45)*(-r^2-2*((2*r+3)^3)^(1/2)-12*r-9)^(1/2)+(12*r+72)*((2*r+3)^3)^(1/2)+4*r^3+126*r^2+432*r+378)/((2*r+30)*((2*r+3)^3)^(1/2)+36*r^2+162*r+162)]

(9)

for i to nops(sol1) do p9 || i := plot([op(Pedal2[i]), r = -20 .. 20], colour = "MediumSeaGreen") end do

 

``

NULL

NULL

``


 

Download Proc_for_Pedal_not_working.mw

Please Wait...