Untitled

261 days ago by jhlee2chn

$x+y=1$

예제 2-6

var('x') # 변수 입력 y = function('y', x) # y = f(x) @interact def _( eq = input_box( 'diff(y,x) - (1-x)/(1+y^3) == 0', 'Diff.eq'), auto_update = False): dq = lambda y : eq # 미분방정식 입력 print dq(y) # 미분방정식 확인 de = desolve ( dq(y), [y, x], contrib_ode=True, show_method=True ) # 일반해 구하기 print de # 미분방정식의 해와 풀이방법 
       

Click to the left again to hide and once more to show the dynamic interactive window

x = var('x') # 변수 입력 y = function('y', x) # y = f(x) @interact def _( eq = input_box('diff(y,x) + 2*x*y == 0', 'Diff.eq'), a1 = slider(-5, 5, 1, 0, 'x of IV'), a2 = slider(-5, 5, 1, 3, 'y of IV'), auto_update = False): dq = lambda y : eq # 미분방정식 입력 de = desolve ( dq(y), [y, x], [a1, a2], contrib_ode=True, show_method=True ) # 특수해(IVP) 구하기 print de # 미분방정식의 해, 풀이방법 show (plot ( de[0], aspect_ratio=1 )) # 특수해-곡선 그리기 
       

Click to the left again to hide and once more to show the dynamic interactive window

x = var('x') # 변수 입력 y = function('y', x) # y = f(x) de = diff(y,x) + 2*x*y == 0 dq= desolve (de, [y, x], [0,3]) # 특수해(IVP) 구하기 print dq plot (dq, (x, -3, 3), aspect_ratio=1 ) # 특수해-곡선 그리기 
       
3*e^(-x^2)
3*e^(-x^2)