-
Notifications
You must be signed in to change notification settings - Fork 0
/
manufactured_solutions.wxm
25 lines (20 loc) · 1001 Bytes
/
manufactured_solutions.wxm
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/
/* [ Created with wxMaxima version 12.04.0 ] */
/* [wxMaxima: input start ] */
f(U,P) := [ diff(U[1],t) + diff(P,x) - diff(U[1],x,2) - diff(U[1],y,2) + U[1]*diff(U[1],x) + U[2]*diff(U[1],y),
diff(U[2],t) + diff(P,y) - diff(U[2],x,2) - diff(U[2],y,2) + U[1]*diff(U[2],x) + U[2]*diff(U[2],y) ];
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
u : [x^2*y*sin(w_*t), -x*y^2*sin(w_*t)]$
p : x*x*cos(w_*t) + y*y*sin(w_*t)$
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
ratsimp(f(u,p))$
trigsimp(%);
div_f : diff(f(u,p)[1],x) + diff(f(u,p)[2],y);
/* [wxMaxima: input end ] */
/* [wxMaxima: input start ] */
at([cos(t*w_)^2*x^3*y^2+(-w_*sin(t*w_)*x^2-2*cos(t*w_))*y+2*cos(t*w_)*x,cos(t*w_)^2*x^2*y^3+w_*sin(t*w_)*x*y^2+2*sin(t*w_)*y+2*cos(t*w_)*x],t=0);
/* [wxMaxima: input end ] */
/* Maxima can't load/batch files which end with a comment! */
"Created with wxMaxima"$