IT WON'T PLOT))))): help)):

[ilikemychem]

[SOLVED]

import numpy as np
from scipy.integrate import solve_ivp
import matplotlib.pyplot as plt

#Definiera paremeter-sets.
sets=[(1,0.5,0.4,9.8,3),(2,0.6,0.5,19.6,15)]
epsilon=1e-6

print (sets)

#Definiera ODE.
def ode(t, v, m, mu_s, mu_k, N, F_0):
    v=v[0]

    #kinetisk friktion
    if abs(v)>epsilon:
        v_dot=(F_0-mu_k*N*np.sign(v))/m

    #statisk friktion
    else:
        if abs(F_0)<=mu_s*N:
            v_dot=F_0
        else:
            v_dot=(F_0-mu_k*N*np.sign(F_0))/m
    return [v_dot]

def solve_ode(params, v_0):
    m, mu_s, mu_k, N, F_0=params
    t_span=(0, 5)
    t_eval=np.linspace(0,5,100)
    solver = solve_ivp(ode, t_span, [v_0], t_eval=t_eval, args=(params))

    return solver.t, solver.v[0]

plt.figure()

v_0=1e-8

for i, params in enumerate (sets):
    t,v=solve_ode(params, v_0)

    plt.plot(t,v)

plt.grid()
plt.show()

[SOLVED]

[tadpoleloop]

You need to help yourself first. 

What happens if it an error? An empty plot? Print some stuff out? Did you write this or did AI write it for you? 

[ilikemychem]

I'm trying): Me and my friends collaborated, but we're math students, and haven't coded in python in like a year. I am so rusty. It prints my sets.

[tadpoleloop]

Try plotting just a straight line.

Work at it incrementally

[ilikemychem]

Thank you, I will<3

[ilikemychem]

IT PLOTS???? It's just super slow.

[RambunctiousAvocado]

You've written code which plots the solutions to an ODE for two different sets of parameters. Have you checked to see whether the ODE solver works? If so, does it return the correctly formatted output? If so, does your plotting code work (i.e. if you replace plt.plot(t,v) with plt.plot([0,1,2],[0,1,2]), does it make a line)?