Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- function[vektor]=proba4(p)
- elso_ertek=0;
- masodik_ertek=0;
- gyokei=roots(p);
- ertek=roots(polyder(p));
- derivalt=polyder(p);
- masodik_derivalt=polyder(derivalt);
- x=ertek-1; %teszt értékek
- y=ertek+1;
- if polyval(derivalt,x)<0
- elso_ertek=-1;
- else
- elso_ertek=1;
- end
- if polyval(derivalt,y)<0
- masodik_ertek=-1;
- else
- masodik_ertek=1;
- end
- if elso_ertek==-1 && masodik_ertek==1
- loc_min=polyval(p,ertek);
- vektor(1)=loc_min;
- vektor(2)=0;
- elseif elso_ertek==1 && masodik_ertek==-1
- loc_max=polyval(p,ertek);
- vektor(1)=0;
- vektor(2)=loc_max;
- else
- vektor(1)=0;
- vektor(2)=0;
- end
- %itt az inflexiós pontot szeretném meghatározni.
- inf_pont_ertek=roots(polyder(derivalt));
- %vizsgálom a nála kisebb és nagyobb elemeket, hogy ismerjem a konvexitást,
- %mert ezekből határozható meg az inflexiós pont
- if polyval(masodik_derivalt,inf_pont_ertek-1)<0
- masodik_elso=-1;
- else
- masodik_elso=1;
- end
- if polyval(masodik_derivalt,inf_pont_ertek+1)<0
- masodik_masodik=-1;
- else
- masodik_masodik=1;
- end
- %ugye ott van inflexiós pont ahol - ból + ba vált vagy fordítva, ha nem
- %vált nincs inflexiós pont
- if masodik_elso==-1 && masodik_masodik==1 || masodik_elso==1 && masodik_masodik==-1
- vektor(3)=polyval(p,inf_pont_ertek);
- else
- vektor(3)=0;
- end
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement