Паскаль, вычисление методом чебышева интегралла

[Anya Lis]

посмотрите сейчас, не нравиться t2. посоветуйте еще куда засунуть e (точность для вычисления интегралла)

program kyrsov;
uses crt;
begin
writeln ('vvod a,b,n,e:');
readln (a,b,n,e);
s:=0;
var
t5:array[1..5]of real:=(0.832498,0.374541,0,-0.374541,-0.832498);
t2:array[1..2]of real:=(0.577350,-0.577350);
t3:array[1..3]of real:=(0.707107,0,-0.707107);
t4:array[1..4]of real:=(0.794654,0.187592,-0.187592,-0.794654);
t6:array[1..6]of real:=(0.866247,0.422519,0.266635,-0.266635,-0.422519,-0.866247);
t7:array[1..7]of real:=(0.883862,0.529657,0.321912,0,-0.321912,-0.529657,-0.883862);
x:array[1..n]of real;
y:array[1..n]of real;
s,b,a:real;
i,n:integer;
function f(x:real):real;
var i:integer;
a,b,s:real;
begin
f:=sin(0.4*x)-3*ln(x);
end;

for i:=1 to n do
begin
x[i]:=(b+a)/2+((b-a)/2)*tn[i];
y[i]:=f(i);
s:=s+y[i];
end;
s:=s*(b-a)/2;
writeln('s=',s:10:7);
Readln;
end.

Anya Lis, вы хоть внимание обращайте на то, что ваши посты правят модераторы, оформляя код тегами. Для правильного оформления кода используйте кнопку # вверху

с нулем я решила вопрос, скобку поставить забыла

[newerow1989]

разве не будет ошибки деления на ноль?! Ведь t[3] равно нулю, если я не ошибаюсь!

Там умножить на 0:

*t[i];

[newerow1989]

Ошибок - целая куча. Вот как правильно писать на Паскале, с учетом синтаксиса:

program kyrsov;
uses crt;

function f(x:real):real;
begin
   f:=sin(0.4*x)-3*ln(x);
end;

const n=10;//!!!!!!!!!!!!!!!!! см. ниже, где употребляется
var i:integer;
    a,b,s,e:real;
    t5:array[1..5]of real=(0.832498,0.374541,0,-0.374541,-0.832498);
    t2:array[1..2]of real=(0.577350,-0.577350);
    t3:array[1..3]of real=(0.707107,0,-0.707107);
    t4:array[1..4]of real=(0.794654,0.187592,-0.187592,-0.794654);
    t6:array[1..6]of real=(0.866247,0.422519,0.266635,-0.266635,-0.422519,-0.866247);
    t7:array[1..7]of real=(0.883862,0.529657,0.321912,0,-0.321912,-0.529657,-0.883862);
    x:array[1..n]of real;//!!!!!!!!!!!!!!!!!!! n - const
    y:array[1..n]of real;

begin
   writeln ('vvod a,b,e:');//// здесь пришлось убрать
   readln (a,b,e);//// тут тоже
   s:=0;
   for i:=1 to n do
   begin
      //x[i]:=(b+a)/2+((b-a)/2)*tn[i];// tn не знает (пришлось отключить эту строку - компилятор ругается), рекомендую по очереди писать t2, t3, t4, t5, t6, t7
      y[i]:=f(i);
      s:=s+y[i];
   end;
   s:=s*(b-a)/2;
   writeln('s=',s:10:7);
   Readln;
end.

[Serge_Bliznykov]

Сообщение от newerow1989

Там умножить на 0:

да, точно, Вы правы! Это меня переклинило что-то...