1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45
| restart;
a := proc (n::nonnegint, r)
option remember;
if n = 0 then
1;
elif n = 1 then
0;
else
-(r+n-1)*(r-n+2)*a(n-2)/(n*(n-1))
end if;
end proc;
Y1 := proc (x, N)
local s, n, r;
s := 0.;
r:=4.5;
for n from 0 to N do
s := s+a(n,r)*x^n
end do
end proc;
with(plots);
p1 := plot(Y1(x, 0), x = -2 .. 2, y = -3 .. 3, color = red);
%;
p2 := plot(Y1(x, 2), x = -2 .. 2, y = -3 .. 3, color = green);
%;
p3 := plot(Y1(x, 4), x = -2 .. 2, y = -3 .. 3, color = blue);
%;
Error, (in a) invalid input: a uses a 2nd argument, r, which is missing
p4 := plot(Y1(x, 6), x = -2 .. 2, y = -3 .. 3, color = magenta);
%;
Error, (in a) invalid input: a uses a 2nd argument, r, which is missing
display(p1, p2, p3, p4);
%;
Error, (in plots:-display) expecting plot structure but received: p3
a(0, 4.5);
print(`output redirected...`); # input placeholder
1
a(2, 4.5);
print(`output redirected...`); # input placeholder
-12.37500000
a(4, 4.5);
%;
Error, (in a) invalid input: a uses a 2nd argument, r, which is missing |
Partager