Метод Монте-Карло

[Sasha124]

Тема: Пошук кореня функції f(x1,x2) методом Монте-Карло
Завдання:
1. Задана функція f(x1,x2).
2. Задані обмеження в вигляді круга радіуса R с координатами центра x1c, x2c.
3. Згенерувати дві послідовності координати {x1}, {x2} довжиною 300.
4. Відібрати з них лише ті точки, які належать заданому кругу.
5. Вся площа круга повинна бути покрита генерованими точками.
6. Визначити значення кореня функції, за критерієм |f(x1,x2)| <E.
7. Зберігати у масиві запису данні пошуку с полями номер точки пошуку і зазначення функції в цій точці.
8. В графічному режимі ввести графік пошуку по горизонтальній вісі номер точки, а по вертикальній значення функції.


Помогите завтра нада сдать курсак ,а прогу которую я написал препод не принял

[NiCola999]

так выкладывайте то, что сами написали. Иначе тему закроют

[Sasha124]

program monte_karlo;
uses crt,graph;
type
result = record
i: integer;
res: real;
end;
var
op: char;
x1c,x2c,rc,i,n,m,gd,gm: integer;
x1,x2,curr: array [1..300] of integer;
res: array [1..300] of result;
k1,k2,e,y: real;
begin
clrscr;
writeln('Vvod funkcii y=k1*x1 operacia k2*x2:');
write('Operacia (+,-,:,*) = ');
readln(op);
write('Koeficient K1: ');
readln(k1);
write('Koeficient K2: ');
readln(k2);
writeln;
writeln('Vvod funkcii y=',k1:1:3,'x1',op,k2:1:3,'x2');
writeln('Koordunatu centra kryga:');
write('Enter X1c: ');
readln(x1c);
write('Enter X2c: ');
readln(x2c);
write('Enter radiys Rc: ');
readln(rc);
write('E: ');
readln(e);
randomize;
for i:=1 to 300 do begin
x1[i]:=random(640);
x2[i]:=random(480);
end;
n:=0;
for i:=1 to 300 do
if ((x1[i] >= x1c-rc) and (x1[i] <= x1c+rc)) and
((x2[i] >= x2c-rc) and (x2[i] <= x2c+rc))
then begin
inc(n);
curr[n]:=i;
end;
m:=0;
for i:=1 to n do begin
case op of
'+': y:=k1*x1[curr[i]]+k2*x2[curr[i]];
'-': y:=k1*x1[curr[i]]-k2*x2[curr[i]];
':': y:=k1*x1[curr[i]]/k2*x2[curr[i]];
'*': y:=k1*x1[curr[i]]*k2*x2[curr[i]];
end;
if y < e then begin
inc(m);
res[m].i:=curr[i];
res[m].res:=y;
end;
end;
readkey;
gd:=Detect;
InitGraph(gd,gm,'');
Line(0, (GetMaxY div 2), GetMaxX, (GetMaxY div 2));
Line(20, 0, 20, GetMaxY);
SetColor(green);
for i:=1 to m-1 do begin
MoveTo(res[i].i+20, (GetMaxY div 2)-round(res[i].res));
LineTo(res[i+1].i+20, (GetMaxY div 2)-round(res[i+1].res));
end;
readkey;
CloseGraph;
end.

[NiCola999]

вставьте код с помощью специального тега

[Sasha124]

Я не знаю как.

[NiCola999]

кнопка #

[Alt]

При наборе сообщения либо вручную вводите тег [CODE][ /CODE], либо пользуетесь кнопкой