Возможно, ваша ваша проблема уже имеет решение на нашем форуме
04.05.2011, 13:05
#2 (permalink)
Извините — ничего не понял! Уравнение вполне себе легко решается аналитически, и выглядит это так:
1. Записываем в виде y’ — y = 1 + x. 2. Решаем однородное уравнение y’ — y = 0. Его решением, очевидно, является C*exp(x). 3. Частное решение ищем методом вариации постоянной, т.е. в виде C(x)*exp(x). Подставляя в исходное уравнение и интегрируя, находим C(x)=C1-exp(-x)*(x+2). 4. Таким образом, окончательно y=C1*exp(x)-x-2.
И при чем тут Паскаль, программирование, точность и всё такое прочее? Кроме того, полагаю, что Ваше предположение о том, что у=1 имеет какое-то отношение к таинственному «шагу», есть полный абсурд. Скорее всего, у=1 — это значение искомой функции при некотором значении аргумента х (например, при х=0), которое Вы пропустили, и которое требуется для определения константы С1.
2.6.Программа на языке Turbo Pascal
for i:=0 to n-1 do
Делись добром 😉
Похожие главы из других работ:
3. ПРОГРАММА И ПРИМЕР РЕАЛИЗАЦИИ АЛГОРИТМА БПФ С ПРОРЕЖИВАНИЕМ ПО ВРЕМЕНИ
Ниже при водится программа вычисления БПФ с прореживанием по времени по способу с замещением и рассматриваются при меры реализации этой программы. Программа 1 — быстрое преобразование Фурье с основанием два и прореживанием по времени.
4. Программа для решения дифференциального уравнения в Visual Basic
Dim y(9) As Single Dim YE(9) As Single Dim YR(9) As Single Dim YT(9) As Single Dim l(9) As Single Private x0 As Single Private Function fun(a As Single, b As Single) As Single f = (b + 2) / (a + 1) fun = f End Function Private Sub Command1_Click() x0 = Val(Text1.Text) xk = Val(Text2.Text) y0 = Val(Text3.Text) h = Val(Text4.Text) N = (xk — x0) / h MSFlexGrid1.Rows = N + 2 MSFlexGrid1.TextMatrix(0, 0) = «x» MSFlexGrid1.
3.1 Програма складена на мові програмування TURBO PASCAL 7.0
program MonteCarlo; uses crt; const k=100; Var a,b,c,d,ng,vg,x,y,s,integral : real; n,i,j : integer; integr:array[1..k]of real; Function f(x.
Программа и её описание. Результаты вычислений
Текст основной программы приведён ниже. Программа 1. #include «stdafx.h» #include #include //необходимые заголовочные файлы #include #include using namespace std; double U(double t,double x,double u) < return u; >double fi(double x) < return pow(x,2)/2; //функция.
program kursovaya; uses crt; const sizemat=10; type mattype=array[1..sizemat,1..sizemat] of double; mattype1=array[1..sizemat] of double; <Процедура для вывода матрицы на экран>procedure writemat (var a:mattype; n,m:byte); var i,j:byte; begin writeln; for i:=1 to n do begin for j:=1 to m do write(a[i,j]:7:3.
2.1. Программа эксперимента
Практическая часть выпускной квалификационной работы представляет собой 3 вида эксперимента: констатирующий, формирующий и контрольный. Констатирующий эксперимент включает в себя контрольную работу.
3.6 Численное решение задачи с использованием Turbo Pascal 7.0
Текст программы: program rungekutt2_v3; var x0,xn,y0,z0,h,x,y,z,k1,k2,k3,k4,l1,l2,l3,l4,c:real; n,i:integer; inp,ou:text; function f1(z:real):real; begin f1:=z end; function f2(x,y,z:real):real; begin f2:=x-x*z-sin(y)+1 end; begin assign(inp,inp2.txt); reset(inp); read(inp,x0,xn,n,y0,z0); close(inp); h:=(xn-x0)/n; x:=x0; y:=y0; z:=z0; assign(ou,ou2.txt); rewrite(ou); writeln(ou, x, ,y, ,z); writeln(ou.
§3. Программа элективного курса
Пояснительная записка Сюжетные задачи — это наиболее древний вид школьных задач. Они всегда широко использовались, и будут использоваться в обучении математике.
Численное решение краевой задачи для обыкновенного дифференциального уравнения второго порядка на Паскале (Pascal) — Лабораторная работа
1. Постановка задачи 3
2. Схема алгоритма. 4
3. Текст программы на Паскале 5
4. Результаты расчёта 8
5. Список литературы 9
Численное решение краевой задачи для обыкновенного дифференциального уравнения второго порядка.
3. Текст программы на Паскале
yx,xy,l,v,p,ff,ay,by,x:array [0.10] of real;
y,a,b:array[0.10,0.1] of real;
writeln(‘введите наивысший порядок производной не больше трех ‘);
if (n=1) and (l[1]=0) or (n=2) and (l[2]=0) or (n=3) and (l[3]=0) then begin
writeln(‘деление на ноль’);
writeln(‘введите коэффициент при x’);
writeln(‘задайте начальные условия y(x)= ‘);
for i:=0 to n-1 do
if n=3 then begin
for i:=0 to o-1 do begin
for i:=0 to o-1 do begin
if n=2 then begin
for i:=0 to o-1 do begin
for i:=0 to o-1 do begin
if n=1 then begin
for i:=0 to o-1 do begin
for i:=0 to o-1 do begin
4. Результаты расчёта
Общее решение однородного уравнения имеет вид
Тогда частное решение при заданных начальных условиях можно записать в виде:
1. Калиткин Н. Н. Численные методы. – М.: Наука, 1978. – 512 с.
2. Самарский А. А. Теория разностных схем. – М.: Наука, 1989. – 616 с.
3. Самарский А. А., Николаев Е. С. Методы решения сеточных уравнений. – М.: Наука, 1978. – 592 с.
4. Рапаков Г. Г., Ржеуцкая С. Ю. Программирование на языке Pascal. . – СПб.: БХВ-Петербург, 2005. – 480 с.
Готовые решение задачи на языке Паскаль
К работе прилагается все исходники (Pascal) и отчет (Word)
Не подошла эта работа? Узнайте стоимость написания работы по Вашему заданию.
