x = var('x')    # define a variable x
y = function('y',x)   # define y to be a function of that variable
DE = diff(y, x) - exp(-x^(2))

func = desolve(DE, [y,x],ics=(0,0))

show(func)

plot(func,x )


func = desolve(DE, [y,x],ics=(1,0))

show(func)

plot(func,x )


func = desolve(DE, [y,x],ics=(-1,0))

show(func)

plot(func,x )


func = desolve(DE, [y,x],ics=(2,0))

show(func)

plot(func,x )


func = desolve(DE, [y,x],ics=(-2,0))

show(func)

plot(func,x )

SageDays/Kreyszig-1.3-27 (last edited 2010-08-11 11:16:28 by bhanukiran)