Производная функции, исправить код

@qwiryll · Регистрация: 14.02.2014

Author24 — интернет-сервис помощи студентам

Дoброгo времени суток. Я вoт попробовал написать программу, берущую производную от прoстенькой функции ( y(x) =x^3 - x^2 +1 для примера) и выбирающей макcимум из погрешностей в точках, но что-то не вышло, видимо, где-то ошибаюсь. Прoшу помощи в пoиске oшибки.

Если чуть подробней описывать задание, то я тут ввожу функцию на отрезке [a,b], нахожу шаг по формуле h=(b-a)/n, нахожу производную в точках отрезка по формулам прaвой (в точку a), лeвой (в точке b), цeнтральнoй(остальные ) рaзностных пр-х и сравниваю с точным значением. Там, где погрешность e самая большая (maxe) , выбираю максимум погрешности (maxpoint) и и указываю в какой точке узла.

Компилятор выдал мне:

prog.pas: In main program:
prog.pas:31: error: syntax error before `else'

Вот моя печатанина:

Pascal

program dff_1;
const a= 1.0; b=3.0;
var n: integer; i,x,df,h,e,maxe,maxpoint: real;
function f(x:real):real;
begin
f:= sqr(x)*x - sqr(x) +1;
end;
function ftrue(x:real):real;
begin
ftrue:= 3* sqr(x) - 2*x; end;
procedure max(var e,i:real);
begin
if e > maxe then 
maxe:= e;
maxpoint:= i; 
if i>=b then
writeln('Максимальная погрешность:', maxe,'; в точке:',i); end;
begin
writeln('Введите, пожалуйста, число n!');
readln(n);
h:=(b-a)/n;
e:=0; maxe:=0; maxpoint:=a;
i:=a;
while i<= b do begin
    if i=a then begin
    df:=(f(a+h)-f(a))/h;
    e:=abs(df- ftrue(a));
    writeln('f`(a)=', df, 'e=', e);
    max(e,i);
    i:=i+h;
    else 
    if i=b then begin
    df:=(f(b)-f(b-h))/h;
    e:=abs(df - ftrue(a));
    writeln('f`(b)=', df, 'e=', e);
    max (e,i);
    else begin df:= (f(i+h)- f(i-h))/2h;
    e:= abs (df - ftrue(i));
    writeln('f`(', i, ')=', df, 'e=', e);
    max(e,i);
    i:=i+h; end;
end;
end.

@Puporev · 14.02.2014, 18:11

Если бы код писал по людски и сам бы ошибки нашел

Pascal

program dff_1;
const a= 1.0; b=3.0;
var n: integer; i,x,df,h,e,maxe,maxpoint: real;
function f(x:real):real;
begin
f:= sqr(x)*x - sqr(x) +1;
end;
function ftrue(x:real):real;
begin
ftrue:= 3* sqr(x) - 2*x; end;
procedure max(var e,i:real);
begin
if e > maxe then
maxe:= e;
maxpoint:= i;
if i>=b then
writeln('Максимальная погрешность:', maxe,'; в точке:',i); end;
begin
writeln('Введите, пожалуйста, число n!');
readln(n);
h:=(b-a)/n;
e:=0; maxe:=0; maxpoint:=a;
i:=a;
while i<= b do
 begin
  if i=a then
   begin
    df:=(f(a+h)-f(a))/h;
    e:=abs(df- ftrue(a));
    writeln('f`(a)=', df, 'e=', e);
    max(e,i);
    i:=i+h;
   end
  else if i=b then
   begin
    df:=(f(b)-f(b-h))/h;
    e:=abs(df - ftrue(a));
    writeln('f`(b)=', df, 'e=', e);
    max (e,i);
   end
  else
   begin
    df:= (f(i+h)- f(i-h))/(2*h);
    e:= abs (df - ftrue(i));
    writeln('f`(', i, ')=', df, 'e=', e);
    max(e,i);
    i:=i+h;
   end;
end;
end.

@qwiryll · 14.02.2014, 19:20 **[ТС]**

Если бы код писал по людски и сам бы ошибки нашел

Спасибо, учту замечание

За исправления тоже благодарен.

Сам я тоже вот, две ошибки в процедуре подправил, из-за которых на экран ерунда выводилась, вместо нужного.

Pascal

program dff_1;
const a =1.0; b=3.0;
var n: integer; i,x,df,h,e,maxe,maxpoint: real;
function f(x:real):real;
begin
f:= sqr(x)*x - sqr(x) +1;
end;
function ftrue(x:real):real;
begin
ftrue:= 3* sqr(x) - 2*x; end;
procedure max(var e,i:real);
begin
if e > maxe then begin
maxe:= e;
maxpoint:= i; end;
if i>=b then
writeln('Максимальная погрешность:', maxe,'; в точке:',maxpoint); end;
begin
writeln('Введите, пожалуйста, число n!');
readln(n);
h:=(b-a)/n;
e:=0; maxe:=0; maxpoint:=a;
i:=a;
while i<= b do
 begin
  if i=a then
   begin
    df:=(f(a+h)-f(a))/h;
    e:=abs(df- ftrue(a));
    writeln('f`(',a:6:2,')=', df:6:2, ',   e=', e:6:2);
    max(e,i);
    i:=i+h;
   end
  else if i=b then
   begin
    df:=(f(b)-f(b-h))/h;
    e:=abs(df - ftrue(a));
    writeln('f`(b)=', df:6:2, ',   e=', e:6:2);
    max (e,i);
   end
  else
   begin
    df:= (f(i+h)- f(i-h))/(2*h);
    e:= abs (df - ftrue(i));
    writeln('f`(', i:6:2, ')=', df:6:2, ',   e=', e:6:2);
    max(e,i);
    i:=i+h;
   end;
end;
writeln('Максимальная погрешность:', maxe:6:2,'; в точке:', maxpoint:6:2);
end.

	14.02.2014, 19:20