这道pascal题目怎么做

这道pascal题目怎么做

function gcd(a,b:int64):int64;
begin
  if b=0 then gcd:=a
    else gcd:=gcd(b,a mod b);
end ;

procedure plus(p1,q1,p2,q2:int64;var p,q:int64);
var
p1_,q1_,temp:int64;
begin
p1_:=p1*q2+p2*q1;
q1_:=q1*q2;
temp:=gcd(p1_,q1_);
p:=p1_ div temp;
q:=q1_ div temp;
end;

var
p1,q1,p2,q2:int64;
p,q:int64;
i:longint;
begin
p1:=0; q1:=1;
for i:=1 to 40 do begin
p2:=1; q2:=i;
plus(p1,q1,p2,q2,p,q);
p1:=p; q1:=q;
writeln(p,'/',q);
end;
end.