|
|||
begin. whilet<=tk do. whilet<=tk do. while(y> 0) do1)кофе program cofee_text; const r=0. 1; T0=22; dt=0. 001; tk=25; var t, Temp: real; procedure Euler; begin whilet< =tk do begin Temp: =Temp-r*(Temp-T0)*dt; t: =t+dt; writeln(t: 2: 1, ' ', Temp: 2: 2); end; end; begin t: =0; Temp: =80; Euler; end.
program Coffe_text; uses GraphABC; const r=0. 395; T0=22; dt=0. 001; tk=15; var Sx, Sy, t, Temp: real; Hu, Hv: Integer; procedure out ( var x, y: real; c: system. Drawing. Color); begin var u, v: integer; u: = Round(Hu+Sx*x); v: = Round(Hv-Sy*y); PutPixel (u, v, c); end;
procedure Euler; begin whilet< =tk do begin Temp: =Temp-r*(Temp-T0)*dt; out (t, Temp, ClRed); t: =t+dt; end; end;
begin t: =0; Temp: =80; Sx: =Window. Width/tk; Sy: =Window. Height/100; Hu: =0; Hv: =Window. Height; Euler; end.
2)одномерное в гравитации program Dvigenie_dvum; uses GraphABC; var n, h, x, y, sk, k, t, dt, ax, ay, m, q, Hu, Hv, rx, ry, Sx, Sy, vx, vy: real; procedure Init; begin h: = 50; n: =12; sk: =5; q: =9. 8; k: =0; m: =0. 5; t: = 0; y: = 0. 1; dt: = 0. 00001; vx: =15; vy: =15; rx: = h; ry: = h; Sx: = Window. Width / rx; Sy: = Window. Height / ry; Hu: = 9; Hv: = Window. Height - 3; Line(0, round(Hv), Window. Width, round(Hv)); Line(round(Hu), Window. Height, round(Hu), 0); //Pen. Color: = clgreen; end; procedure OutRes(x, y: real; c: color); var v, u: word; begin u: = round(x * Sx + Hu); v: = round(-y * Sy + Hv); PutPixel(u, v, c); end; procedure Euler; begin while(y> 0) do begin sk: =sqrt(sqr(vx))+sqr(vy); ax: = (-k)/m *vx*sk+n; // +- n ветер ay: = -q - k/m * vy* sk; vx: = vx + ax * dt; vy: = vy + ay * dt; x: = x + vx * dt; y: = y + vy * dt; t: = t + dt; OutRes(x, y, clred); end; end; begin Init; OutRes(t, y, clred); Euler; end.
program lab21; uses graphABC; var a, dt, v, y, t, rx, ry, Sy, Sx, Hu, Hv, g, k2, m: real;
procedure OutRes(x, y: real; c: color); var u, v: word; begin u: =round(x*Sx + Hu); v: =round(-y*Sy + Hv); PutPixel(u, v, c); end;
procedure Init; begin v: =0; y: =50; t: =0; m: =20; k2: =1. 2; g: =9. 81; dt: =0. 001; rx: =10; ry: =100; Hu: =3; Hv: =Window. Height-3; Sx: =Window. Width/rx; Sy: =Window. Height/ry; Line(round(Hu), 0, round(Hu), Window. Height, clblack); Line(round(Hu), round(Hv), Window. Width, round(Hv), clbrown); end;
procedure Euler; begin
|
|||
|