Составить пользовательскую процедуру или функцию для вычисления n-ого члена последовательности: b0=cos21; b1=sin21; bk=2 bk-1+ bk-2, где k=2,3,4, ··· n — Pascal(Паскаль)

uses crt;
function Bn(n:integer):real;//функция
var b0,b1,b:real;
    i:integer;
begin
b0:=cos(21);//нулевой
b1:=sin(21);//первый
for i:=2 to n do//от 2 и дальше
 begin
  b:=2*b1+b0;//по формуле
  b0:=b1;//первый станет нулевым
  b1:=b;//второй-первым и все повторяем
 end;
Bn:=b;//значение функции
end;
var n:integer;
begin
clrscr;
write('n=');
readln(n);
write('b[',n,']=',Bn(n):0:3);
readln
end.