1

(2 odpowiedzi, napisanych Mathematica: numeryka i grafika)

Ponieważ rzecz jest bardziej zawiła odpowiadam w notebooku:

http://mth.drac.oa.uj.edu.pl/mufiles/00 … _square.nb

2

(1 odpowiedzi, napisanych Jak to zrobić w Mathematice?)

Funkcja N (Numeric) zwraca liczby przybliżone (przybliżenia dziesietne). Odpowiadają im liczby zmiennoprzecinkowe w Fortranie, Pascalu czy C. Takie liczby zawsze mają kropkę (czyli po polsku przecinek). Liczba 1 i liczba 1. to dwa zupełnie odmienne typy. Np:

In[1]:= 1\[Element] Integers
Out[1]= True

In[2]:= 1.\[Element] Integers
Out[2]= False

In[3]:= 1.\[Element] Reals
Out[3]= True

Kropkę można usunąć poleceniem IntegerPart. Ale należy pamietać, ze każdy wynik eksperymentu (odczyt aparatury pomiarowej) jest liczbą zmiennoprzecinkową, bo taka jest natura pomiaru obarczonego błędem. Brak kropki w takim przypadku jest mylący.

3

(0 odpowiedzi, napisanych Mathematica: numeryka i grafika)

Solve w Mathematica 7.0 nie chce rozwiązać układu, który bez kłopotów rozwiązywała wersja 6.0
(przykładowy notebook: http://mth.drac.oa.uj.edu.pl/mufiles/00 … t_Solve.nb)

4

(1 odpowiedzi, napisanych Jak to zrobić w Mathematice?)

Na początku notebooka trzeba zadeklarować Needs["nazwapakietu`"]. Polecenie Needs opisane jest w helpie. Np. Needs["ComputerArithmetic`"]. Pewien problem może powstać jeśli pakiet nie leży we właściwej kartotece instalacyjnej Mathematic'i. Trzeba go albo najpierw tam włożyć (C:\Program Fles\Wolfram Research\Mathematica\7.0\AddOns\Packages) albo ustawić w odpowiedni sposób kartotekę domyślną (przykłady w http://mth.drac.oa.uj.edu.pl/mufiles/00 … lady/os.nb).

5

(1 odpowiedzi, napisanych Jak to zrobić w Mathematice?)

To rzeczywiście jest jakiś problem wersji 6 i 7. Mathematica 5 poprawnie eksportowała polskie znaki diakrytyczne do html  (konwersji do pdf jeszcze nie było).  W systemie Windows problem można obejść następująco:
1. Trzeba doinstalować drukarkę postscriptową (np HP LaserJet 5P/5MP).
2. Wydrukować notebooka do pliku ps.
3. Plik ps przekonwertować do pdf destillerem Adobe.

Rozwiązanie nie jest dobre, jest co najwyżej dostateczne.
Proszę spróbować podobnego w Linuksie. (Jeśli HP LaserJet 5P/5MP nie pomoże eksperymentować z innymi postscriptowymi). Prosze dać znać czy działa.

6

(4 odpowiedzi, napisanych Jak to zrobić w Mathematice?)

M={1,3,5,7,10,12,13,15};

Select[M,(#<11)&]   (* zachowuje liczby do 10 *)

Select[M, (Divisible[#,3]==False)& ] (* zachowuje liczby niepodzielne przez 3*)

Składnia: Select [lista, operator dyskryminyjacy]
Znak # pełni podobną rolę jak x w zadaniach szkolnych (oznacza niewiadomą).
Znak & kończy definicję operatora.
Operator może być zdefiniowany odrębną dowolnie skomplikowaną procedurą.

7

(3 odpowiedzi, napisanych Mathematica: obliczenia symboliczne)

Zgoda!
Ten sam problem ma funkcja Sin[\[Pi] x]/(x ):

FullSimplify[Sin[\[Pi] x]/(x ),Assumptions->x\[Element]Integers]
daje zero.

8

(4 odpowiedzi, napisanych Jak to zrobić w Mathematice?)

Do sortowania np. w porządku rosnącym służy funkcja Sort. Pogrupować można przy pomocy Split.
Wyrzucić wyrazy też można, ale nie potrafie doradzić jak, bo problem nie jest wystarczająco jasno opisany. Prosze o przykład.

9

(3 odpowiedzi, napisanych Mathematica: obliczenia symboliczne)

Nie widzę sprzeczności w tym, że odwrócenie porządku przejść granicznych prowadzi do odmiennych rezultatów.

Symetrie Lie a całkowalnośc nieliniowych równań różniczkowych zwyczajnych
Wiesław Zajiczek
2 grudnia 2009
(transparencje referatu wygłoszonego na seminarium Zakładu Ogólnej Teorii Wzglednosci)

http://mth.drac.oa.uj.edu.pl/mufiles/00 … ek_Lie.pdf

Mimo wszystko nie rozumiem, jaką role ma do spełnienia całka oznaczona, skoro całkę nieoznaczoną mozna policzyć natychmiast i odjać wartosc w zerze.

f[t_]=Integrate[1/Sqrt[(t+1)^2*(t*(1/3)+1)-t*(t+2)*(2/3)],t]//FullSimplify
f0=Simplify[Limit[f[t],t->0],0<t];
g[z_]=f[z]-f0
Series[g[t],{t,0,3}]//N//Chop

Widać teraz, że problem ulokowany jest w przejściu granicznym do zera z pewną funkcją eliptyczną. Wyskakuje alert, ale kernel liczy. Jak się to potem narysuje wszystko jest  OK.

Relatywistyczne równanie struktury gwiazdy

Notatnik przedstawia rozwiązanie równań Einsteina dla konfiguracji sferycznie symetrycznej, nierotującej,  pozostajacej  w równowadze. Rachunek prowadzony wg. podręcznika Stevena Weinberga: "Gravitation and cosmology" (1972).

http://mth.drac.oa.uj.edu.pl/mufiles/00 … Gwiazdy.nb

13

(0 odpowiedzi, napisanych Jak to zrobić w Mathematice?)

Jak konstruować procedury: trzy proste przykłady.

http://mth.drac.oa.uj.edu.pl/mufiles/00 … ocedury.nb

słowa kluczowe: Module

Proponuję najpierw policzyc całkę nieoznaczoną:

Integrate[1/Sqrt[(t+1)^2*(t*(1/3)+1)-t*(t+2)*(2/3)],t]//FullSimplify

Przenosząc wyrażenia na jedną stronę, każde równanie możemy napisać w postaci "wyrażenie=0" .
Czasami jednak potrzebujemy zachować oryginalny podział na prawą i lewą stronę równania, nie tracąc przy tym możliwości jego przekształcania. Zamieszczony poniżej notatnik podaje przykłady przekształcania równań stronami.

http://mth.drac.oa.uj.edu.pl/mufiles/00 … tronami.nb

słowa kluczowe:
First, Last
List
Equal
/@, Map
@@, Apply

16

(0 odpowiedzi, napisanych Jak to zrobić w Mathematice?)

Przykłady manipulacji na wyrazeniach:
http://mth.drac.oa.uj.edu.pl/mufiles/00 … ssions2.nb

%
Part
Position
Współczynnik (Coefficient)
Licznik i mianownik
Zastąp (ReplaceAll, /.)  i  Wskaż (Rule)  (->)
Wzorce (pattern matching ): podstawy
Clear

17

(0 odpowiedzi, napisanych Jak to zrobić w Mathematice?)

Przykłady manipulacji na wyrazeniach:
http://mth.drac.oa.uj.edu.pl/mufiles/00 … ssions1.nb

Factor/Expand
Together
Collect/FactorTerms
Trygonometryczne
Simplify/FullSimplify
RootReduce
PowerExpand
Refine

18

(2 odpowiedzi, napisanych Mathematica: numeryka i grafika)

Mathematica 6 robi to samo.

Jak widzę problem ma później funkcja Solve. Jest podejrzenie, że Solve traktuje tu jako zespolone zmienne, które w istocie są rzeczywiste. Próbowałem użyc pakietu ReIm (z Mathematica 5.0), ale z mizernym skutkiem.  Sugerowałbym poeksperymentowac z Simplify/FullSimplify z opcjami  Element[x, Reals] ,  Element[y, Reals].