«Сандық әдістер»

	0.025
2	0.2 0.25 0.25 0.3	-0.94978 -0.93650 -0.93569 -0.92164	0.026553 0.028176 0.028138 0.030236	0.02824	0.023
3	0.3 0.35 0.35 0.4	-0.92154 -0.90642 -0.90514 -0.88881	0.030231 0.032790 0.032732 0.035743	0.03284	0.023
4	0.4 0.45 0.45 0.5	-0.88870 -0.87083 -0.86910 -0.84957	0.035745 0.039209 0.039134 0.04307	0.03925	0.022
5	0.5	-0.84945

Программасы:

Program Runge_Kutta;

Const a=0;b=1; n=5;

Var x,y,k1,k2,k3,k4:array[0..n] of real;

I:integer;

H:real;

Function f(x,y:real):real;

f:=y-(2*x)/y;

End;

Begin

Readln(x[0],y[0]);

x[n]:=b;

k1[0]:=h;

for i:=0 to n do

Begin

K2[i]:=h*f(x[i]+h/2,y[i]+k1[i]/2);

K3i]:=h*f(x[i]+h/2,y[i]+k2[i]/2);

K4[i]:=h*f(x[i]+h,y[i]+k3[i]);

Y1[i]:=1/6*(k1[i]+2*k2[i]+2*k3[i]+k4[i]);

y[i+1]:=y[i]+y1[i];

Writeln(x[i]:5:4,’ ‘,y[i]:5:4);

X[i+1]:=x[i]+h/2;

End;

End.
6-тақырып: «Қарапайым дифференциалдық теңдеулер үшін шектік есептер» (4сағ)

1.Қуалау және ақырлы-айырымдық, вариациялық әдістерді қолданып төмендегі қарапайым дифференциалдық теңдеулер үшін шектік есепті шешу:

• 
  , y(0)=1; y(1)=e^-1+1=1.367; e=10^-2
  
• 
  , y(0)=1; y(1)=e^-1+1=1.367; e=10^-2
  
• 
  , y(0)=1; y(1)=e^-1+1=1.367; e=10^-2
  
• 
  , y(0)= y(1)=0, a=1+0,4k, k=0,1,2.,
  
• 
  b=2,5+0,5n, n=0,1,2,3,4,5.
  
• 
  
  
  жүктеу/скачать 34,89 Mb.
  
  Достарыңызбен бөлісу: