criterioderivada.m

% PROGRAMA QUE NOS DICE EL TIPO DE CARACTERISTICAS DE PUNTO FIJO
% VERSION 1.0
% AUTOR:EDGAR BERNAL FLORES
% 01032009
%email:[e-mail]

function[eval tipo]=criterioderivada(funcion,puntosfijos,varsyms)

% ENTRADAS DEL PROGRAMA:
% 
%             FUNCION-> ES UN ARREGLO QUE CONTIENE LAS FUNCIONES DE DONDE SE OBTUBIERON LOS PUNTOS FIJO
%             
%             PUNTOSFIJOS-> VALORES EN DONDE HAY PUTNOS FIJOS
% SALIDA DEL PROGRAMA:
% 
%             TIPO -> ES un arreglo de numeros donde nos dice:
%                                 -1 -> es atractor
%                                 0-> es nuetro
%                                 1-> es repulsor
% PROCEDIMIENTO: SE APLICA EL CRITERIO : |f'(X)| EVALUADA EN EL PUNTO FIJO SI >1 REPULSOR
%                                                                             =1 NEUTRO
%                                                                             <1 ATRACTOR

% DEFINIMOS LAS VARIABLES SIMBOLICAS PARA PODER DERIVAR LAS FUNCIONES
 syms(varsyms);
 
%  REALIZAMOS LA DERIVADA DE LA FUNCION Y LAS EVALUAMOS EN LOS PUTNOS FIJOS
%  DEL SISTEMA
for i2=1:length(funcion(:,1))
    ecuacionpf=funcion(i2,:);
    for i3=1:length(puntosfijos(i2,:))
        eval(i2,i3)=subs(diff(ecuacionpf),x,double(puntosfijos(i2,i3)));
%         eval(i2,i3)=double(eval(i2,i3));
        if(abs(eval(i2,i3))<1)%atractor
           tipo(i3)=-1;
        elseif(abs(eval(i2,i3))>1)%repulsor
           tipo(i3)=1;
        else%neutro
           tipo(i3)=0;
        end
    end
end