Хелпикс

Главная

Контакты

Случайная статья





begin. whilet<=tk do. whilet<=tk do. while(y> 0) do



1)кофе

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



  

© helpiks.su При использовании или копировании материалов прямая ссылка на сайт обязательна.