Can Someone help me solve what is wrong with this code
1 view (last 30 days)
Show older comments
%Use symbolic processing with the variable x
syms x
syms fp(x)
syms Increasing
syms Decreasing
syms Minimum
syms Maximum
syms Inflection
syms CP
%Enconde the function f(x)
f(x) =(1+x^2)/(1-x^2)
%Set the distance to the left/right of the critical point
h=0.01
%Find the first derivative of the function
fp(x) = diff(1+(x^2))/(1-(x^2))
%Solve for the roots of fp
root1 = solve(fp,x)
%Consider the first critical point only
c=root1(1)
%Determine if the function is "Increasing" or "Decreasing" at the right of the critical point
if (fp(c+h)>0)
IoDR = "Increasing"
elseif (fp(c+h)<0)
IoDR = "Decreasing"
end
%Determine if the function is "Increasing" or "Decreasing" at the right of the critical point
if (fp(c-h)>0)
IoDL = "Decreasing"
elseif (fp(c-h)<0)
IoDL = "Increasing"
end
%Use first derivative Test to Determine if the critical point is a "Maximum" point or "Minimum" point.
if IoDL=="Increasing" & IoDR =="Decreasing"
CP= "Maximum"
elseif IoDL=="Decreasing" & IoDR =="Increasing"
CP= "Minimum"
end
%Find the second derivative of the function
fpp(x)= diff(f,2);
%Find the points of inflection of the function by equating the second derivative of the function to zero.
cc = diff(fp,x)
%Apply Second Derivative Test to check whether the critical point is a "Maximum" point, a "Minimum" point or "Point of Inflection".
if fpp(c) >0
CP2 = "Maximum"
elseif fpp(c)<0
CP2 = "Minimum"
else
CP2 = "Point of Inflection"
end
%GRAPH THE FUNCTION
clf();
g1= ezplot(f);
hold on
grid on
plot(c,f(c), 'r*')
title("Curve Tracing")
text(c+.5,f(c),["("+string(c)+","+string(f(c))+") "+CP ])
0 Comments
Answers (1)
VBBV
on 4 Oct 2022
Edited: VBBV
on 4 Oct 2022
%syms CP
Comment this line and try as CP is symbolic. Then you get know whether it's maximum or minimum
3 Comments
VBBV
on 6 Oct 2022
In the given code, the if-else condition for IoDL and IoDR variables are not satisfied since they are both increasing. so, output of CP will be unknown. In order to input CP, c and h needs to computed in such a way that condition is satisfied so that CP can be determined either as minimum or maximum strings.
See Also
Categories
Find more on Assumptions in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!