Интерполяционный полином

[Chardash 307]

Я понял в чем была моя ошибка, не через DSolve решал, а через NDSolve. Но теперь непонятно, почему оно не решает численно задачу? И вручную неверно решил) Буду разбираться.

Clear[f, g]
f[x_, y_] = y;g[x_, y_] = 4 x^3;Solve[{f[x, y] == 0, g[x, y] == 0}]
xxx = {{D[f[x, y], x], D[f[x, y], y]}, {D[g[x, y], x], D[g[x, y], y]}};MatrixForm[xxx]
xxx /. {x -> 0, y -> 0} // Eigenvalues

Вручную быстрее решается) 

Не выдержал, обновился до вин 10ки, поинтереснее пошло дело.

Но пропала сеть на виртуальных машинах... Спасла переустановка

[Chardash 307]

А я, признаться, не очень понимаю о чем вообще речь. Можно полную постановку задачи ?

Все тоже. переделал на 100 раз, не рисует график, шаг 0,00001 - в ответ NDSolve::ndsz: At t == 0.8900537843152546`, step size is effectively zero; singularity or stiff system suspected; At t == 0.5562962428361147`, step size is effectively zero; singularity or stiff system suspected. "Further output of NDSolve will be suppressed during this calculation. И справка на этом же зацикливается, показывает ошибку но решения не показывает, показывает стрелку при нажатии на которую снова возвращаешься к этой ошибке. И график:

 

Система {1 - y, x^2 - y^2}, например, прекрасно считается.

[Fedor 1 613]

"не очень понимаю о чем вообще речь" - человечество открыв мир автоматизированных вычислений снова оказалось в ситуации когда люди могут больше чем понимают и приходится принимать решения в условиях неопределенности. Это нормально. И раньше такое частенько бывало. Просто начался очередной героический этап развития . Такое бывало в математике например, когда изобрели дифференциальное и интегральное исчисление, метод рядов и многое другое ... :)

"сеть на виртуальных машинах" - крутая же молодежь пошла. Скоро и на пенсию пора, окуней ловить :)

[Борман 2 389]

"сеть на виртуальных машинах" - крутая же молодежь пошла.

Крутая, потому что ловит виртуальных окуней.

[Fedor 1 613]

Приятно, что пошли другие инженеры, не те кто только баблосы в экселе считать умеет, да набирать элементы для мкэшек  . А то необходимость увеличения возраста выхода на пенсию министры мотивируют тем, что некому реально работать будет. Не юристики же с экономистиками :)

[Chardash 307]

сеть на виртуальных машинах"

 

это совсем не сложно. на основной системе стоит только самое необходимое, все остальное на виртуальных машинах. Вин10 ка хороша, нужный софт тоже работает нормально. Обновление на про было бесплатным, сменил свою старую 7-ку. Читал о шпионаже 10ки, но есть решения и скрывать тоже особо нечего, много работаю на виртуальных. За молодежь спасибо, хотя во всех этих вопросах это так)

[Fedor 1 613]

Виртуалки отдельные это да, у меня тоже так на переносном винте все, так как фирма хочет быть белой и пушистой для любых контролей :)   Но объединять их еще в сеть, такого не слышал  :) 

Купил ноутбук для дачи, так на него кучу драйверов пришлось ставить. Изначально он для экономии был с досом . Поставил семерку, накормил кучей драйверов, но не уверен, что десятка их имеет тоже...

Наверное стоит поставить десятку на виртуальную машину :)

[Chardash 307]

объединять их еще в сеть

целую инфраструктуру разворачивают и удаленные сотрудники на серверах компании.

 

 

не уверен, что десятка их имеет тоже

скорее всего будет нормально, но ручаться лучше не буду)

[Fedor 1 613]

"целую инфраструктуру разворачивают и удаленные сотрудники на серверах компании" - на удаленных это не проблема обычно. Интересно когда на одном системная виртуализация. Удобно программы отлаживать сетевые...

[Chardash 307]

"целую инфраструктуру разворачивают и удаленные сотрудники на серверах компании" - на удаленных это не проблема обычно. Интересно когда на одном системная виртуализация. Удобно программы отлаживать сетевые...

это вообще не проблема, мост в настройках. у меня он и полетел после установки 10ки, Но восстановил без проблем/

Особенно консоль нравится, похоже как на убунту стала.

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

http://befcee53.ngrok.io - выключил, мешает, кому интересно включу

Не я делал, Илья Никитин с группой, я поднял и допиливаю

[Fedor 1 613]

"сайтострояния" - а ошибки будет мастер заборостроения нынче исправлять ? :)

[Chardash 307]

да, есть, там еще есть работа, просьба сильно не придираться)

[Fedor 1 613]

"x''[t] - 4 x[t]^3 == 0" - вспомнил. Попробуйте в численном решении вместо 4 использовать возмущенную величину, например, 4.00001 или вообще с машинным эпсилон. Иногда в численных делах помогает... :)

При коэф. Пуассона 0.5 в теории упругости возникает деление на 0, а при 0.5 -+ 0.00001 все нормально решается :)

[Chardash 307]

к сожалению, не помогло.   ошибки те же. что то считает, но ошибки вываливает

[Fedor 1 613]

Попробуйте 4 менять на другое число... Судя по аналитике там комплексность есть. Ну и степень попробуйте менять. Краевые условия. Надо найти ближайшее приближение к уравнению и потом посмотреть с какого начинает портится. В каком то интервале задача должна неплохо линейной аппроксимироваться, а там решение экспоненты, но могут и комплексные быть... Сейчас некогда, а в понедельник может будет время поразбираться с причиной. Насколько помню должны условия Липшица соблюдаться...

[Chardash 307]

Федор, спасибо, уже в общем то разобрался как тут все работает, потом еще более детально посмотрю. пока хватит и векторного, вручную сделал

 

Чудо, пока писал тут сообщение, Математика в это время работала в подвисшем состоянии, заглянул потом, что то посчитала.

Ошибки остались, но все равно лучше, может вообще неверно что то сделал. Благодарю еще раз за решение.

На самом деле кроме дифур еще много всего. Надолго.

 

Можете поставить 10ку и откатиться обратно в течении месяца http://www.winline.ru/articles/kak_udalit_windows_10_vernutsya_na_windows_7_8.php

 

еще непонятно, рисует седло, тогда как по корням получается вырожденная.

[execut1oner 11]

раз тема зашла о математике с полиномами как в математике сделать интерполяцию полиномом?

 

написано так -

 

polynomial of degree 6 fit {(1,1),(2.5,8),(3.5,2),(4.5,5),(6,-3),(7.5,4),(8,-3)}

 

но оно не работает , вернее на сайте работает а в программе - нет

[piden 2 874]

написано так -   polynomial of degree 6 fit {(1,1),(2.5,8),(3.5,2),(4.5,5),(6,-3),(7.5,4),(8,-3)}   но оно не работает , вернее на сайте работает а в программе - нет

Да вроде работает

 

[Chardash 307]

или так еще

Expand[interpolatingPolynomial[{{1, 1}, {2.5, 8}, {3.5, 2}, {4.5, 5}, {6, -3}, {7.5, 4}, {8, -3}}, x]]

[execut1oner 11]

а как сделать чтобы я мог задавать степень полинома?

ой не увидел сообщение PIDEN'a