Tomasz Samotyjak, PWSZ writes:
> f[0] = 1/L1 * ( 400*sin(314*t) - f[1] - R1 * f[0] );
f[] shouldn't depend on itself... perhaps the RHS should use y[] instead.
--
Brian Gough
Network Theory Ltd,
Publishing Free Software Manuals --- http://www.network-theory.co.uk/
_
Hi,
I've tried to use gsl fo ode's. The jac and func is:
int
func (double t, const double y[], double f[], void *params)
{
double* mu = (double *)params;
double R1 = mu[0], R2 = R1;
double L1 = mu[1], L2 = L1;
double C1 = mu[2];
f[0] = 1/L1 * ( 400*sin(314*t) - f[1] - R1 * f[0] );
f[