# NAME    :- IBRAR AHMAD 
TUTORIAL:- GROUP 5TH   
ROLL NO :- 11232225    
DATE    :- 26-2-2013   

kill(all)$
y:sin(x)$
is(diff(y,x,3)+2*diff(y,x,2)+diff(y,x)+2*y=0);

kill(all)$
f:2*x^2+y^2-2*x*y+5*y=0$
depends(y,x)$/* Declares functional dependencies among 
variables for purpose of computing derivatives. */
derivabbrev:false$ 
solve(diff(f,x),diff(y,x));

kill(all)$
f:2*x^2+y^2-2*x*y+5*y=0$
depends(y,x)$/* Declares functional dependencies among 
variables for purpose of computing derivatives. */
derivabbrev:true$ 
solve(diff(f,x),diff(y,x));

kill(all)$
de:'diff(y,x,2)+y^2*'diff(y,x)^3=0;
gsoln:ode2(de,y,x);
bc2(gsoln,x=0,y=1,x=1,y=2);
wxdraw2d(implicit(%,x,-5,5,y,-5,5));

kill(all)$
load(implicit_plot)$
de:'diff(y,x,2)+y^2*'diff(y,x)^3=0;
gsoln:ode2(de,y,x);
bc2(gsoln,x=0,y=1,x=1,y=2);
wximplicit_plot(%,[x,-1,5],[y,-3,3],[nticks,1000]);

kill(all)$
points:rk(exp(-t)+u,u,-0.1,[t,2,7,0.01])$
wxplot2d([discrete,points],[t,0,7],[style,[lines,2,4]],
[ylabel," "],[xlabel," t0=2,u0=-0.1,du/dt=exp(-t)+u "])$