Модераторы: Poseidon
  

Поиск:

Ответ в темуСоздание новой темы Создание опроса
> [Java] График 
:(
    Опции темы
Merhaba
Дата 30.3.2014, 10:23 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Шустрый
*


Профиль
Группа: Участник
Сообщений: 93
Регистрация: 23.4.2011

Репутация: нет
Всего: нет



Добрый день! Помогите Пожалуйста написать код, чтобы строился график для метода Рунге-Кутта 4-ого порядка!

Код

package metodrungekutta;

import java.io.*;
import static java.lang.Math.*;

public class Main {

PrintWriter out_rk;

    public static void main(String[] args) throws Exception {
    new Main().run();
    }
     int n = 40;

    public void run() throws Exception {
        out_rk = new PrintWriter(new FileWriter("RK.txt"));
        solve();
        out_rk.flush();
        out_rk.close();
    }

    public void solve() throws Exception {
        RK();
    }

    public double F(double x) {
        return 0.125*exp(-cos(x+1));
    }

    public double f(double x, double y) {
        return F(x) - g(x) * y;
    }

    public double g(double x) {
        return -sin(x+1);
    }

    public double RK_4(double x, double y, double h) {
        double k1 = f(x, y);
        double k2 = f(x + h / 2, y + h*k1 / 2);
        double k3 = f(x + h / 2, y + h*k2 / 2);
        double k4 = f(x + h, y + h*k3);
        return y + h*(k1 + 2 * k2 + 2 * k3 + k4) / 6;
    }

public void RK(){
        double x0 = -1;
        double b = 2*PI-1;
        double y0 = 8;
        double y1 = y0;
                double x1 = x0;
        out_rk.printf("x = %.7f   y = %.7f\n", x0, y0);
        double h = (b - x0) / n;
        for (int i = 0; i < n; i++) {
            double y_nplus = RK_4(x1, y1, h);
            double x_nplus = x1 + h;
                        x1 = x_nplus;
                        y1 = y_nplus;
            out_rk.printf("x = %.7f   y = %.7f\n", x_nplus, y_nplus);
        }
    }
}


PM MAIL   Вверх
Mirkes
Дата 3.4.2014, 20:56 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Опытный
**


Профиль
Группа: Участник
Сообщений: 586
Регистрация: 18.8.2011
Где: Красноярск

Репутация: 4
Всего: 17



А в чем проблема?


--------------------
Mirkes
PM MAIL   Вверх
  
Ответ в темуСоздание новой темы Создание опроса
Правила форума "Центр помощи"

ВНИМАНИЕ! Прежде чем создавать темы, или писать сообщения в данный раздел, ознакомьтесь, пожалуйста, с Правилами форума и конкретно этого раздела.
Несоблюдение правил может повлечь за собой самые строгие меры от закрытия/удаления темы до бана пользователя!


  • Название темы должно отражать её суть! (Не следует добавлять туда слова "помогите", "срочно" и т.п.)
  • При создании темы, первым делом в квадратных скобках укажите область, из которой исходит вопрос (язык, дисциплина, диплом). Пример: [C++].
  • В названии темы не нужно указывать происхождение задачи (например "школьная задача", "задача из учебника" и т.п.), не нужно указывать ее сложность ("простая задача", "легкий вопрос" и т.п.). Все это можно писать в тексте самой задачи.
  • Если Вы ошиблись при вводе названия темы, отправьте письмо любому из модераторов раздела (через личные сообщения или report).
  • Для подсветки кода пользуйтесь тегами [code][/code] (выделяйте код и нажимаете на кнопку "Код"). Не забывайте выбирать при этом соответствующий язык.
  • Помните: один топик - один вопрос!
  • В данном разделе запрещено поднимать темы, т.е. при отсутствии ответов на Ваш вопрос добавлять новые ответы к теме, тем самым поднимая тему на верх списка.
  • Если вы хотите, чтобы вашу проблему решили при помощи определенного алгоритма, то не забудьте описать его!
  • Если вопрос решён, то воспользуйтесь ссылкой "Пометить как решённый", которая находится под кнопками создания темы или специальным флажком при ответе.

Более подробно с правилами данного раздела Вы можете ознакомится в этой теме.

Если Вам помогли и атмосфера форума Вам понравилась, то заходите к нам чаще! С уважением, Poseidon, Rodman

 
0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема | Центр помощи | Следующая тема »


 




[ Время генерации скрипта: 0.0934 ]   [ Использовано запросов: 22 ]   [ GZIP включён ]


Реклама на сайте     Информационное спонсорство

 
По вопросам размещения рекламы пишите на vladimir(sobaka)vingrad.ru
Отказ от ответственности     Powered by Invision Power Board(R) 1.3 © 2003  IPS, Inc.