Convert this Matlab code to Python: close all, clear all %Analytical solution f(

Question

Convert this Matlab code to Python: close all, clear all %Analytical solution f(t) = t + 5*exp (-t/2) - 2 f2 = dsolve('Dy = (t-y) /2', y(0) = 3 ', ' t'); f inline (f2); &Differential; equation f1: inline ('t-y) /2'; a = 0 b = 1; %Initial condition y0 = 3; %Number of subintervals N = 100; h = (b - a)/N; x = [a:h:b] '; [yana] = f (x); [yEu] = euler (f1, a y0, b, N); hold on plot (x, yana, 'b'); plot (x, yEu, 'r'); title('Comparison of solutions with N=100'); xlabel('axi t') ylabel ('axi Y'); legend('Analytical Solution', 'Euler'); grid

Explanation / Answer

close(mstring('all')) clear(mstring('all')) #Analytical solution f(t)=t + 5*exp(-t/2) - 2 f = inline(f2) #Differential equation fl = inline(mstring('(t-y)/2')) a = 0 b = 1; print b #Initial condition y0 = 3 #Number of subintervals N = 100 h = (b - a) / N x = mcat([mslice[a:h:b]]).cT [y] = f(x) [yFu] = euler(fl, a, y0, b, N) hold(mstring('on')) mslice[:] y F u mslice[:] title(mstring('Comparison of solutions with N=100')) ylabel(mstring('axi Y')) grid()

Related Questions

Navigate

• Browse (All)
• Browse C
• Subjects

---