Árpád apánk

Blog. ennyi. semmi több. megtalalsz itt jó pár infót rólam, az életemről, meg néhány érdekes dolgot is, csak épp fejlesszem ki. addig is türelem kispajtás :)

2008/05/20

harmadik

> f:=x->sin(3*x)-cos(4*x);

>

>

f := proc (x) options operator, arrow; sin(3*x)-cos(4*x) end proc

> fsolve(f(x),x=-2..-1);

>

-1.570796327

> fsolve(f(x),x=-1..0);

-.6731984258

> fsolve(f(x),x=0..1);

.2243994753

> fsolve(f(x),x=1..2);

1.121997376

> plot(f(x),x=-2..2);

[Plot]

> felezo:=proc(f,bal,jobb,eps)
local a,b,kozep;

a:=bal;

b:=jobb;

while b-a>2*eps do

kozep:=(a+b)/2;

if evalf(f(a)*f(kozep))<0>
b:=kozep

else

a:=kozep

end if;

end do;

RETURN (kozep)

end:

> felezo(f,-2.,-1.,0.0001);

-1.000122070

> felezo(f,-1.,0.,0.0001);

-.6732177733

> felezo(f,0.,1.,0.0001);

.2244873047

> felezo(f,1.,2.,0.0001);

1.121948242

> g:=x->3*sqrt(x+2)-3;

g := proc (x) options operator, arrow; 3*sqrt(x+2)-3 end proc

> plot({g(x),x},x=-2..5);

[Plot]

> fsolve(g(x)=x,x=-2..5);

4.854101966

> plot(g(x)-x,x=-2..5);

[Plot]

> c:=tg(alpha)=D(g)(4);

c := tg(alpha) = 6^(1/2)/4

> v:=evalf(c);

v := tg(alpha) = .6123724358

> d:=arctan(rhs(v));

d := .5494672448

> evalf(convert(d,degrees));

31.48215410*degrees

>