Неявная схема и метод прогонки. Ур-ие теплопроводности.

[Matveykinpk]

$ Переведите пожалуйста код с Fortran на C/C++.\\ Программа описывает метод прогонки для смешанной краевой задачи для уравнения теплопроводности($f(x,t) = 0,\varphi = 15\sin5x+5,\mu_1 = 5,\mu_2 = 0,a=0.2,l = \frac{\pi}{2}$):\\ \begin{equation}\label{1} \left\{ \begin{array}{llll} u_{t} &=& a^2 u_{xx} + f(x,t)&, 0 < x < l, 0 < t < T ; \\ u(x,0) &=& \varphi(x)&, 0 \leq x \leq l ; \\ u(0,t) &=& \mu_{1}(t), u_{x}(l,t) = \mu_{2}(t)&, 0 \leq t \leq T . \end{array} \right. \end{equation} $

PROGRAM HEAT
IMPLICIT DOUBLE PRECISION (A-H,O-Z)
INTEGER,PARAMETER::NJ=100
INTEGER,PARAMETER::NN=100
DIMENSION U(NJ,NN)
DIMENSION A(NJ),B(NJ)
OPEN(UNIT=1,FILE='DATA.DAT')
! ВВОД КОНСТАНТ
HH=0.5*3.14/(NJ-1)
TT=0.1/(NN-1)
AA=0.2
GG=0.2*0.2*TT/HH/HH
! НАЧАЛЬНЫЕ УСЛОВИЯ
DO J=1,NJ
X=HH*(J-1)
U(J,1)=15.0*SIN(5*X)+5.0
ENDDO
! ГРАНИЧНЫЕ УСЛОЫИЯ
A(2)=0.0
B(2)=5.0    ! MU
A(NJ)=1.0
B(NJ)=0.0
! РАСЧЕТ
DO I=2,NN
U(1,I)=5.0  !MU
U(NJ,I)=U(NJ,I-1)   !???????
! ПРЯМОЙ ХОД ПРОГОНКИ
DO J=3,NJ-1
A(J)=GG/(1.0+2.0*GG-GG*A(J-1))
B(J)=(U(J,I-1)+GG*B(J-1))/(1.0+2.0*GG-GG*A(J-1))
!PRINT*,I,J,A(J),B(J)
!PAUSE
ENDDO
! ОБРАТНЫЙ ХОД ПРОГОНКИ
DO J=1,NJ-2
U(NJ-J,I)=A(NJ+1-J)*U(NJ+1-J,I)+B(NJ+1-J)
ENDDO
ENDDO

! ВЫВОД ДАННЫХ
DO I=1,NN
TIME=TT*(I-1)
DO J=1,NJ
X=HH*(J-1)
WRITE(1,*)TIME,X,U(J,I)
ENDDO
ENDDO

CLOSE(101)
END PROGRAM

[mclaudt]

500$

[Matveykinpk]

Здраво мыслишь для нашего времени!=)

[Munin]

Вам не кажется, что вот эта просьба: "Переведите пожалуйста код с Fortran на C/C++." - сама по себе не здравая?

[mclaudt]

Matveykinpk в сообщении #461267 писал(а):

Здраво мыслишь для нашего времени!=)

Не надо ля-ля. Видно же, что раздолбай-студень ниасилил лабу, насерфил в интернете кусок чьего-то кода и хочет прокатиться на шее.

Если задача реальная, и видна твоя заинтересованность, то тут (и не только) помогут и обсудят всегда с удовольствием. А вот помогать неучам - дураков плодить с дипломами только.

[Matveykinpk]

Задание было таким: была смешанная краевая задача уравнения теплопроводности. Я уже писал, как выглядит. Так вот, краевую задачу решал сам (Тихонов,Самарский - уравнения мат.физики), разностную схему, а то есть неявную схему строил сам, метод прогонки для данного случая выводил сам, отчет вбивал в LaTeX сам!Проблема в программировании. Так что прошу не кидаться громкими словами про диплом и неуча и не хамить. Если не хотите помогать не надо, я что заставляю! Единственное что мне помогали, так это написание программы на Fortran на другом компьютере. К сожалению фортран мне не удалось нормально настроить на своем компьютере,все время при прогонке вылезает ошибка о том, что не хватает определенного файла. Видел данную проблему в интернете, но решения так и не нашел. Поэтому и попросил перевести на C/C++, так как сам в программировании не силен и никогда не любил!

Кому интересно выкладываю свой отчет!
http://www.fayloobmennik.net/743245

[Vince Diesel]

Есть разные компиляторы фортрана. Можно какой-нибудь другой попробовать. Есть еще транстляторы, которые переводят с фортрана на С.

[ewert]

Matveykinpk в сообщении #461302 писал(а):

Поэтому и попросил перевести на C/C++,

Видите ли, Вы либо знаете Си -- либо не знаете. Если знаете, то переведёте на него код с любого другого языка (который хоть чуть понимаете) мгновенно. Если не знаете -- то всё вообще бесполезно.

[Gortaur]

ewert
Фортран тоже надо знать, видимо (по-моему, это такой лысый усатый дядька с котом, гусеницей и воробьем).

[ИСН]

(Оффтоп)

Да-да. Кота звали Икс.

Ну, читать проще, чем писать.

[TOTAL]

Matveykinpk в сообщении #461302 писал(а):

Единственное что мне помогали, так это написание программы на Fortran на другом компьютере. К сожалению фортран мне не удалось нормально настроить на своем компьютере,все время при прогонке вылезает ошибка о том, что не хватает определенного файла. Видел данную проблему в интернете, но решения так и не нашел. Поэтому и попросил перевести на C/C++, так как сам в программировании не силен и никогда не любил!

От Вас требуют самостоятельно написать программу или разрешают подсунуть чужую?

U(NJ,I)=U(NJ,I-1) !???????

[Munin]

Matveykinpk в сообщении #461302 писал(а):

неявную схему строил сам, метод прогонки для данного случая выводил сам... Проблема в программировании.

Звучит как "Евгения Онегина сочинил сам, проблема буковки в столбик записать".

И почему все относятся с недоверием?..

[Matveykinpk]

тема закрыта. Думал поможете, сам все сделал уже.

[Gortaur]

Простите