Mã:
(*Giai phuong trinh vi phan cap 1*)
Clear[y, x, f, k, T, x0, y0];
x0 = Input["Nhap x0 = ?"];
y0 = Input["Nhap y[x0] = ?"];
y[x0] = y0;
y'[x_] = f[x_, y_] = Input["Nhap y'[x] = ?"];
k = Input["Nhap vao cap cua dao ham"];
T = {y0, y'[x0]};
For[i = 1, i < k, i++,
    	y''[x_] = D[f[x, y], {x, i}];
    	T = Append[T, y''[x0]];
    ];
y = Sum[(T[[i]] (x - x0)^(i - 1))/(i - 1)!, {i, 1, k + 1}];
Print["GIAI PT VI PHAN CAP 1 BANG KHAI TRIEN TAYLO"];
Print["y(", x0, ") = ", y0];
Print["y'[x] = ", y'[x]];
Print["Tinh dao ham den cap k = ", k];
Print["Nghiem y(x) = ", y];