построение графика функции Опции

[yar11]

program Grafik_funkcii;
uses crt,graph;
var
 {координаты графика функции}
 x,y,a,b:real;
 {Номера типа графического дрaйвера и номера графического режима}
 driver,mode:integer;
 {cx,cy - координаты центра системы координат на экране,
 mx,my - масштаб (количество пикселей на единицу),
 ex,ey - координаты точки графика в масштабе экрана}
 cx,cy,mx,my,ex,ey,n:integer;
 {Описание функции, график которой затем выводится на экран}
 function f(x:real):real;
   begin
            f:=sin(3*x)+cos(x)/exp(x)
    end;
     procedure DeleniyaX(a:real);
     var
     stroka:string;x,s:real;
     begin
     {Деления выставляются в 10-ти точках от -a до a c шагом s}
      x:=-a;s:=2*a/10;
      repeat
        ex:=trunc(cx+x*mx)-10;
        ey:=cy+5;
        str(x:3:1,stroka);
        outtextxy(ex,ey,stroka);
        x:=x+s;
      until x>a;
     end;
     procedure DeleniyaY(b:real);
      var
       stroka:string;y:real;s:integer;
        begin
        {Деления по оси Y выставляются в 10 точках от -b до b}
         s:=trunc((getmaxy-30)/10);ey:=getmaxy-30;
         repeat
           ex:=cx+5;
           y:=b-ey/my;
           str(y:3:1,stroka);
           if abs(y)>0.5 then outtextxy(ex,ey+5,stroka);
           ey:=ey-s;
         until ey<10;
        end;
        {Подпрограмма построения графика}
        procedure MyGrafik(a,b:real);
        {Процедура строит график от -a до a c шагом s}
         var p,n:integer;s,t:real;
             strx:string;
          begin
           {Положение начала координат на экране}
           cx:=trunc(getmaxx/2);cy:=trunc(getmaxy/2);
           {Цвета фона и текста}
           SetBkColor(yellow);SetColor(blue);
           line(10,cy,getmaxx-50,cy);{Ось X}
           {Стрелочка на оси X}
           line(getmaxx-60,cy-5,getmaxx-50,cy);
           line(getmaxx-60,cy+5,getmaxx-50,cy);
           line(cx,10,cx,getmaxy-20);  {Ось Y}
           {Рисование стрелочки на оси y}
           line(cx-5,20,cx,10);line(cx+5,20,cx,10);
           setcolor(magenta);
           {Подписи к осям}
           outtextxy(cx+10,10,'y');
           outtextxy(getmaxx-50,cy-10,'x');
           {Определение масштаба по осям X и Y исходя из того,
            что X от -a до a}
           mx:=trunc((getmaxx)/(2*a));
           {Y от -b до b}
           my:=trunc((getmaxy)/(2*b));
           {Вывод значений по оси X от -a до a}
           DeleniyaX(a);
           {Вывод значений по оси Y от -b до b}
           DeleniyaY(b);
           outtextxy(10,50,'График функции');
           {x меняется от -a до a c шагом s в n+1 точке}
           x:=-a;n:=500;s:=2*a/n;
           moveto(trunc(cx+x*mx),trunc(cy-f(x)*my));
           Repeat
            y:=f(x);
            ex:=trunc(cx+x*mx);
            ey:=trunc(cy-y*my);
            lineto(ex,ey);
            x:=x+s;
           Until ex>getmaxx-70;
          end;
 Begin
  clrscr;
  write('a=');readln(a);
  write('b=');readln(b);
  {Открытие графического режима}
  driver:=detect;
  initgraph(driver,mode,'d:\tp\bgi');
  {очистка экрана в графическом режиме}
  cleardevice;
  {выполнение процедуры вывода графика}
  MyGrafik(a,b);
  {задержка}
  repeat
  {до нажатия клавиши Enter}
  until keypressed;
 End.

это программа построения графика функции по точкам
почему-то она не работает при данной функции f:=sin(3*x)+cos(x)/exp(x)
при значениях b больше 19
пишет ошибку invalid floating point operation
может кто подскажет

Сообщение отредактировано: volvo - 5.11.2006 19:00

[volvo]

Другими словами:
F = 10 *
((-1)⁰ * 19! * (4*x²)¹⁰) / (0! * 20!) +
((-1)¹ * 18! * (4*x²)⁹) / (1! * 18!) +
((-1)² * 17! * (4*x²)⁸) / (2! * 16!) +
((-1)³ * 16! * (4*x²)⁷) / (3! * 14!) +
((-1)⁴ * 15! * (4*x²)⁶) / (4! * 12!) +
((-1)⁵ * 14! * (4*x²)⁵) / (5! * 10!) +
((-1)⁶ * 13! * (4*x²)⁴) / (6! * 8!) +
((-1)⁷ * 12! * (4*x²)³) / (7! * 6!) +
((-1)⁸ * 11! * (4*x²)²) / (8! * 4!) +
((-1)⁹ * 10! * (4*x²)¹) / (9! * 2!) +
((-1)¹⁰ * 9! * (4*x²)⁰) / (10! * 0!)
?

Тогда ты будешь смеяться, но тебе вообще не нужна НИ функция вычисления факториала, НИ функция возведения в степень... Все, что нужно знать - это чему равен первый член ряда, и как изменяется относительно предыдущего последующий... Вот так будет выглядеть твоя функция:

function f(x: real): real;
var
  s, item: real;
  p, k: integer;
begin
  item := 1/10; s := item;
  for k := 9 downto 0 do begin

    if k = 9 then p := 1 else p := k + 1;
    item := - item * ((19 - k) * p * 4 * sqr(x)) /
                     (pred(2 * (10 - k)) * (2 * (10 - k)));
    s := s + item;
  end;
  f := 10 * s;
end;

которая при х=1 равна 1'200'000 с хвостиком

Это ты преувеличила, при X = 1 значение функции равно 1, но вот растет оно очень быстро - уже F(1.1) становится равным 357.15, а уж при Х = 2 вообще 13 737 919 114.6, то есть больше 13 миллиардов...

Сообщение отредактировано: volvo - 5.07.2008 12:36