ranicu: n"); («% f% f», amp; t0, amp; tn); («Vvedite shag h: n"); («% f», amp; h ); (« nVvedite y1_0:»); («% f», amp; y1); (« nVvedite y2_0:»); («% f», amp; y2);
=(tn-t0)/h +1;/// kilkist vidrizkiv=t0;=(float *) calloc (N, sizeof (float)); (! X) exit (1);=(float **) calloc (K, sizeof (float)); (! Y) exit (1); (i=0; i lt; K; i ++)
{[i]=(float *) calloc (N, sizeof (float)); (! Y [i]) exit (1);
}=(float *) calloc (K, sizeof (float)); (! y) exit (1); [0]=y1; [1]=y2;
//************* (i=0; i lt; N; i ++)
{[0] [i]=y [0]; [1] [i]=y [1]; [i]=x; _Kutt (K, x, y, h); =x + h;
}
//******************** (« n # X Y1 Y2 n n"); (i=0; i lt; N; i ++)
{(«% d t% .4f t t% .4f t t% .4f n», i, X [i], Y [0] [i] , Y [1] [i]);
}
//********************
_ getch (); (X); (i=0; i lt; K; i ++) (Y [i]); (Y); (y); 0;
}
//**************** R_Kutt (int N, float x, float * y, float h)
{i=0; * k1=NULL, * k2=NULL, * k3=NULL, * k4=NULL, * y_temp=NULL, * F=NULL;
=(float *) calloc (N, sizeof (float));=(float *) calloc (N, sizeof (float));=(float *) calloc (N, sizeof (float ));=(float *) calloc (N, sizeof (float)); _ temp=(float *) calloc (N, sizeof (float));=(float *) calloc (N, sizeof (float)); ( ! k1 ||! k2 ||! k3 ||! k4 ||! y_temp ||! F) exit (1); (x, y, F); (i=0; i lt; N; i ++)
{[i]=h * F [i]; _ temp [i]=y [i] + 0.5 * k1 [i];
}=x + h * 0.5; (x, y_temp, F); (i=0; i lt; N; i ++)
{[i]=h * F [i]; _ temp [i]=y [i] + 0.5 * k2 [i];
} (x, y_temp, F); (i=0; i lt; N; i ++)
{[i]=h * F [i]; _ temp [i]=y [i] + k3 [i];
}=x + h * 0.5; (x, y_temp, F); (i=0; i lt; N; i ++)
{[i]=h * F [i]; [i]=y [i] + (k1 [i] + 2.0 * k2 [i] + 2.0 * k3 [i] + k4 [i])/6;
}
//*** Free memory *** (k1); (k2); (k3); (k4); (y_temp); (F);
} System (float x, float * y, float * F)
{[0]=y [1]; [1]=- 0.01 * exp (-x);
}
3. Тестування програми
Результат виконан програми зображено на малюнку 2.
Малюнок 2 - Результат виконан програми
Диференціальний Рівняння алгоритм
4. Виконання програми в середовіщі MathCad
Результат виконан програми в середовіщі MathCad зображено на малюнку 3.
Малюнок 3 - Результат виконан програми в середовіщі MathCad
Висновки : Вівче алгоритми набліженіх методів розв язання систем лінійніх алгебраїчніх рівнянь на ЕОМ.
