second verion maxima and minima points of a function
Subject: second verion maxima and minima points of a function
From: Luigi Marino
Date: Sun, 8 Feb 2009 10:38:26 +0100
Hi
My block (simplest function how said Mr. Fateman)
have the purpose to test really capability of CAS Maxima.
Solve is very hard, Maple and other CAS have any limititation.
It is a new version that found also inflection point only
as a zero of second derivative.
load(descriptive);
stud_funz(f):=block([p,fs,fp],
p:solve(diff(f,x,1),x),
fs:diff(f,x,2),
fp:ev(fs,x:p),
for i:1 thru length(float(p)) do
if floor(rhs(fp[i]))#0 and rhs(fp[i])>0 then print("minimo",[p[i],rhs(ev(f,x:p[i]))])
else if floor(rhs(fp[i]))#0 then print("massimo",[p[i],rhs(ev(f,x:p[i]))])
else print("flesso",[p[i],rhs(ev(f,x:p[i]))]))$
Example:
f:x^3/(1-x^2);
minimo[x=-sqrt(3),(3*sqrt(3))/2]
massimo[x=sqrt(3),-(3*sqrt(3))/2]
flesso[x=0,0]
Best regards.
Luigi Marino