Uebung 12with(plots):
#with(linalg):
#with(stats):1)a)rec1a:=proc(n,xx)
local f0, f1, i:
#f0:=xx:
#f1:=xx/2;
if n=0 then return xx:
elif n=1 then return xx/2:
else return (3*rec1a(n-1,xx)+rec1a(n-2,xx)/2):
end if:
return xx;
end proc:
rec1a(3,2);f17_7:=rec1a(17,7);
evalf(%);b)fib:=proc(n)
option remember:
local f0, f1, i:
#f0:=xx:
#f1:=xx/2;
if n=1 then return 1:
elif n=2 then return 2:
else return (fib(n-1)+fib(n-2)):
end if:
return xx;
end proc:fib23:=fib(23);
fib29:=fib(29);
with(combinat, fibonacci);
kontrolle23:=seq(fibonacci(i), i=1..24);
kontrolle29:=fibonacci(30);2)a)restart:
with(linalg);
with(plots):f:=x->x*ln(x);
a:=1;
b:=4;#hilbert-Matrix (kubisch):
m4:=matrix([[b-a, (b^2-a^2)/2, (b^3-a^3)/3, (b^4-a^4)/4], [(b^2-a^2)/2, (b^3-a^3)/3, (b^4-a^4)/4, (b^5-a^5)/5], [(b^3-a^3)/3, (b^4-a^4)/4, (b^5-a^5)/5, (b^6-a^6)/6],[(b^4-a^4)/4, (b^5-a^5)/5, (b^6-a^6)/6, (b^7-a^7)/7] ]);#L\366sungsvektor:
v4:=vector([int(f(x),x=a..b), int(f(x)*x,x=a..b), int(f(x)*x^2,x=a..b), int(f(x)*x^3,x=a..b)]);loesung4:=linsolve(m4,v4);f_a:=loesung4[1]+loesung4[2]*x+loesung4[3]*x^2+loesung4[4]*x^3;
f_a:=unapply(f_a,x);plot1:=plot(f(x), x=-0..5, color=black):
plot2:=plot(f_a, -0..5, color=red):
display(plot1,plot2);b)a:='a';b:='b';
with(CurveFitting);
with(stats[fit]):x_i:=[1.0,2.0,3.0,4.0];
f_x:=map(f,x_i);ff := fit[leastsquare[[x,y], y=a*x^3+b*x^2+c*x+d, {a,b,c,d}]]([x_i,f_x]);
ff_curve:=sort(LeastSquares(x_i, f_x, x, curve= a*x^3+b*x^2+c*x+d ));f_b :=unapply(rhs(ff), x);plot1:=plot(f, -0..5,color=black):
plot2:=plot(f_b, -0..5, color=red):
display(plot1,plot2);c)c:=abs(f(x)-f_a(x)):
d:=abs(f(x)-f_b(x)):
c:=unapply(c,x);
d:=unapply(d,x);extremal:=extrema(d(x),y);
extrem_x:=fsolve(diff(d(x),x)=0, x=1..4);
extrem_y:=d(extrem_x);plot([c,d],1..4,color=[red,blue]);3)a)Length := proc(ll) return nops(ll); end proc:
Last := proc(ll) return op(-1, ll); end proc:
Append := proc(ll, x) return [op(ll), x]; end proc:
Member := proc(ll, p) return op(p, ll); end proc:
Take := proc(ll, n)
if n > 0 then return [ op(1..n, ll) ]
elif n < 0 then return [ op(n..-1, ll) ]
end if;
end proc:lagrPol := proc(xl, fl)
local i, j, l, ll, pol;
ll := []:
pol := 0:
i := 1:
while i <= Length(xl) do
l := 1;
j := 1;
while j <= Length(xl) do
if j <> i then
l := l*(x-Member(xl, j))/(Member(xl, i) - Member(xl, j));
end if;
j := j+1;
end do;
ll := Append(ll, l);
pol := pol + l*Member(fl, i);
i := i+1;
end do:
return( sort( expand( pol ) ) );
end proc:l:=lagrPol(x_i,f_x);poly:=unapply(l,x);plot(poly,1..4);b)nst4:=k->( cos( (2*k-1)*Pi/8) );
ltr:=x->0.5*( (4-1)*x +4+1);procA:restart:
with(stats):
x_i:=[1940,1950,1960,1970,1980,1990,2000];
f_x:=[132165,151326,179323,203302,226542,249633,274882];ff := fit[leastsquare[[x,y], y=a*x^6+ b*x^5 +c*x^4 +d*x^3+ e*x^2 +f*x +g, {a,b,c,d,e,f,g}]]([x_i,f_x]);
ff1 := fit[leastsquare[[x,y], y=a*x+b, {a,b}]]([x_i,f_x]);
evalf[5](ff);
func:=unapply(ff,x);
smartplot(func);
bew_2005:=evalf(func(2005));
#ff_curve:=sort(LeastSquares(x_i, f_x, x, curve= a*x^3+b*x^2+c*x+d ));procA:=proc()
local ls, i:
ls:=[1]:
i:=100:
while nops(ls) <101 do
ls:=[op(ls),(op(-1,ls)+2)];
end do:
return ls;
nops(ls);
end proc:procA();nops(ls);procB:=proc()
local ls, n, f, m, i:
ls:=[]:
f:=x->1-abs(cos(x)):
n:=0:
while n <10 do
ls:=[op(ls),sqrt(n)];
n:=n+1:
end do:
m:=map(f,ls):
return m;
end proc:procB();Sum('i^(3/2)','i'=10..50);
evalf(sum('i^(3/2)','i'=10..50));