Здраствуйте. Выложу сначала код:
Код Pascal
program Runge_Kutt;
uses crt;
const
h = 0.1;
n = 50;
var
file1: text;
i: integer;
xi, yi, dy, dx: real;
k1, k2, k3, k4, k11, k12, k13, k14: real;
function FotXY(x, y: real): real;
begin
FotXY := (70-x*0.031-y*1.2)/0.01;
end;
function FotYX(x, y: real): real;
begin
FotYX := (y*1.2-148)/13;
end;
begin
clrscr;
assign(file1,'file1.txt');
rewrite(file1);
writeln(file1,' k | xk | yk |');
write(file1,i:3, ' ', xi:2:3, ' ', yi:2:3);
xi := 0;
yi := 1;
write('Расчёт методом Рунге-Кутты 4го порядка точности...');
for i := 0 to n do
begin
k1 := h * fotxy(xi, yi);
k2 := h * fotxy(xi + h/2, yi + k1/2);
k3 := h * fotxy(xi + h/2, yi + k2/2);
k4 := h * fotxy(xi + h, yi + k3);
k11 := h * fotyx(xi, yi);
k12 := h * fotyx(xi + h/2, yi + k1/2);
k13 := h * fotyx(xi + h/2, yi + k2/2);
k14 := h * fotyx(xi + h, yi + k3);
dy := (k1 + 2*k2 + 2*k3 + k4) / 6;
dx := (k11 + 2*k12 + 2*k13 + k14) / 6;
yi := yi + dy;
xi := xi + dx;
write(file1,i:3, ' ', xi:12:5, ' ', yi:12:5);
end;writeln('Расчёт окончен!');
close(file1);
end.
Вообщем не могу выдать нормальную структуру в файле file1, получается всё сдвинутым а для того что бы Mathcad прочитал файл с этого файла надо что б всё было ровно в несколько столбцов. То есть сейчас в file1 так:
i | xi | yi |
0 0.000 0.000 0 27.44933 -36463.13950 1 17893.92421 -23263669.00857 2 11398994.12827 -14818700261.42470 3 7261016365.99329 -9439324818154.14400 4 4625175658649.36896 -6012730610702397.44000 5 2946178439478830.08000 -3.830033407381392E18 6 1.87667843079908576E18 -2.43968287479018656E21 7 1.1954204420997432E21 -1.55404715741473088E24 8 7.6146771334792128E23 -9.89908398515320064E26 9 4.85045309625826368E26 -6.30559138939728512E29 10 3.08967732007459776E29 -4.0165820221016992E32 11 1.96808540413418432E32 -2.55851198468041312E35 12 1.2536455288711168E35 -1.62973979859823328E38 13 7.98556357746042624E37 -1.0381236543110352E41 14 5.08670307364199616E40 -6.61271647515171584E43 15 3.2401655698329024E43 -4.21221682014014208E46 16 2.06394451729097984E46 -2.68312887850895712E49 17 1.3147065724407344E49 -1.7091191850017856E52 18 8.37451470782527616E51 -1.088687320217152E55 19 5.3344600279423456E54 -6.93480063650660736E57 20 3.39798361845666176E57 -4.41738036027644928E60 21 2.1644726196877248E60 -2.8138154606252464E63 22 1.37874170314738224E63 -1.79236488613314016E66 23 8.78241039737429376E65 -1.14171377974062768E69 24 5.59428442701302144E68 -7.27257248194497024E71 25 3.56348847688523904E71 -4.63253675690580992E74 26 2.26989712278072608E74 -2.95086736603332736E77 27 1.44589577921459616E77 -1.87966521775350176E80 28 9.21017337468355456E79 -1.19732298764132784E83 29 5.86676403729510784E82 -7.62679610812670976E85 30 3.73705454491314368E85 -4.85817272994358144E88 31 2.38045651450724928E88 -3.09459462916788928E91 32 1.51632071444429568E91 -1.97121766787941216E94 33 9.65877131147172096E93 -1.25564074128984704E97 34 6.1525152534435264E96 -7.99827282840367488E99 35 3.91907445814528448E99 -5.0947987058675264E102 36 2.496400898785312E102 -3.24532238524426752E105 37 1.59017582187133024E105 -2.06722934353383136E108 38 1.0129219011636088E108 -1.31679896524229744E111 39 6.45218449271529728E110 -8.3878429855250752E113 40 4.1099599762046688E113 -5.34294997238825792E116 41 2.61799256129076608E116 -3.40339160576891072E119 42 1.66762817415632992E119 -2.1679174392570144E122 43 1.06225807069094544E122 -1.38093600967581696E125 44 6.76644965727395584E124 -8.79638785263393792E127 45 4.31014291420182976E127 -5.60318788935995072E130 46 2.74550656279154336E130 -3.56915986987420672E133 47 1.74885298153212512E133 -2.27350972843704224E136 48 1.11399724643310784E136 -1.44819696335990944E139 49 7.09602166771815552E138 -9.22483162685503488E141 50 4.520076119573168E141 -5.876101158670656E144
А надо так:
*****i***************Xi********** Yi
***0.00100***0.99019***0.02785
***0.00200***0.98075***0.05540
***0.00300***0.97168***0.08266
(вместо звёзд пробелы)
RPI.su - самая большая русскоязычная база вопросов и ответов. Наш проект был реализован как продолжение популярного сервиса otvety.google.ru, который был закрыт и удален 30 апреля 2015 года. Мы решили воскресить полезный сервис Ответы Гугл, чтобы любой человек смог публично узнать ответ на свой вопрос у интернет сообщества.
Все вопросы, добавленные на сайт ответов Google, мы скопировали и сохранили здесь. Имена старых пользователей также отображены в том виде, в котором они существовали ранее. Только нужно заново пройти регистрацию, чтобы иметь возможность задавать вопросы, или отвечать другим.
Чтобы связаться с нами по любому вопросу О САЙТЕ (реклама, сотрудничество, отзыв о сервисе), пишите на почту [email protected]. Только все общие вопросы размещайте на сайте, на них ответ по почте не предоставляется.